Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	115,797	0	C	A	43.3%	93.6 / 38.8	67	E435* (GAA→TAA)	hofB	T2SE secretion family protein; P‑loop ATPase superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (28/10);  new base A (0/29);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.15e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTA  >  NC_000913/115652‑115943                                                                                                                                                   |                                                                                                                                                   cccATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCaaaa                                                                                                                                                >  2:424824/1‑150 (MQ=255)  ccATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAg                                                                                                                                              <  2:740452/151‑1 (MQ=255)   cATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGc                                                                                                                                             <  2:522226/151‑1 (MQ=255)       gCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTAc                                                                                                                                         <  1:673396/151‑1 (MQ=255)           cccATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCAt                                                                                                                                     >  1:300676/1‑151 (MQ=255)             cATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATAc                                                                                                                                   >  1:532631/1‑151 (MQ=255)               tCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATAccc                                                                                                                                 >  2:447093/1‑151 (MQ=255)                cTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCg                                                                                                                                >  2:758051/1‑151 (MQ=255)                  tGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATAccc                                                                                                                                 <  2:221614/148‑1 (MQ=255)                   gcgcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATAccc                                                                                                                                 <  1:499857/147‑1 (MQ=255)                      cATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                                                                                          >  1:670150/1‑151 (MQ=255)                                     ttATGCCATCCCAGCGCAAGAGTTGCTAACCCGCCATGCGCCATCCCCAGTCCGCGGATTACCTCTTAAAA‑GGTGGTTAAGCCTTGCTCTCCGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                                                                                              <  1:849953/148‑1 (MQ=255)                                       atgccatgccAGCGCCAGAGTTGCTTACACGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGc                                                                                                         >  2:714214/1‑151 (MQ=255)                                        tgcctggccaGCGCCAGAGTGGCTTACCCCCCATGCGGCCCCCCCAGGCCGCGGATAACCTCTAAAAA‑GGTGGTAAAGCCTTGCTCTCCGGCCGGCCGGCCTTTTTAAAAAAGCGTACGAACACCCGCCTGTCGTGCGTGCGTTTCCAg                                                                                                          <  2:331747/143‑1 (MQ=255)                                          ccatgccgGCGCAGGGTTTGCTTCCTCGCCAGGCGGCACCCCCAGTCCGCGGATTACCTCTTAAAG‑GGTGTTAAAGCCTTGTTCTACGGCCGGGCGGCCGTTTTAAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCg                                                                                                        <  2:922992/149‑1 (MQ=255)                                            atgccaGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTc                                                                                                    >  2:866655/1‑151 (MQ=255)                                                caGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACg                                                                                                >  2:127433/1‑151 (MQ=255)                                                   cgcCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCg                                                                                             >  2:416820/1‑151 (MQ=255)                                                     ccAGATTTGCTTACTCGCCATGCGGCACCCCCAGTCCGCGGATTACCTCTTAAAA‑GGTGGTTAAGCCTTCCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCgg                                                                                            <  1:654510/150‑1 (MQ=255)                                                     ccAGAGTTGCTTACTCCCCTTGCGGCCCCCCCAGTCCGCGGATTATCT‑TTCAAAGGGTGGTTAAGCTTTGCTCTCCGGCCAGGCACCGTTTTTAAAAAAGCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCg                                                                                             <  1:757604/149‑1 (MQ=255)                                                       agatTTGCTTACTCGCCATGCGGCCCCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                                                                                          <  1:22168/150‑1 (MQ=255)                                                       agagTGCCTTCCCCGCCAGGGGGCACCCCAATTCCGCGAATTAACTCTTAAAA‑GGGGTTAAAGCTTTGCCCGCGGGCCGGGCGGCCGTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGt                                                                                           <  1:253355/149‑1 (MQ=255)                                                        gatttGCTTACTCGCCATGCGCCACCCCCAGTACGGGATTTACCTCTTAAAA‑GGTGTTAAACCCTTCCTCTACGGCCGGGCGGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                                                                                          <  2:45741/146‑1 (MQ=255)                                                        gagTTGCTTACTCGCCATGCGGCACCCCCAGTCCGGGAATTAACTCTTAAAA‑GGGGTTTAAGCTTTGCTCTACGCCCAGGCGGCCTTTTTAAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCgg                                                                                            <  1:461781/147‑1 (MQ=255)                                                          gttgttTACTCGCCATGCGGCCCCCCCAGTACGCGGATTACCTCTTCAAA‑GGTGGTTAGGCTTTGCTCTACGGCCAGGCACCCTTTTTAAAAAAGCGTACACATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTa                                                                                       <  1:2722/150‑1 (MQ=255)                                                           ttgcttcccccacATGCGGCACCCCCAGTACGCGGTTAAACTCTACAAG‑GGTGGTAACGCCTTGCTCTCCGCCCAGCCAGCCGTTTTAAAAAAGCGTACCAATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAtt                                                                                        <  2:652466/136‑1 (MQ=255)                                                           ttGCTTACTCGCCAGGCGGCACCCCCAGTCCGCGGATTAACTCTTAAAA‑GGTGTTAAAGCCTTGCTCTCCGGCCGGGCAGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAg                                                                                      <  1:276421/150‑1 (MQ=255)                                                                tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaa                                                                                 >  2:904743/1‑150 (MQ=255)                                                                    cGCCATGCGGCACCCCCAGTACGCGGATTACCTCTTAAAA‑GGTGTTTAAGCCTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAg                                                                            <  1:677462/151‑1 (MQ=255)                                                                     gcaaTGCGGCACCCCCAGTACGCGGATAAACTCTAAAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAg                                                                            <  2:747470/147‑1 (MQ=255)                                                                      ccATGCGGCATCCCCAGTACGCGATTTATCTCTTAAAG‑GGTGGTTAAGCTTTGCTCTGCGGCGAGGCACCGGTTTTCAAAAGGCGTCCGAATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt                                                                          <  1:134736/151‑1 (MQ=255)                                                                      ccATGCGGCACCCCCAGTCCGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGc                                                                           <  1:906866/150‑1 (MQ=255)                                                                         tGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGAc                                                                       >  1:123926/1‑151 (MQ=255)                                                                            ggcccccccAGTCCGCGAATTACCTCTTCAAG‑GGTGTTTAAGCCTTGCTCTCCGGCCGGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGAc                                                                       <  1:304198/143‑1 (MQ=255)                                                                            ggCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGa                                                                     >  1:485769/1‑150 (MQ=255)                                                                            ggCACCCCCGGTGCGCGGTTTATTTCTCCAAG‑GGTGGTTAACCTTTGCTCTCCGGCCGGCCAGCTGTTTTCAAAAGGCGTCCGCATAACCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACg                                                                      <  2:553790/149‑1 (MQ=255)                                                                            ggCACCCCAAGTCCGGGTATTACCCCTTCAAA‑GGTGTTTAGGCTTTGCTCTGCGGCAAGGCGGCCTTTTTCAAAAGGCGTGCGTATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACg                                                                      <  1:384029/149‑1 (MQ=255)                                                                               aTCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATg                                                                  >  2:688943/1‑150 (MQ=255)                                                                                 ccccAGTACGGGGATTACCTCTTCAAA‑GGTGGTTAACCTTTGCTCTACGGCCAGGCGGCCTTTTTAAAAAGGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGa                                                                 <  1:210235/149‑1 (MQ=255)                                                                                     aGTACGCGAATTAACTCTTAAAG‑GGTGGTTAAGCTTTGCTCTCCGGCCAGGCGGCCGTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCg                                                           <  1:13462/151‑1 (MQ=255)                                                                                      gtccgcgGTTTACCTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGt                                                          <  1:86509/148‑1 (MQ=255)                                                                                      gTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGt                                                          >  1:58584/1‑151 (MQ=255)                                                                                           cgGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATgg                                                     >  2:344404/1‑151 (MQ=255)                                                                                           cgGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATgg                                                     >  2:700833/1‑151 (MQ=255)                                                                                            ggATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATggg                                                    >  2:141633/1‑151 (MQ=255)                                                                                               ttAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGc                                                   >  1:948074/1‑149 (MQ=255)                                                                                                 aaCTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGa                                                >  1:71275/1‑150 (MQ=255)                                                                                                         aaa‑GGTGGTTAAGCCTTGCTCTCCGGCCGGGCGGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCa                                          <  1:423611/148‑1 (MQ=255)                                                                                                         aaa‑GGTGGTTAACCTTTCCTTTCCGGCAAGGCAGCGTTTTTAAAAAGGCGTCCGATTACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAAt                                       <  1:8433/151‑1 (MQ=255)                                                                                                           aggggggtTAACCCTTGCTCTACGCCCAGCCAGCCTTTTACAAAAAGCGTACGCAGACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATa                                      <  1:436362/146‑1 (MQ=255)                                                                                                               tgtttAACCCTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATaa                                     <  2:165964/145‑1 (MQ=255)                                                                                                                gtttAAGCCTTCCTCTCCGGCCAGGCGGCCTTTTTAAAAAAGCCTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGc                                  <  1:216810/148‑1 (MQ=255)                                                                                                                 gtTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCg                                >  2:107840/1‑151 (MQ=255)                                                                                                                   tAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGt                               >  2:504498/1‑150 (MQ=255)                                                                                                                    aaGCCTTGCTCCCCGGCCAGGCAGCCGTTTTAAAAAAGCTTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTAc                             <  2:782357/151‑1 (MQ=255)                                                                                                                       ccTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGAc                          >  2:450950/1‑151 (MQ=255)                                                                                                                             tctcgggaaacgctgCCGTTTCAAAAAATCGTGCGAATAGCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACg                            <  1:945607/129‑1 (MQ=255)                                                                                                                                  ggcagggCAGCCGTTTTAAAAAAGCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAAcc                <  1:520678/145‑1 (MQ=255)                                                                                                                                    ccgggCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGt              <  1:159292/147‑1 (MQ=255)                                                                                                                                      aGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGt           >  2:292379/1‑151 (MQ=255)                                                                                                                                       ggCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGt           >  1:629607/1‑150 (MQ=255)                                                                                                                                        gCAGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTg             <  2:672665/147‑1 (MQ=255)                                                                                                                                         cAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAg         >  2:487246/1‑150 (MQ=255)                                                                                                                                         cAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGc        >  2:614458/1‑151 (MQ=255)                                                                                                                                            cctttttAAAAAGGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCa       <  2:855535/146‑1 (MQ=255)                                                                                                                                               ttttAAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGt   <  1:9229/150‑1 (MQ=255)                                                                                                                                                tttCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTa  >  2:607810/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAA‑GGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTA  >  NC_000913/115652‑115943

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