Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	28.8%	104.7 / 10.6	52	G247G (GGT→GGG)	gpmA	2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (32/5);  new base C (0/15);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 3.46e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTT  >  NC_000913/786718‑786994                                                                                                                                          |                                                                                                                                             gTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAcc                                                                                                                                            <  1:1614116/140‑1 (MQ=255)  tATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAccc                                                                                                                                           <  2:243483/140‑1 (MQ=255)         tGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTg                                                                                                                                   >  2:909626/1‑140 (MQ=255)                     ggCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCgctg                                                                                                                       >  1:581694/1‑140 (MQ=255)                      gCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCgctgc                                                                                                                      >  2:498023/1‑140 (MQ=255)                      gCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCgctgc                                                                                                                      >  2:956324/1‑140 (MQ=255)                        cTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCtt                                                                                                                    >  1:1595998/1‑140 (MQ=255)                           aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCg                                                                                                                 >  1:1278405/1‑140 (MQ=255)                           aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCg                                                                                                                 >  2:183397/1‑140 (MQ=255)                                tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGTTTTCGCTGCg                                                                                                            >  2:532934/1‑140 (MQ=255)                                tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                            >  2:268172/1‑140 (MQ=255)                                tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                            >  2:403271/1‑140 (MQ=255)                                tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                            >  2:229816/1‑140 (MQ=255)                                tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                            >  2:1154467/1‑140 (MQ=255)                                tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                            >  1:1158995/1‑140 (MQ=255)                                  gcgcAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTCCTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGAt                                                                                                          <  2:611683/140‑1 (MQ=255)                                  gcgcAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGAt                                                                                                          <  1:105632/140‑1 (MQ=255)                                     cAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                       >  1:74660/1‑140 (MQ=255)                                       gCGTCGCATCAGGCAATGTGCTCCATTTTTAGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTCCCCTG‑TTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                     <  2:159364/140‑1 (MQ=255)                                       gCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGTTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                     <  1:147128/140‑1 (MQ=255)                                             cATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt                                                                                               >  1:1058001/1‑140 (MQ=255)                                                cAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                            <  2:581694/140‑1 (MQ=255)                                                 aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTG‑TTTTGCAACCGCCGCTGCTTCCGCTGCGATCTCGTCAGCATTAcc                                                                                           <  1:275447/140‑1 (MQ=255)                                                       tgtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                     >  1:479899/1‑140 (MQ=255)                                                        gtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataa                                                                                    >  1:332531/1‑140 (MQ=255)                                                           cTCCTTTGTTAGAAAAAAAAAAGCCGACTCCTTTGCAGTCGGTTTTCTCTTTTTAAACGAATGACTTTTACTTCCCTTTCCCCTG‑TTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAg                                                                                 <  2:243977/140‑1 (MQ=255)                                                           cTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAg                                                                                 <  1:183397/140‑1 (MQ=255)                                                            tCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                >  1:432052/1‑140 (MQ=255)                                                            tCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                >  2:589353/1‑140 (MQ=255)                                                             ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                               >  2:50449/1‑140 (MQ=255)                                                             ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                               >  2:52021/1‑140 (MQ=255)                                                               aTTGTTAGCAAAAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGtt                                                                             <  1:488671/140‑1 (MQ=255)                                                                             aaaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaaa                                                               >  2:786275/1‑140 (MQ=255)                                                                                    gACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                        <  2:962205/140‑1 (MQ=255)                                                                                     aCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                       >  1:698830/1‑140 (MQ=255)                                                                                      cTCCCTTGCAGTCGGTTTTCTCATTTAAAACGAATGACTTTTCCTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGTGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                      <  1:123235/140‑1 (MQ=255)                                                                                      cTCACTTGCAGTCGGCTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                      <  1:551133/140‑1 (MQ=255)                                                                                       tcccTTGCAGTCGGTTTTCTCATTTTAAACGAATGCCTTTTCCTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGa                                                     <  2:74660/137‑1 (MQ=255)                                                                                            tGCAGTCGGTTTTCTCATTTAAAACGAATGACTTTTCCTTCGCTTTCCCCTG‑TTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCa                                                <  1:403271/140‑1 (MQ=255)                                                                                            tGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTG‑TTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCa                                                <  2:933900/140‑1 (MQ=255)                                                                                                gTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATaca                                            >  1:348176/1‑140 (MQ=255)                                                                                                gTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATaca                                            >  1:1466984/1‑140 (MQ=255)                                                                                                 tCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                           >  2:1235294/1‑140 (MQ=255)                                                                                                 tCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                           >  1:337562/1‑140 (MQ=255)                                                                                                 tCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACCCATacac                                           >  2:1000819/1‑140 (MQ=255)                                                                                                   ggCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCa                                         >  1:600782/1‑140 (MQ=255)                                                                                                       ttCTCTTTTTAAACAAATGACTTTTACTTCCCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                     <  1:564364/140‑1 (MQ=255)                                                                                                        tctcATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGc                                    >  1:121035/1‑140 (MQ=255)                                                                                                             tttaAACCGAATGACTTTTCTTTCGCTTTACCCTG‑TTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                               <  1:1537684/140‑1 (MQ=255)                                                                                                             tttaAAACGAATGACTTTTACTTCGCTTTCCCTTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                               <  1:229816/140‑1 (MQ=255)                                                                                                               ttAAACGAATGACTTTTACTTCGTTTTACCTTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                             <  1:1154467/140‑1 (MQ=255)                                                                                                               ttAAACGAATGACTTTTACTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                             <  1:1220667/140‑1 (MQ=255)                                                                                                                aaaaCGAATGACTTTTACTTCGCTTTCCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                            <  1:1607899/139‑1 (MQ=255)                                                                                                                  aaCGAATGACGTTTACTTCGCTTTCCCCTG‑TTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                          <  1:585359/140‑1 (MQ=255)                                                                                                                         gACGTTTACTTCGCTTTCCCCTG‑TTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                   <  1:1345016/140‑1 (MQ=255)                                                                                                                         gACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                   <  2:1058001/140‑1 (MQ=255)                                                                                                                            ttttCCTCCGTTTTCCCTTG‑TTTTCAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                <  2:600782/139‑1 (MQ=255)                                                                                                                            ttttACTTCGCTTCCCCCTG‑TTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                <  2:432052/139‑1 (MQ=255)                                                                                                                            gTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                <  1:956324/140‑1 (MQ=255)                                                                                                                               tACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc             <  2:121035/140‑1 (MQ=255)                                                                                                                                ccttcccttccCCCTG‑CTTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa            <  2:337562/139‑1 (MQ=255)                                                                                                                                ccTTTGCTTT‑CCCTGCGTTTTCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTCCAGCGGTTTGAAATTCTCGTCGAACTCATACACCACCGGCACGCCATGCGGGATATTAAGCTCa            <  2:1278405/139‑1 (MQ=255)                                                                                                                                      cTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc      >  1:1246468/1‑140 (MQ=255)                                                                                                                                      cTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc      >  1:109012/1‑140 (MQ=255)                                                                                                                                          aCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt  >  1:203779/1‑140 (MQ=255)                                                                                                                                          aCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt  >  2:842352/1‑140 (MQ=255)                                                                                                                                          aCCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt  <  2:332531/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                             GTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTG‑GTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTT  >  NC_000913/786718‑786994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑