Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,871	0	T	A	43.9%	83.9 / 43.9	67	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (34/3);  new base A (0/29);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.11e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTT  >  NC_000913/1755724‑1756015                                                                                                                                                      |                                                                                                                                                  cATCGTACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTaaa                                                                                                                                                <  1:524212/151‑1 (MQ=255) cATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTaa                                                                                                                                                 >  2:567045/1‑150 (MQ=255)         cGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑Ccgccgc                                                                                                                                       >  2:454579/1‑151 (MQ=255)          gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCt                                                                                                                                      >  1:322050/1‑151 (MQ=255)          gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCt                                                                                                                                      >  1:456341/1‑151 (MQ=255)           aaaaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCt                                                                                                                                      >  2:240540/1‑150 (MQ=255)                        aaGAGTTGTTAGCAAAATGGCGGGCCGTGGGAATGACCTTTATGCGTCGAACCTTCTCTCATGGTAACTATGAAAAACCCGGC‑CAGCGATTTAAAATTCGGCGCACCGTAATGAGAAAAA‑CTGGAAAAA‑CCGCCGCTATCCTGCTTGAt                                                                                                                          <  2:28253/149‑1 (MQ=255)                        aaGAGATTTTAGAAAAAATGGTGGCCGCGGGCAGAACCGTTATGCGTCTAACCTTCCCTAAGGGGAACTATGCAAAACAGGGA‑CAGCGCATACAAAATCGGCGAAAGGTGATGAGAAAAC‑CGGAAAAAA‑CCGCCGCTATCCTGCTTGAt                                                                                                                          <  1:111964/149‑1 (MQ=255)                            gaTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCaa                                                                                                                     >  2:429562/1‑150 (MQ=255)                              tGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                  >  2:76296/1‑151 (MQ=255)                              tGTTAGCAAAAATGCGGGCCGCGGGCATGAACGTTAGGCTTCTGACCTTCTCTAATGGTGACTATGAAAAACGCGGC‑CAGCGCATTCAGAATCTGCCCACGGTGATGAGAAAAA‑CTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAg                                                                                                                   <  2:349220/150‑1 (MQ=255)                              tGTTACCAAAAATGCGGCGCGCTGTCAGAAGCGTTTTGCTTTAAAATTTCTCTTGTGGTAATTATGCAAAACGCGTT‑CACCCCATTAAGATTTGGCGCACCGTGATGAAAAAAC‑CTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCaaa                                                                                                                    <  1:570242/149‑1 (MQ=255)                               gTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGt                                                                                                                 >  2:90693/1‑151 (MQ=255)                                taagaaaaaaGGGTGGACGCGGGAATGAACGTTAGGCTTCTAACCTTCTCTCAGGGTAATTAGGAAGAACACGGT‑CACCGCATTAAAAATCTGCGAACCGTAAGGAAAAAAA‑CGGAAAAAC‑CCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                  <  1:115838/143‑1 (MQ=255)                                 tactaaaatTGCTGGCCGCTGCCAGGACCGTTTTGCTTCTAACCTTCCCTCATGGTAACTAGGCAAACCCCGGT‑CAGCGCATAAAAATTCTCCCCACGGTGATGACAAAAA‑TTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                                <  2:853309/145‑1 (MQ=255)                                     aaaaaGGCGGGACGCGGACATAAATGTTATGCTTTTAATCTTCTCTAGGGGTAACTATCCAGAACACGCT‑CCGCCCTTTAAGATTCTGCGCACGGTGAGGAAAAAAA‑CGGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                                <  1:891112/145‑1 (MQ=255)                                       aaaTGCTGCACGCTGATATGAGCGTTATGTGTCTGATCTTTTCTGATGGTTATTATGAAAAACGCGGT‑CGGCGCATCGAAAACCCGCG‑AATTGGAAGAAAAAAATGTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                                <  2:1033444/144‑1 (MQ=255)                                           gCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGc                                                                                                     >  1:817092/1‑151 (MQ=255)                                            cTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGc                                                                                                     >  2:15410/1‑150 (MQ=255)                                             tGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAc                                                                                                   >  1:964142/1‑151 (MQ=255)                                              ggACGCTGGCATAAACTTTATGCGTCTAAACTTCTCTATGGGTGACTATGAAAAACCCGGC‑CAGCCCATTAAAATTCTGCGAAACGTGAGGAAAAAAA‑CGGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAc                                                                                                   <  1:838348/150‑1 (MQ=255)                                               gACGCTGGCTTGAACGTTTGGCTTTTGAATTTCTCTAGTGGTGACTATGCAGAACAGGGT‑CGGCGCTTTCAAACTCGGCGCACCGTGATGAGAAAAA‑CTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCa                                                                                                 <  2:222174/151‑1 (MQ=255)                                                      gCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACt                                                                                          >  2:907918/1‑151 (MQ=255)                                                            aCGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAgg                                                                                    >  2:562114/1‑151 (MQ=255)                                                                tATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                >  1:190069/1‑151 (MQ=255)                                                                 aTGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                >  2:780866/1‑150 (MQ=255)                                                                      tttaatcTTCTCTCATGGTGATTATGAAGAACACGGC‑CAGCGCATGCAGATTCTGCGCAACGTGAGGAGAAAAA‑CTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTaa                                                                              <  2:1096912/141‑1 (MQ=255)                                                                       cTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgacg                                                                         >  2:513066/1‑151 (MQ=255)                                                                          acctcctctcATGGTGACTAGGAAGACCGCGGC‑CAGCCCATACAAAATCTCCGAAACGTGATAAAAAAAA‑CGGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGtt                                                                       <  1:64586/145‑1 (MQ=255)                                                                          accTTCTCTCAGGGTGACTAGCAAGAACACGGT‑CAGCCCATTAAGAATCTGCCAACCGGAAGGAAAAAAA‑GTAGTAAAC‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGttt                                                                      <  1:727734/149‑1 (MQ=255)                                                                          aaCTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGttt                                                                      >  2:53869/1‑151 (MQ=255)                                                                          aaCTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGttt                                                                      >  1:1083878/1‑151 (MQ=255)                                                                             ttCCCTCAGGGTGACTATGCAAAACACGGT‑CAGCGCATTAAAAATCGGCGCAACGTAAGGAAAAAAA‑CTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGt                                                                        <  2:892320/146‑1 (MQ=255)                                                                              tctctcATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                  >  1:508207/1‑151 (MQ=255)                                                                              tctctcATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                  >  1:914614/1‑151 (MQ=255)                                                                              tctctcATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                  >  1:63289/1‑151 (MQ=255)                                                                               ctctcGTGGTGAATAGGAAGAGC‑CGGTCCAGCGCATTCAGAATCTGCGCAGGGTGATGAGAAAAA‑ATTGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGtt                                                                       <  1:657647/145‑1 (MQ=255)                                                                                gctcATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATAAGCAAAA‑CTGGGAAAA‑CCGCCGCTATCCTGCTTGATACGAAAGGTACGGAAATCCGCACCATGGAACTGGAAGGCGGTAACGACGTTTCTCTg                                                                 >  1:197775/2‑150 (MQ=255)                                                                                         gattaggCAAAACACGGC‑CAGCGATTTAAAAACCGGCGCAACGTAAGGAAAAAAA‑TTGTAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTc                                                       <  2:579836/145‑1 (MQ=255)                                                                                            tATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTc                                                       >  2:1075069/1‑148 (MQ=255)                                                                                             aTGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAc                                                   >  2:996407/1‑151 (MQ=255)                                                                                              tGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAcc                                                  >  1:1070123/1‑151 (MQ=255)                                                                                                 agaccccGGC‑CAGCGTATCAAGAACCTCCGAAAGGTGAGGAAAAAAA‑GTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                 <  1:762051/143‑1 (MQ=255)                                                                                                 aaaaCCCGGT‑CAGCGAATTAAAAATCTGCGCAACGTGAGGAAAAAAA‑CTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                 <  2:1096241/147‑1 (MQ=255)                                                                                                 aGAACCCGGT‑CAGCGCATTCAGAATCTGCGCAGGGTGATGAAAAAAA‑CTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCtt                                                <  2:157542/150‑1 (MQ=255)                                                                                                  gAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                               >  1:825933/1‑150 (MQ=255)                                                                                                     cacGGT‑CAGCGCATTAAAATCCTGCGAAACGTGATGAAAAAAA‑CTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTa                                              <  2:309686/148‑1 (MQ=255)                                                                                                     cacGGC‑CAGCGCATTCAGAACCTGCGCAAGGTGATGAGAAAAA‑CTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                               <  1:610903/147‑1 (MQ=255)                                                                                                            cAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcacc                                      >  2:288863/1‑150 (MQ=255)                                                                                                              gcgaATTCAAAATCTGCGAAACGTGATGAAAAAAA‑CTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                    <  2:217284/150‑1 (MQ=255)                                                                                                               cgcATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTg                                  >  2:763116/1‑151 (MQ=255)                                                                                                                 cATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                >  1:470283/1‑151 (MQ=255)                                                                                                                    tCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATaa                              >  1:104638/1‑150 (MQ=255)                                                                                                                      aGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                           >  2:305869/1‑151 (MQ=255)                                                                                                                      aGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                           >  1:794354/1‑151 (MQ=255)                                                                                                                      aGAATCTGCGCAACGGGATGAAAAAAA‑CTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATaaa                             <  2:182555/149‑1 (MQ=255)                                                                                                                      aGAATCTGCGAAAGGTGATGAAAAAAT‑ATGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATaa                              <  2:508231/148‑1 (MQ=255)                                                                                                                              cgcacggTGATGAAAAAAA‑CTGGAAAAG‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTa                      <  2:817133/142‑1 (MQ=255)                                                                                                                               gcAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCgg                  >  1:198428/1‑151 (MQ=255)                                                                                                                               ccaaccgtaataaacaaa‑aTGGTAAAAGCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCg                   <  2:9119/136‑1 (MQ=255)                                                                                                                                 aaCGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                 >  1:464528/1‑150 (MQ=255)                                                                                                                                 aaCGTGAGGAAAAAAA‑CTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCa                <  2:282766/151‑1 (MQ=255)                                                                                                                                          agaaaaa‑CTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGaa        <  1:35395/147‑1 (MQ=255)                                                                                                                                          aGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGaa        >  2:311615/1‑150 (MQ=255)                                                                                                                                            aaaaa‑CTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCaa               <  1:624070/140‑1 (MQ=255)                                                                                                                                               aa‑CTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt  <  1:269185/151‑1 (MQ=255)                                                                                                                                                aagtGGAAAAG‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt  <  1:531532/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAA‑CTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTT  >  NC_000913/1755724‑1756015

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