Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	50.0%	103.6 / 30.2	82	G566G (GGT→GGG)	gyrA	DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (38/3);  new base C (0/41);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 6.24e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTAC  >  NC_000913/2337574‑2337868                                                                                                                                                      |                                                                                                                                                  ccGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTa                                                                                                                                                   >  2:1151265/1‑150 (MQ=255)  cgcgcgcgCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAcc                                                                                                                                                 <  1:1099975/151‑1 (MQ=255)      cgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTc                                                                                                                                             <  2:1056180/151‑1 (MQ=255)        cgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTcc                                                                                                                                            >  1:1457710/1‑150 (MQ=255)            aCGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCa                                                                                                                                       >  2:1280195/1‑151 (MQ=255)                       cggccAATTGAAAAACTTTCACAGAATGGAGGCCCCCAGGGCGGGAGAAGCACAAAATAGGGTGGGGAGTGTCCGCCACCAGCAGTCGGTCAAAAAAGTCTTCTTCTTTAAGGCGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTgcgc                                                                                                                            <  1:1406200/147‑1 (MQ=255)                          gCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCtt                                                                                                                          >  2:1427179/1‑150 (MQ=255)                            aattaaaaaaCTTTCATCGAAAAGACGCGCCCCCGGCTGGAAAAGCAAAAAATATGGTCGGGAGTTTTCCCCCCCAGCAGTCGGTCGAAAAAGTCTTCTTTTTTAATACGTGCGGCAGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGt                                                                                                                       <  2:1247721/144‑1 (MQ=255)                             aCTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTa                                                                                                                      >  1:203955/1‑151 (MQ=255)                              cTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTAt                                                                                                                     >  2:1336993/1‑151 (MQ=255)                              cTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTAt                                                                                                                     >  2:1307359/1‑151 (MQ=255)                              cTGAAAAACTTTCATCGAATAGACGCGCCCACGGCGGGAAAAGCACAGAATATGGTGGTGAGTTTTCGCCACCAGAAGTGGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTAt                                                                                                                     <  2:1109847/151‑1 (MQ=255)                                gATAACCTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGACTTACCTTTCCCGCCACGACGCTGCGCTTCGTAtt                                                                                                                    >  1:1221907/1‑150 (MQ=255)                                aaTAACCTTCCACCAAATGGACGCCCCCACGGCTGGAGAAGCAAAGAATATGGTCGTGGGTTTCCGCCACCAGCAGTCGGTCAATAAGGTCTTCTTCTTTAATACGGGCGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                   <  2:1087970/150‑1 (MQ=255)                                   aaccTTTCATGGAATAGCCGGACCCCCGGCTGGGGANNNNANNAATGTGGGCGTGAGTGTCCCCCCCCACCAGTCGTTCGATAAATTCTTTTTCTTTATTACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGa                                                                                                                <  1:1457116/148‑1 (MQ=255)                                              gAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGAt                                                                                                     >  2:535551/1‑151 (MQ=255)                                              gAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGAt                                                                                                     >  1:1510379/1‑151 (MQ=255)                                                aTAGACGCCACCGCGGGGGGAAAAGAACAGAATATGGTGGGGAGTGTCCGCCACCAGAGGTCGGTCAAAAAATTCTTCTTTTTAAATCCGGGCGGAAAATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATAc                                                                                                   <  1:229649/151‑1 (MQ=255)                                                   gACGCGACCACGGCTGGCGAAGCACAGAATATGGTCGTGAGTGTTCCCCACCAGCAGTCGGTCGCTAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTa                                                                                                >  1:863195/1‑151 (MQ=255)                                                   gACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTa                                                                                                >  2:485114/1‑151 (MQ=255)                                                   gACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGGCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTa                                                                                                >  1:1395449/1‑151 (MQ=255)                                                     cgcgCCCACGGCTGGAGAAGAAAAAAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGAAAAAGTCTTCTTTTTAAAGACGTGCGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                                                                              <  1:88191/151‑1 (MQ=255)                                                       cgcccccggcggGAGAACCACAAAATAGGGTCGGGGGGGTTCGCCCCCAGCAGCCGGTCGAAAAAGTTTTTTTTTTTAATACGGGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                            <  2:611541/145‑1 (MQ=255)                                                       cgcccACGGCGGGAGAAACACAGAATTTGGTCGTGAGTTTTCGCCCCCAGCAGTCGGTCGATAAAGTTTTTTTTTTAAATCCGGGGGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                            <  2:369616/148‑1 (MQ=255)                                                       cgACCCCGGCGGGAAAAGCACAGAATATGGTGGGGAGTGTTCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGGGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                            <  2:1025135/151‑1 (MQ=255)                                                       cacccAGGGGTGGAGAAGAACAAATTATGTTCGTGAGTGTGCGCCACCAGCAGTCGGTAGATAAGGTCTTCTTTTTTAATACGGGCGGCAGTTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                            <  1:486950/148‑1 (MQ=255)                                                              ggCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGcc                                                                                        >  1:1031962/1‑148 (MQ=255)                                                                  ggAGAAGCAGAGAATTTGTGGGGGGGGGTCCCCCCCCAGCTGTCGGTCAATAAGGTTTTTTTTTTTAAGCCGGGGGGAGTTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                  <  1:186676/150‑1 (MQ=255)                                                                  ggAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTga                                                                                 >  1:160867/1‑151 (MQ=255)                                                                  ggAGAAGCAAAAAATTTGGTCGTGAGTGTTCCCCCCCAGCAGTCGGTCGATAAAGTTTTCTTCTTTAATGCGGGCGGCAAATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                  <  2:344874/150‑1 (MQ=255)                                                                  ggAAAAGCAAAAAATAGGGTGGGGAGTTTTCCCCCCCAGCAGTCGGTCAATAAAGTCTTTTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                  <  2:989168/150‑1 (MQ=255)                                                                    agaAGCAAAGAATAGGGTGGTGATTGTCCCCCCCCAGCGTTCGGTCGATAAAGTTTTCTTTTTTAATACGGGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgag                                                                                <  2:1386982/150‑1 (MQ=255)                                                                    aaaaGCACAAAATAGGGCCGTGAGTGTCCCCCCCCAGCAGTCGTTCGAAAAAGTTTTTTTTTTTAATCCGGGGGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgag                                                                                <  2:926012/148‑1 (MQ=255)                                                                     aaaGCACAGAATATGGTCGGGAGTGTTCCCCCCCAGCAGTCGGTCGAAAAAGTCTTCTTTTTTAATACGTGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagag                                                                              <  2:452097/150‑1 (MQ=255)                                                                      aaGCGCAGAATATGGTCGGGTGTGTTCGCCCCAAGCGTCCGGTTATAAAAGTCTTTTTTTTTATTTCGGGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagag                                                                              <  1:391051/150‑1 (MQ=255)                                                                        gCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag                                                                            >  1:917011/1‑150 (MQ=255)                                                                        gCACAGAATATGGTCGGGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGCCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag                                                                            >  2:610727/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAGGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                       >  2:365764/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                       >  2:1342371/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                       >  1:492082/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                       >  2:1025958/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                       >  2:277689/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                       >  1:41110/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                       >  1:1423689/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                       >  1:360695/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGGGAGAGAGCGTCa                                                                       >  1:1083186/1‑150 (MQ=255)                                                                             gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGCGAGAGAGCGTCa                                                                       >  1:294365/1‑150 (MQ=255)                                                                               atatGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGa                                                                    >  2:530743/1‑151 (MQ=255)                                                                               atatGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGa                                                                    >  1:1053213/1‑151 (MQ=255)                                                                               atatGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGa                                                                    >  1:866837/1‑151 (MQ=255)                                                                                   ggTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACca                                                                 >  1:1450692/1‑150 (MQ=255)                                                                                    ttcgggggttttCCCCCCCGGCAGTCGGTCAAAAAAGTCTTTTTTTTAAATCCGGGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcac                                                                <  2:1400756/140‑1 (MQ=255)                                                                                     tCGTGAGTTTCCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATACGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACAt                                                              <  2:41136/151‑1 (MQ=255)                                                                                       gTGAGTGTCCCCCCCCAGCAGTGGGTCAAAAAAGTCTTTTTCTTTAATACGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATc                                                             <  2:371255/150‑1 (MQ=255)                                                                                          agtttccccccccaGCGGTCGGTCGATAAATTCTTTTTTTTTAATCCGTGGGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTc                                                          <  2:142096/139‑1 (MQ=255)                                                                                           gtgtTCCCCCCCAGCAGTGGGTGAAAAAAGTCTTCTTCTTAAATACGTGCGGCAGTTTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTcc                                                         <  2:1386408/150‑1 (MQ=255)                                                                                            tgtTCGCCCCCAGCAGTCGGTCAAAAAAGTCTTCTTTTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt                                                        <  2:360802/150‑1 (MQ=255)                                                                                              ttCGCCACCCGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCATggg                                                     >  1:275273/1‑151 (MQ=255)                                                                                              tcccccccgagcagTCGGTCGAAAAATCTTTCTTCTTAAATGCGTGCGAAATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                     <  1:71240/142‑1 (MQ=255)                                                                                              tccccccccaGCGGTGGGTCAATAAAGTCTTTTTCTTTAATCCGGGCGGCAATTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                     <  2:1032470/144‑1 (MQ=255)                                                                                              tccGCCCCCGGCAGCCGGTCAAAAAAGTTTTCTTCTTAAATACGGGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                     <  1:1279583/149‑1 (MQ=255)                                                                                              tccGCCACCAGCGGTGGGTCGAAAAATTTTTTTTCTTTAAGACGTGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                     <  1:1414364/149‑1 (MQ=255)                                                                                                   ccccggcagTCGTTGAATAAAGTCTTCTTTTTAAATCCGGGGGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                <  1:594745/146‑1 (MQ=255)                                                                                                    cccaGCTGTCGGTCGATAAATTCTTTTTCTTTAATCCGTGCGGCAATTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCa                                               <  1:679251/150‑1 (MQ=255)                                                                                                         cggtcggtcaATAAATTTTTCTTTTTTAATACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATc                                           <  2:934995/148‑1 (MQ=255)                                                                                                            tcgttcgagaagctattctttttTATTGCGTCCGGCATATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTc                                        <  1:724063/135‑1 (MQ=255)                                                                                                            tcggtcgAAAAAGTTTTTTTTTTTAATACGGGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                       <  2:450099/151‑1 (MQ=255)                                                                                                            tcggtcgAAAAAGTCTTCTTCTTTAATGCGTGGGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                       <  2:1086676/151‑1 (MQ=255)                                                                                                             cggtcgATAAAGTCTTTTTTTTTAATCCGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCa                                      <  2:1262746/151‑1 (MQ=255)                                                                                                                tcgATAAAGTCTTTTTTTTTAATACGTGCGGAAGTTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                    <  1:1336335/150‑1 (MQ=255)                                                                                                                tcaaTAAAGTCTTTTTTTTTAATACGGGCGGAAAATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGt                                   <  2:239398/148‑1 (MQ=255)                                                                                                                     aaaGTCTTTTTTTTTAATCCGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg                              <  2:720780/151‑1 (MQ=255)                                                                                                                        gTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGCCGCTGCGCTTCGTATTCAGAAAGCGGCTGCTCCCTACCGTAGCCCTGGTGCGAGAGCGTCACGACCACATCTTCCTGGCCGATCAGATCTTCCAGGTTGATGt                             >  1:301468/1‑149 (MQ=255)                                                                                                                              ttttttAATACGTGCGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTg                     <  2:268766/148‑1 (MQ=255)                                                                                                                              ttttttAATACGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTg                     <  2:340605/148‑1 (MQ=255)                                                                                                                               tcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg                        >  2:1212002/1‑147 (MQ=255)                                                                                                                               tcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg                        >  1:529973/1‑147 (MQ=255)                                                                                                                                ctctAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGtt                   >  2:1200035/4‑151 (MQ=255)                                                                                                                                                 taTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTac  <  1:1313994/150‑1 (MQ=255)                                                                                                                                                 gaattACCTTTCCCGCGATGACGAAGCGCTTCGTATTGAGAAAGCGGCTGATACTTAACGTAGCACTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGCGAGTTCAGTa   >  1:460341/4‑150 (MQ=255)                                                                                                                                                 gATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTa   >  2:822785/1‑150 (MQ=255)                                                                                                                                                 gATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTa   >  1:167669/1‑150 (MQ=255)                                                                                                                                                 aaTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTac  <  2:1083710/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTAC  >  NC_000913/2337574‑2337868

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