Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	57.7%	99.3 / 121.7	78	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (28/5);  new base C (0/45);  total (28/50)
Fisher's exact test for biased strand distribution p-value = 1.94e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-16

CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786706‑787005                                                                                                                                                      |                                                                                                                                                       cGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAc                                                                                                                                                       >  1:1074448/1‑151 (MQ=255)     agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                                   >  2:1341446/1‑151 (MQ=255)        gTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGtt                                                                                                                                                <  1:467116/151‑1 (MQ=255)          ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTg                                                                                                                                              <  2:929359/151‑1 (MQ=255)          ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTg                                                                                                                                              <  2:271316/151‑1 (MQ=255)              tATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAc                                                                                                                                          >  2:1148276/1‑151 (MQ=255)               aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAcc                                                                                                                                         >  2:1158891/1‑151 (MQ=255)                tCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccg                                                                                                                                        >  1:1239941/1‑151 (MQ=255)                     tGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCgct                                                                                                                                   >  1:1348070/1‑151 (MQ=255)                              acaGGCCTGATAAGCCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTTTTAACAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTAAAACGAATGACTTTTTCTTCTCTTTCCCCTGGTTGGAAACCGCCGCTGCTTTCGCt                                                                                                                          <  2:865787/151‑1 (MQ=255)                                 ggCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTTTTAGCAAAAAAAAACCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                                       <  1:267773/151‑1 (MQ=255)                                 ggCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCACTGCg                                                                                                                       >  2:908850/1‑151 (MQ=255)                                  cctttaTAAGTCGCCCGCGGTCGCATCAGGAAATGTCCTCCATTGTTAGCAAAAAAAAAGCCATCTCACTTGCTGTGGCCTTTCTCATTTTAAACGAATGACTTTTTCTTCGTTTTCCCCTGTTTTGAACCCGCCGCTGCTTTCGCTGCGa                                                                                                                      <  2:1298059/146‑1 (MQ=255)                                         aaGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                               >  1:154810/1‑151 (MQ=255)                                           gTCGCGCAGCGTCGAATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGCCTCACTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAg                                                                                                             <  1:705917/151‑1 (MQ=255)                                                gcAGCGTCGCATCAGGCAATGTGCTCCATTTTTAGCAACAAAAAAGCCGCCTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTa                                                                                                        <  2:595166/151‑1 (MQ=255)                                                   gCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGCCTCACTTGCGGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     <  2:1265866/151‑1 (MQ=255)                                                    cGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     >  1:273981/1‑150 (MQ=255)                                                     gTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     >  1:70696/1‑149 (MQ=255)                                                     gTCGCATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAg                                                                                                   <  1:1242138/151‑1 (MQ=255)                                                       cGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGat                                                                                                 <  1:1222113/151‑1 (MQ=255)                                                       cGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     >  1:965487/1‑147 (MQ=255)                                                       cGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     >  1:572062/1‑147 (MQ=255)                                                          aTCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              <  1:1056683/151‑1 (MQ=255)                                                          aTCAGACAATGGCCTCCATTGTTAGAAACAAAAAACCCGCCTCACTTGCAGTCGTCTTTCTAATTTTAAACGAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              <  2:519025/151‑1 (MQ=255)                                                            cAGGCATTGTGCTCCATTGTTAGCAAAAAAAAAGCCGCCTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAg                                                                                            <  1:1323231/151‑1 (MQ=255)                                                            cAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              >  1:773417/1‑149 (MQ=255)                                                            cAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              >  1:1228905/1‑149 (MQ=255)                                                              ggCAATGTGCTCCATTTTTAACAAAAAAAAAGCCGCCTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGCCTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                          <  1:222721/151‑1 (MQ=255)                                                              ggCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                          <  1:59243/151‑1 (MQ=255)                                                               ccattgtgCTCCATTGTTACCAACAAAAAAGCCGACCCCCTTCCATTCGGCTTTCTCTTTTTAAACGAATGACTTTTCCTTCGCTTTCCCCTGGTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                         <  2:637837/147‑1 (MQ=255)                                                                    gtgCTCCATTGTTAGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCa                                                                                     <  2:184885/150‑1 (MQ=255)                                                                      ccTCCATTTTTAGCAAAAAAAAAGCCGCCTCACTTGCAGTCGTTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg                                                                                  <  2:622512/150‑1 (MQ=255)                                                                         ccATTGTTACAAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACAAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt                                                                               <  2:346692/151‑1 (MQ=255)                                                                           aTTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                             <  1:725597/151‑1 (MQ=255)                                                                             tGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGa                                                                            <  1:1060013/150‑1 (MQ=255)                                                                              gTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaaa                                                                          <  1:509024/151‑1 (MQ=255)                                                                                tAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAtt                                                                        >  1:713826/1‑151 (MQ=255)                                                                                     aaaaaaaaGCCGCCTCACTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                   <  2:208920/149‑1 (MQ=255)                                                                                      aaaaaaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                  <  1:225844/150‑1 (MQ=255)                                                                                       aaaaaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                                 >  2:1026395/1‑151 (MQ=255)                                                                                        aaaaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                                 >  1:1121266/1‑150 (MQ=255)                                                                                        aaaaaGCCGACTCACTTGCAGGCGGCTTTCTCTTTTTAACCGAATGTCTTTTACTTCGTTTTCCCCTGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGa                                                                <  2:442894/151‑1 (MQ=255)                                                                                             gCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTc                                                            >  2:356984/1‑150 (MQ=255)                                                                                              ccGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                          <  1:11317/151‑1 (MQ=255)                                                                                                gccTCACTTGCGGTCGTTTTTCTCATTTAAAACGAATGACGTTTCCTTCTCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATac                                                        <  2:640861/149‑1 (MQ=255)                                                                                                gACTCTCTTGCAGTCGGTTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATac                                                        <  2:511357/151‑1 (MQ=255)                                                                                                   tCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     >  1:1057736/1‑151 (MQ=255)                                                                                                    cACTTGCAGTCGGCTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     <  1:1075563/150‑1 (MQ=255)                                                                                                     aCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     >  1:903405/1‑149 (MQ=255)                                                                                                        tGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGc                                                  <  1:947036/149‑1 (MQ=255)                                                                                                         gcggTCGGCTTTCTCATTTTAAACGAATGACTTTTCCTTCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                <  1:207872/147‑1 (MQ=255)                                                                                                         gCAGTCGGTTTTCTCATTTAAAACGAATGACTTTTCCTCCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                <  1:619006/150‑1 (MQ=255)                                                                                                            gTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACg                                            <  1:1334609/151‑1 (MQ=255)                                                                                                              cGGTTTTCTCATTTTAAACGAATGACGTTTACTTCTCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                          <  2:1311503/151‑1 (MQ=255)                                                                                                                gCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                                        <  2:1239669/151‑1 (MQ=255)                                                                                                                  tttCTCATTTTAAACGAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                       <  1:335154/150‑1 (MQ=255)                                                                                                                    tctcTTTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                     <  1:825124/150‑1 (MQ=255)                                                                                                                    tctcATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                     <  1:1232144/150‑1 (MQ=255)                                                                                                                    tctcATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                     <  1:1082640/150‑1 (MQ=255)                                                                                                                     ctcATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                                   >  1:1349946/1‑151 (MQ=255)                                                                                                                       cATTTTAAACGAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                                    <  1:1001551/148‑1 (MQ=255)                                                                                                                       cATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                                 <  1:971268/151‑1 (MQ=255)                                                                                                                            aaaaCGAATGACGTTTACTTCTCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTaa                            <  2:161090/150‑1 (MQ=255)                                                                                                                             aaaCGAATGACGTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                           <  1:544421/151‑1 (MQ=255)                                                                                                                              aaCGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                          >  2:1258187/1‑151 (MQ=255)                                                                                                                               aCGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCt                         >  1:1231184/1‑151 (MQ=255)                                                                                                                                  aaTGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                       >  1:1244146/1‑150 (MQ=255)                                                                                                                                  aaTGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                       >  1:420761/1‑150 (MQ=255)                                                                                                                                     gACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                    <  1:113741/150‑1 (MQ=255)                                                                                                                                     gACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                    >  1:819119/1‑150 (MQ=255)                                                                                                                                     gACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                    >  1:556522/1‑150 (MQ=255)                                                                                                                                        gTTTCCTCCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc                 <  2:1300074/150‑1 (MQ=255)                                                                                                                                         tttACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct                <  1:142100/150‑1 (MQ=255)                                                                                                                                            ccTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttc            <  1:709091/150‑1 (MQ=255)                                                                                                                                                cGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCg        >  2:89150/1‑151 (MQ=255)                                                                                                                                                  ttttCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc       <  1:401965/149‑1 (MQ=255)                                                                                                                                                   tttCCCCTGTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc       <  2:572040/149‑1 (MQ=255)                                                                                                                                                      aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  2:244862/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786706‑787005

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