Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	43.0%	161.4 / 18.7	95	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (41/12);  new base C (0/40);  total (41/54)
Fisher's exact test for biased strand distribution p-value = 1.24e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786705‑787005                                                                                                                                                         |                                                                                                                                                        ccGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑a                                                                                                                                                         <  2:656411/151‑1 (MQ=255)   gTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑Acc                                                                                                                                                       >  2:1225016/1‑151 (MQ=255)   gTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑Acc                                                                                                                                                       >  2:689178/1‑151 (MQ=255)     cAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑Accc                                                                                                                                                      <  1:385772/150‑1 (MQ=255)      agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTg                                                                                                                                                    >  1:1100459/1‑151 (MQ=255)      agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGCCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTg                                                                                                                                                    >  2:313491/1‑151 (MQ=255)           ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGttt                                                                                                                                                <  2:80172/150‑1 (MQ=255)           ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTg                                                                                                                                               <  1:290670/151‑1 (MQ=255)             cGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑Ca                                                                                                                                            <  1:1151172/151‑1 (MQ=255)                aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAAc                                                                                                                                          <  1:1152255/150‑1 (MQ=255)                    gATGAATCACCACAGTCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAAccgcc                                                                                                                                      >  2:1324345/1‑150 (MQ=255)                     aTGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACcgccgc                                                                                                                                    >  2:1509601/1‑151 (MQ=255)                      tGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCgct                                                                                                                                   >  1:1406517/1‑151 (MQ=255)                       gAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCgctg                                                                                                                                  >  1:991039/1‑151 (MQ=255)                             ccacGGCCCTGATAAGTCGCGCAGCGTCGCATCAGGCATTGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCCCTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGCCTTTTACTTCGCTTT‑CCCCTGTTTTG‑CAACCGCCGCTGCtt                                                                                                                               <  2:205675/148‑1 (MQ=255)                                  ggCCTGATAGGTCGCGCGGCGTCGCATCAGGCAATGTCCTCCATTTTAAGCAACAAAAAAGCCGTCTCCCTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTCCTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGc                                                                                                                        <  1:270894/150‑1 (MQ=255)                                  ggCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACTTTTACTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGc                                                                                                                        <  1:469159/150‑1 (MQ=255)                                  ggCCTGATAAGGCGCGCGGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCg                                                                                                                       <  1:131510/151‑1 (MQ=255)                                  ggCCGGATAAGTCGCGCAGCGTCGCATCAGGCAAGGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCATTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGTTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGc                                                                                                                        <  2:705168/150‑1 (MQ=255)                                     cTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGAtc                                                                                                                    >  1:15091/1‑151 (MQ=255)                                       gATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                  >  1:536605/1‑151 (MQ=255)                                          aaGTCGCGCAGCGTCGCATCGGGAAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGTTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGAtct                                                                                                                   <  1:212160/147‑1 (MQ=255)                                          aaGTCGCGCAGCGCCGCATCGGGAAATGTGTTCCTTTGTTAGCAACAAAAAAGCCGACTCCCTTGCAGTCGTTTTCCCCATTTAAAA‑CGAATGACGTTTACTCCGCTTT‑CCCCTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                  <  2:740832/148‑1 (MQ=255)                                            gTCGCGCAGCGTCGCATCAGGAATTGTGCTCCTTTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTGGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑CCCCTGTTTTG‑AAACCGCCGCTGCTTTCGCTGCGATCTCGTCAg                                                                                                             <  2:162387/151‑1 (MQ=255)                                              cgcgCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                            >  2:280054/1‑150 (MQ=255)                                               gcgcGGTGTCGCATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGCCTTTTCCTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                            <  1:573290/149‑1 (MQ=255)                                               gcgaggCGTCGCATCAGGCATTGTGCTCCTTTTTTAGCAACAAAAAAGCCGACTCCCTGGCAGTCGTCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGTTTTG‑CACCCGCCGCTGCTTTCGCTGCGATCTCGTCa                                                                                                              <  2:140167/142‑1 (MQ=255)                                                  cAGCGTCGCATCGGGCATTGTGCTCCATTGTTAACAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTAAAA‑CGAATGACGTTTATTTCCCTTT‑CCCCTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTa                                                                                                        <  1:1477962/150‑1 (MQ=255)                                                    gCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGACTTTTACTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt                                                                                                          <  1:1554248/146‑1 (MQ=255)                                                           aTCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              >  1:1188232/1‑151 (MQ=255)                                                              ggccAATGTGCTCCTTTGTTGGAAAAAAAAAAGCCGCCGCTCTTCCGGTGGGCTTTCTCATTTAAAA‑CAAATGCCGTTTCCTCCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                           <  2:925247/148‑1 (MQ=255)                                                                gCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                         >  1:134071/1‑151 (MQ=255)                                                                acAATCTCCTCTATTGTTACCAAAAAAAAGGCCGCCTCCCTTTCGGTCGGCTTTCTCATTTTAAA‑CGAATGCCTTTTCCTGCCCTTT‑CCCCGGTTTTC‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                          <  1:30362/149‑1 (MQ=255)                                                                  atttttcTCCATTTTAAGCAACAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CAAATTACGTTTCCTTCGCTTT‑CCCGTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGttt                                                                                       <  1:59786/145‑1 (MQ=255)                                                                  aaTGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCACTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑CCCCTGTTTGG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGttt                                                                                       <  1:1216166/151‑1 (MQ=255)                                                                    tgtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                      >  1:1036124/1‑150 (MQ=255)                                                                     gtgCTCCATTGTTAGCAAC‑AAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACTTTTACTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCa                                                                                     <  1:1109486/149‑1 (MQ=255)                                                                        cTCCATTTTTACCAACAAAAAAGCCGCCCCCCTTGCAGCCGGCTTTC‑CATTTTAAACCGAATGCCTTTTCCTCCGCTTT‑CCCCTGTTTTG‑AAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc                                                                                   <  2:206944/149‑1 (MQ=255)                                                                         tctaTTGTTAACAAAAAAAAAGCCCACTCACTTGATGTCGGTTTTCTCATTTAAAA‑CGAATGACGTTTTCTTCGTTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc                                                                                   <  2:287777/145‑1 (MQ=255)                                                                          ccATTTTTAGCAAAAAAAAACCCGACTCCCTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGCCGTTTACTTCGTTTT‑CCCCGGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAg                                                                                    <  1:885926/146‑1 (MQ=255)                                                                          ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt                                                                               >  2:295764/1‑151 (MQ=255)                                                                           cATTGTTAGCCACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                              >  1:614618/1‑151 (MQ=255)                                                                            aTTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                             >  1:1230532/1‑151 (MQ=255)                                                                              tGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                             >  2:291292/1‑149 (MQ=255)                                                                                ttAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAt                                                                         >  2:954028/1‑151 (MQ=255)                                                                                ttAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAt                                                                         >  2:10780/1‑151 (MQ=255)                                                                                taaGCAAAAAAAAAGCCGACTCTCTTGCAGGCGGTTTTCTCATTTAAAA‑CGAATGACGTTTACTTCGCTTT‑CCCCTGTTTTG‑AACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                           <  2:859372/147‑1 (MQ=255)                                                                                  aGCAAAAAAAAAGCCGCCCCCCTTGAAGTCGTCTTTCTCATTTAAAA‑CGATTGCCTTTTCCTTCTTTTT‑ACCTTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAtt                                                                        <  2:448904/150‑1 (MQ=255)                                                                                     aacaaAAAGGCCGACTCACTTGCAGCGGACTGCCCCATTTCAAA‑CGAATGACTTTAACTGCGCTTT‑GCCCTGTTTGC‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaaa                                                                          <  1:106950/145‑1 (MQ=255)                                                                                      aaaaaaaaCCCGACTCACTTGCATTCGGTTTTCTCATTTTAAA‑CGAATGACGTTTCCTTCGCTTT‑CCCTTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                      <  2:923227/146‑1 (MQ=255)                                                                                        aaaaaaGCCGACTCCCTTGCGGTCGTTTTTCTCATTTTAAA‑CGAATGACGTTAACTCCGTTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                                 <  1:1131398/151‑1 (MQ=255)                                                                                        aaaaaaGCCGACTCCCTTGCAGGCGGCTTTCTCATTTAAAA‑CGAATGCCTTTTCCTTCGCTTT‑CCCCTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                                 <  2:769999/151‑1 (MQ=255)                                                                                        aaaaaaCCCCACTCTCTTCCACGTGGCTTTCTAATTTTAAA‑AGAAAGCCGTTT‑CATTGCTTTAACCCTGTTTCC‑GAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                      <  2:1203597/146‑1 (MQ=255)                                                                                         aaaaaGCCGACTCCCTTCCAGTCGGTTTTCTCATTTTAAA‑CGATTGACTTTTACTTCGCTTT‑CCCCTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGa                                                                <  1:1574/151‑1 (MQ=255)                                                                                         aaaaaGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGCCTTTTACTTCGTTTT‑CCCCTGGTTTG‑AAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGa                                                                <  2:693757/151‑1 (MQ=255)                                                                                           aaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                              >  1:90455/1‑151 (MQ=255)                                                                                           aaaGCCGACTCACTGACGGTCGTTTTTCTCATTTAAAA‑CAAATGACTTTTACTTCGTTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                              <  2:1546741/151‑1 (MQ=255)                                                                                            aaGCCGACTCATTTGCAGTCGGCTTTCTCTTTTTAAA‑CGAATGACGTTTCCTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                             <  1:183446/151‑1 (MQ=255)                                                                                              gCCGACTCACATGCAGTGGGCGTTCTCATTTTAAA‑CGAATGCCGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCACCATTACCCAGATAATAGCGTTTCAGCGGTTTGACATTCTCGTCGAACTc                                                            >  2:393229/1‑150 (MQ=255)                                                                                                  aCTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATac                                                        >  2:1443800/1‑150 (MQ=255)                                                                                                   cTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                      >  1:1416987/1‑151 (MQ=255)                                                                                                     cACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     >  1:1455655/1‑150 (MQ=255)                                                                                                     cACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     >  2:673684/1‑150 (MQ=255)                                                                                                       cTTGCAGTCGGTTTTCTCATTTTAAA‑CGAATGACTTTTCCTTCGTTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATaca                                                       <  1:604876/146‑1 (MQ=255)                                                                                                        ttGCAGTCGGTTTTCTCATTTTAAA‑CGAATGACTTTTCCTTCGCTTT‑CCCCTGTTTTG‑AAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATaca                                                       <  1:1462744/145‑1 (MQ=255)                                                                                                           cAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGc                                               >  1:11983/1‑150 (MQ=255)                                                                                                            aGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                             >  2:333412/1‑151 (MQ=255)                                                                                                              tCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGc                                           >  2:1064167/1‑151 (MQ=255)                                                                                                                ggttttCTCATTTTAAC‑CGACTGACTTTTACTTCGTTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                          <  2:222235/147‑1 (MQ=255)                                                                                                                ggttttCTCATTTTAAA‑CGAATGACGTTTCCTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCa                                         <  1:883758/148‑1 (MQ=255)                                                                                                                ggCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                          >  1:742432/1‑150 (MQ=255)                                                                                                                   tttCTCATTTAAAA‑CGAATGACTTTTCCTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                       <  1:666666/150‑1 (MQ=255)                                                                                                                    ttCTCATTTTAAA‑CGAATGACTTTTCCTTCGTTTT‑CCCCTGGTTTG‑AAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                      <  1:1351369/150‑1 (MQ=255)                                                                                                                    ttCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                     >  2:59829/1‑151 (MQ=255)                                                                                                                      ctcATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                                   >  2:1395233/1‑151 (MQ=255)                                                                                                                       tcATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                                   >  1:351883/1‑150 (MQ=255)                                                                                                                          ttttAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGatat                               >  2:1285419/1‑151 (MQ=255)                                                                                                                          tttaAAA‑CGAATGCCGTTTACTTCGCTTT‑CCCTTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                      <  2:129834/144‑1 (MQ=255)                                                                                                                           tttAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                              >  2:630298/1‑151 (MQ=255)                                                                                                                            ttAAA‑CGAATGACTTTTACTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                              <  1:910927/150‑1 (MQ=255)                                                                                                                             taaaccaaaTGCCTTTTACTTCTCTTT‑ACCCTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                              <  2:1555491/143‑1 (MQ=255)                                                                                                                               aaccaaaTGACGTTTACTTCGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                           <  2:597284/146‑1 (MQ=255)                                                                                                                                   gAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                        >  1:834240/1‑150 (MQ=255)                                                                                                                                    aaTGACGTTTACTTCGCTTT‑ACTCTGGCTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACAAGCGGCACGCCAGTCGGGATATTAAGCCCa                       >  2:124541/1‑150 (MQ=255)                                                                                                                                    aaTGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAGGCTCa                       >  1:751599/1‑150 (MQ=255)                                                                                                                                       gACGTTTACTTCGTTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                    <  2:239151/150‑1 (MQ=255)                                                                                                                                        ccttttACTCCGCTTT‑‑CCCCGGTTTGCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAg                     <  2:1164705/145‑1 (MQ=255)                                                                                                                                           tttACTTCGTTTT‑CCCCTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc                 <  1:114603/149‑1 (MQ=255)                                                                                                                                              aCTTCTCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt             <  2:704389/150‑1 (MQ=255)                                                                                                                                               ctccGCTTT‑CCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATctct              <  1:185327/145‑1 (MQ=255)                                                                                                                                               cTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct           >  1:470093/1‑151 (MQ=255)                                                                                                                                                 tCGCTTT‑CCCCTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt          <  2:937360/150‑1 (MQ=255)                                                                                                                                                    cTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc       >  1:1289214/1‑150 (MQ=255)                                                                                                                                                     ttt‑CCCCTGTTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt      <  1:1530929/150‑1 (MQ=255)                                                                                                                                                     ttt‑ACCCTGGTTGG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc         <  1:1376645/147‑1 (MQ=255)                                                                                                                                                       t‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCAt   >  2:780557/1‑151 (MQ=255)                                                                                                                                                         aCCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  1:153146/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                        CCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAA‑CGAATGACGTTTACTTCGCTTT‑ACCCTGGTTTG‑CAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786705‑787005

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