Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	44.8%	113.4 / 21.0	68	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (34/3);  new base C (0/30);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 5.46e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATT  >  NC_000913/4537960‑4538252                                                                                                                                                            |                                                                                                                                                 aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CAc                                                                                                                                                 >  1:506244/1‑151 (MQ=255) aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CAc                                                                                                                                                 >  2:307036/1‑151 (MQ=255)     ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑ttt                                                                                                                                            <  2:947724/151‑1 (MQ=255)      cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑ttt                                                                                                                                            <  2:12310/150‑1 (MQ=255)         aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTc                                                                                                                                        >  2:1061729/1‑151 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTc                                                                                                                                        >  1:359768/1‑150 (MQ=255)                      aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCAc                                                                                                                           <  1:14303/151‑1 (MQ=255)                      aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCAc                                                                                                                           >  2:822938/1‑151 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCAcc                                                                                                                          >  2:595495/1‑150 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCAcc                                                                                                                          >  2:372183/1‑150 (MQ=255)                         ccATCTGGTGATGAGACGGGAGCAAGCTTTTCCCATTTTAAATTATTCCCGGTAAAATCAAGTTA‑AAATACGGCATCCTTCCGCAACCCTGTTTTGG‑CTT‑CGCCATAAATAAG‑CCAGT‑TTTTAC‑CCCC‑TTTATCCACAACCGCCGCACCa                                                                                                                         <  2:995433/150‑1 (MQ=255)                           aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGcc                                                                                                                      >  1:571122/1‑151 (MQ=255)                                     gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                            >  2:47607/1‑151 (MQ=255)                                       gaCGGAAGCAGGCTTATTCCATTTTAAATTATTCCCGGTAAAACCAAGTTC‑AAATACGGCATCCTTTCGCAATCCGGGTTTGC‑CTT‑CGCCATTAATAAC‑CCTGT‑TTTTAT‑CCCC‑TTTATCCACACCCGCCGCACCAGCCGTTCCGTACCa                                                                                                           <  2:263578/150‑1 (MQ=255)                                          gggACCAAGTTTATTCATTTTTAATTTTTTCCCGGTAAAATCAATTTA‑AAAACCGCCACCCGTTCGCACCCCTTTTTTGG‑CTC‑CGCCATAAAAAAG‑CCAGG‑TTTTAC‑CCCT‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                        <  2:871960/150‑1 (MQ=255)                                              gCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                   >  1:273443/1‑151 (MQ=255)                                               cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                  >  2:1043937/1‑151 (MQ=255)                                                   cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                >  2:481137/1‑149 (MQ=255)                                                     ttttCCATTTTAAATTATTCCCGGTAAAATCAATTTA‑AAATACGGCATCCTTTCGCAACCCTGTTTTGG‑TTT‑CGCAATTAAAAAG‑CCAGT‑TTTTAT‑CCCT‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg                                                                                             <  1:999124/148‑1 (MQ=255)                                                      aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                           >  1:130110/1‑151 (MQ=255)                                                       ttCCTTTTAAAATTATTCCCGGTAAAATC‑AGTTCAAAATCCGGCCCCCGTCCGC‑ATTCTTGTTTGG‑TTTCCCCCATTAAAACC‑CCGGG‑TTTTAC‑CCCC‑TTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                           <  1:1007402/150‑1 (MQ=255)                                                        tCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                          >  1:771082/1‑150 (MQ=255)                                                            tttaAAATTATTACGGGTGAAATCAGGTTA‑AAATACGCCATCCGTTCCAAATCCTGGTTTGG‑CTC‑CGCCATTAATAGG‑CCAGG‑TTTTAT‑CCCC‑TTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa                                                                                        <  1:136246/148‑1 (MQ=255)                                                              ttAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                    >  2:497326/1‑150 (MQ=255)                                                              ttAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                    >  1:1040447/1‑150 (MQ=255)                                                               tAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                  >  1:294692/1‑151 (MQ=255)                                                                aaaTTTTTCCCGGTAAAATCAATTTA‑AAATCCGCCATCGGTCCGCATTCCTTTTTTGC‑TTT‑CGCCATAAATAAC‑CCGGG‑TTTTAC‑CCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa                                                                                   <  2:249249/149‑1 (MQ=255)                                                                   ttattcCGGGTAAAATAAAGTTA‑AAAAACGCCTCCCTTTCGAAACCCTGTTTGGG‑CTC‑CGCAATTAATAAG‑CCAGT‑TTTTAC‑CCCC‑TTTACCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                               <  2:579974/150‑1 (MQ=255)                                                                   ttattcCCGGTAAAACCAATTTA‑AAATACGGCTTCCTTCCGCAA‑CCTTGTTTGGCCTT‑CGCCATTAAAAAG‑CCAGT‑TTTTAC‑CCCT‑TTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                               <  2:730839/150‑1 (MQ=255)                                                                    tattaCCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                             >  1:253297/1‑151 (MQ=255)                                                                    tattaCCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                             >  2:117674/1‑151 (MQ=255)                                                                      ttcccGGTGAAATAAATTTA‑AAATCCGCCTTCCGTTCGCAATCCTGTTTGGG‑CTC‑CCCCATTAAAAAG‑CCAGT‑TTTTAT‑CCCT‑TTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                             <  1:807487/146‑1 (MQ=255)                                                                               aaaTCAATTTA‑AAATCCGGCATCCGTTCCCAACCCTGTTTTGG‑CTT‑CCCCATTAAAAAC‑CCGGT‑TTTTAT‑CCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                    <  1:339597/149‑1 (MQ=255)                                                                               aaaTCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                  >  1:991638/1‑151 (MQ=255)                                                                                       ttC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                          >  1:726688/1‑151 (MQ=255)                                                                                         c‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                         >  1:534401/1‑150 (MQ=255)                                                                                           aaagccGGCATCCGTTCCAAACCCGGGTTTGG‑CTG‑CG‑CAATAAAAAA‑CCCGG‑TTTTAA‑CCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                          <  2:797962/143‑1 (MQ=255)                                                                                           aaaTAGGCCTTCCTTTCGCAATCCGGGTTGGC‑CTC‑CGCAATAAAAAGG‑CCGGT‑TTTTAT‑CCCC‑TTAATTCACATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                         <  1:34267/149‑1 (MQ=255)                                                                                               aCGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                    >  2:1011016/1‑150 (MQ=255)                                                                                                cGGCATCCTTCCGCATTCCTGGTTTGG‑CTC‑CCCCATTAAAAAG‑CCAGT‑TTTTAT‑CCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                    <  1:261123/149‑1 (MQ=255)                                                                                                  gCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                >  1:340928/1‑151 (MQ=255)                                                                                                      ccGTCCGCATTCTTGTTTTGG‑TTT‑CCCCATAAAAAAC‑CCGGT‑TTTAAT‑CCCT‑TTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                              <  2:970913/149‑1 (MQ=255)                                                                                                        gTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                           >  1:262034/1‑150 (MQ=255)                                                                                                            gCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                       >  2:457285/1‑150 (MQ=255)                                                                                                            gCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                       >  1:821912/1‑150 (MQ=255)                                                                                                              aaTCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                    >  1:844362/1‑151 (MQ=255)                                                                                                                     gTTTGG‑CTTCCCCCATAAATAGG‑CCGGT‑TTTTTT‑CCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                 <  2:48309/148‑1 (MQ=255)                                                                                                                     gTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                             >  2:498435/1‑151 (MQ=255)                                                                                                                      tttGG‑CTC‑CCCAATAAAAAAG‑CCGGG‑TTTTAT‑CCCC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                  <  1:680449/145‑1 (MQ=255)                                                                                                                      ttggg‑CTC‑CCCCATAAATAAC‑CCAGT‑TTTTAC‑CCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                <  2:616039/144‑1 (MQ=255)                                                                                                                         gg‑CTT‑CCCCATTAAAACG‑CCAGT‑TTTTCT‑CCCC‑TTTATCAAAAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                           <  2:40546/149‑1 (MQ=255)                                                                                                                               tccccAATAAAAAGCCCAGT‑TTTTAC‑CCCT‑TTTATACAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                          <  2:968172/144‑1 (MQ=255)                                                                                                                                ccccATTAAAAAG‑CCAGT‑TTTTCT‑CCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                    <  2:827376/149‑1 (MQ=255)                                                                                                                                cGCCATAAATAAC‑CCGGT‑TTTTAC‑CCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                    <  2:80118/151‑1 (MQ=255)                                                                                                                                cGCAATAAAAAAC‑CCAGT‑TTTTAC‑CCCT‑TTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                    <  1:830172/151‑1 (MQ=255)                                                                                                                                   cATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGTGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                 >  1:645487/1‑151 (MQ=255)                                                                                                                                     ttAATAAG‑CCGGT‑TTTTAT‑CCCC‑TTTATCAACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                     <  2:530655/145‑1 (MQ=255)                                                                                                                                     ttAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                >  1:9033/1‑150 (MQ=255)                                                                                                                                      taataaG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa               >  1:720629/1‑150 (MQ=255)                                                                                                                                      taaaaaG‑CCCGGTTTTTAC‑CCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                  <  2:363145/148‑1 (MQ=255)                                                                                                                                       aataaG‑CCAGT‑TTTTAT‑CCCT‑TTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                 <  1:494926/147‑1 (MQ=255)                                                                                                                                       aataaC‑CCGGT‑TTTTGC‑CTCC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                          <  1:789787/138‑1 (MQ=255)                                                                                                                                        ataaG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt            >  2:815848/1‑151 (MQ=255)                                                                                                                                        ataaC‑CCAGG‑TTTAGT‑CCCC‑TTCATACACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt            <  2:490533/151‑1 (MQ=255)                                                                                                                                              ccAGG‑TTTTAC‑CCCT‑TTTTGTCGCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt         <  1:247167/149‑1 (MQ=255)                                                                                                                                               ccag‑TTTTATCCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg        <  2:242267/147‑1 (MQ=255)                                                                                                                                                aGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc     >  1:18628/1‑151 (MQ=255)                                                                                                                                                     tttAT‑CCCT‑TTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt  <  1:70433/150‑1 (MQ=255)                                                                                                                                                            |                                                                                                                                                 ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTC‑AAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATT  >  NC_000913/4537960‑4538252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑