Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,147,733	0	T	G	42.5%	142.9 / 10.6	80	G20G (GGT→GGG)	yegD	Hsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base T (10/36);  new base G (34/0);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 5.76e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCG  >  NC_000913/2147583‑2147879                                                                                                                                                       |                                                                                                                                                   gACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGt                                                                                                                                                    <  2:1272771/151‑1 (MQ=255)  aCGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTa                                                                                                                                                   <  2:1150812/151‑1 (MQ=255)    gCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGaa                                                                                                                                                  >  1:86470/1‑150 (MQ=255)      tGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGaaa                                                                                                                                                 >  2:99424/1‑149 (MQ=255)      tGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAAcc                                                                                                                                               >  1:426105/1‑151 (MQ=255)       gcgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAAc                                                                                                                                                <  2:487851/149‑1 (MQ=255)       gcgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCg                                                                                                                                              <  1:1329043/151‑1 (MQ=255)        cgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGGGCTACGGGTAACCGc                                                                                                                                             >  1:397736/1‑151 (MQ=255)         gctttgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGGAAAACGc                                                                                                                                             >  2:755065/1‑150 (MQ=255)          ctttgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTAAGTGTGACGGGAAAACGCa                                                                                                                                            >  1:730821/1‑150 (MQ=255)             tgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCAttt                                                                                                                                         <  1:142099/150‑1 (MQ=255)                tttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCATCGGGACGGGaaaaagccttttt                                                                                                                                       >  2:551464/1‑139 (MQ=255)                  tCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGGCGGTAAAACGCATTTTCtta                                                                                                                                    >  1:833868/1‑148 (MQ=255)                    atatTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTaaaaa                                                                                                                                 <  1:451929/151‑1 (MQ=255)                       ttCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGGAAAACCCATTTTTTaaaaa                                                                                                                                 >  2:568108/1‑148 (MQ=255)                        tCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGGCGGGAAAACCCATTTTCTAAAAATgg                                                                                                                              >  2:582278/1‑150 (MQ=255)                         ccGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTATGCGGGAAGGGAAAACGCCTTTGCTtaaaa                                                                                                                                 >  2:86215/1‑146 (MQ=255)                          cGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGCCATTAGGGGCGGGAAACCCCAGTTGCTtaaaaag                                                                                                                               >  1:1086860/1‑145 (MQ=255)                              tGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACg                                                                                                                        <  2:433304/150‑1 (MQ=255)                              tGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACg                                                                                                                        <  2:457063/150‑1 (MQ=255)                                  gcgACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCCATTTCTAAAAATAGaaaa                                                                                                                          >  2:851711/1‑144 (MQ=255)                                    gACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCAAGGGGGGCGGGAAAACGCATTTTCTTAAAAAGGaaaaaggaa                                                                                                                     >  1:195022/1‑142 (MQ=255)                                      cGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGTGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAGCACgc                                                                                                               >  1:364840/1‑151 (MQ=255)                                      cGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACgc                                                                                                               <  2:1057351/151‑1 (MQ=255)                                       ggCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACgct                                                                                                              <  2:1132233/151‑1 (MQ=255)                                         cAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGTCTACGTGTTGGTGAACCGCCAGTTCTAAAAAAGGTAAAcaaaaagagcgccc                                                                                                            >  2:455205/1‑138 (MQ=255)                                          aaCATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTATTCGTGGAGGGAAAACGCATTTTCTTAAAAAGGAAAACGACAGCACgctgct                                                                                                           >  1:1193172/1‑151 (MQ=255)                                          aaCATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACgctgct                                                                                                           <  2:235466/151‑1 (MQ=255)                                            cATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAAACGCATTTTCTAAAAAAGGAAAAAGACAGCACgc                                                                                                               >  1:420226/1‑145 (MQ=255)                                              tttCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCCTTTTCTAAAAATGGAAAACGACAACACCCTGCTGCCt                                                                                                       >  1:232897/1‑151 (MQ=255)                                                 cGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGTAAAACGCAATTGCTTAAAAAGGAAAAAGAAAGCACGCTGCTGCCTTc                                                                                                     >  1:1076291/1‑150 (MQ=255)                                                  gCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGGAGGGAAAACCCATTTTTTCAAAAAAGAAAACGACAGCAAGccgctg                                                                                                          >  2:873816/1‑144 (MQ=255)                                                     aaaGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAAt                                                                                                  <  1:887668/149‑1 (MQ=255)                                                       aGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCtt                                                                                              <  2:1332155/151‑1 (MQ=255)                                                             gCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGGGGCGGGGAAACCCCTTTTCTAAAAATGGAAAAAGACAACACGCTGCTGCCTTCCATGCttt                                                                                             >  2:1405743/1‑146 (MQ=255)                                                              ccccTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGGGGGGCGGGAAAACGCCTTTTCTAAAAAATGAAAAAGACAGCACGCCGCTTCCCTCAATGCTTTgcgcgc                                                                                       >  1:1214736/1‑151 (MQ=255)                                                                   tcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAAc                                                                                   <  2:1176505/150‑1 (MQ=255)                                                                     ttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGGGACGGGAAAACGCCTTTTTTAAAAAAAGAAAACGACAGCACGCCGCTGCCTTTCATGCTTTTCGCGccca                                                                                    >  2:958516/1‑145 (MQ=255)                                                                      tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACCCCTTTTCTAAAAATGGAAAACGACAGCACGCCGCCGCCTTCAATGCTTTTCGCGCCCAcccgg                                                                               >  2:697169/1‑150 (MQ=255)                                                                      tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAAcgcg                                                                                <  2:1252265/150‑1 (MQ=255)                                                                         ccGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGaa                                                                            <  2:1113198/151‑1 (MQ=255)                                                                          cGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGTGACGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAACACGCTGCCGCCTTCAATGCTTTGCGCGCCAACGCGTg                                                                              >  2:1248901/1‑148 (MQ=255)                                                                          cGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAg                                                                           <  2:1321928/151‑1 (MQ=255)                                                                            gCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGTGACGGGAAAACGCATTTTCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGaa                                                                            >  1:865884/1‑148 (MQ=255)                                                                             cATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGGAAAACGCCTTTTCTAAAAATGGGAAAAGACAGCCCCCTGCTGCCCTCAATGCTTTTCGCGCCCACGCGGGAAACgg                                                                        >  1:1076124/1‑151 (MQ=255)                                                                                     ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAACCGCATTTGCTAAAAAAGGAAAACAACAGCACGCTGCTGCCCTCAATGCTTTTCGCGCCCACgcgggaagcgggaa                                                                     >  2:600388/1‑146 (MQ=255)                                                                                      tattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGGGACGGGGAACCGCATTTTCTAAAAATGGAAAAAGACAGCACGCTGCTGGCTTCAATGCTTTTCGGGCCCAAGCGGGAAGCgggta                                                                     >  2:945513/1‑142 (MQ=255)                                                                                      tattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGaa                                                                <  2:651992/150‑1 (MQ=255)                                                                                           cGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCt                                                           <  2:599111/150‑1 (MQ=255)                                                                                            gtgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGTAAACCGCAATTGCTAAAAATGGAAAACGACAGCACGCCGCTGCCTTCCATGCTTTGCGCGCCAACGCGTGAAGAGGGAAGAGAATGGCTGt                                                         >  1:963229/1‑151 (MQ=255)                                                                                             tgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTCCATGCTTTTCGCGCCCACGCGTGAAGCGGGAAGCGAAAGGCTg                                                          >  1:229153/1‑149 (MQ=255)                                                                                              gtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTa                                                        <  1:1232009/150‑1 (MQ=255)                                                                                              gtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAACGCGTGAAGCGGGAAGCGAATGGCTGTAc                                                       >  1:368028/1‑151 (MQ=255)                                                                                                  aTTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACCCATTTGCTAAAAAAGGAAAACCACAGCCCCCAGCTGTCCTCTATGTTTTGCGCGCCCACGCGCGGAGGGGTAAGGGAATGGCTGTgccgcc                                                   >  1:186122/1‑151 (MQ=255)                                                                                                  aTTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGc                                                    <  2:1133481/150‑1 (MQ=255)                                                                                                     ggTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATAAGGGGCGGGAAACCGCGATTTCTAAAAAAGGGAAAAGGAAGACCCCGGCTGCCTTCAATTCTCTGCGCGCCACCGCGTGGATCGGGTAGCGAAAGGGCGTACCGc                                                    >  1:1181494/1‑147 (MQ=255)                                                                                                        tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCcatcat                                              <  2:442400/150‑1 (MQ=255)                                                                                                        tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCcatcat                                              <  1:50540/150‑1 (MQ=255)                                                                                                           gATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgatg                                          <  2:1233971/151‑1 (MQ=255)                                                                                                            aTTACGGTCCAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgat                                           <  1:1416452/149‑1 (MQ=255)                                                                                                            aTTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGGCGGGAAAACGCATTTTCTAAAAATGGAAAAAGGCAGCACGCTGCTGCCTTTAAAGCTTTGTGCGGCAAAGCGTGAAGCGGGAAGAGGATGGGTGTACCGCCATCatgaag                                          >  2:615966/1‑148 (MQ=255)                                                                                                                cGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCCATGCTTTTCGCGCCCACGCGTGAAGCGGGAAGCGGAAGGCTGTACCGCCCTCATGATGTTcc                                      >  2:937005/1‑150 (MQ=255)                                                                                                                  gTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGc                                   <  2:1327380/151‑1 (MQ=255)                                                                                                                  gTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGc                                   <  1:755084/151‑1 (MQ=255)                                                                                                                    aCAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCTTTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCa                                  <  2:839326/150‑1 (MQ=255)                                                                                                                    aCAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCa                                  <  2:1053883/150‑1 (MQ=255)                                                                                                                        cAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTTTAAAAAAGGAAAACGACAACACCCCGCTGCCTTTAATGCTTTTCGCGCCCACGCGGGAAGCGGGAAGCGAATGGCTGTACCCCCCTCATGATGtt                                        >  2:215349/1‑140 (MQ=255)                                                                                                                          aaCTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGacgacg                           <  1:877614/151‑1 (MQ=255)                                                                                                                             tGTTCAGTGGCGGTCATGCGTGGCGGTAAAACGCATTTGCTAAAAATGGAAAAAGACAGCACGCTGCCTCCTTCAATGCTTTTTGCGCCAACGCGGGAAGCGGTAAGCGAATGGCTGTACCGCCATCCTGATGTTTCGGCAgacgac                            >  2:993736/1‑147 (MQ=255)                                                                                                                             tGTTCAGTGGCGGTCATGCGTGACGGGAAAACCCATTTTCTAAAAAAGGAAAAAGACAGCACGCCGCTGCCTTCAAAGCTTTGCGCGCCCACGCGGGAAGCGGGAAGCGAAAGGCTGTACCG‑CATCATGATGtt                                        >  2:523860/1‑134 (MQ=255)                                                                                                                              gTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCTGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGAt                         <  1:711541/149‑1 (MQ=255)                                                                                                                              gTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATg                        <  1:208359/150‑1 (MQ=255)                                                                                                                                 cAGTGGCGGTCATGCGTGGCGGGAAACCGCATTTGCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCCACGCGGGAAGCGGGAAGCGAAAGGCTGTACCGCCCTCATGAAGTTTCCGCAGACGACGATGaaa                     >  2:998950/1‑150 (MQ=255)                                                                                                                                     ggCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAg                <  2:1004147/151‑1 (MQ=255)                                                                                                                                      gCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAgg               <  1:568121/151‑1 (MQ=255)                                                                                                                                       cGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGgc              >  2:411968/1‑151 (MQ=255)                                                                                                                                            aTGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGt         <  2:86469/151‑1 (MQ=255)                                                                                                                                               cGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTAc      <  2:420215/151‑1 (MQ=255)                                                                                                                                                gTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACg     <  1:99426/151‑1 (MQ=255)                                                                                                                                                  gACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTc   >  1:796131/1‑151 (MQ=255)                                                                                                                                                   aCGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCg  >  2:300779/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                   GACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCG  >  NC_000913/2147583‑2147879

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: ‑