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

CTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTAC  >  NC_000913/115649‑115944                                                                                                                                                     |                                                                                                                                                    ctGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTc                                                                                                                                                     <  1:946571/149‑1 (MQ=255)   gCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCa                                                                                                                                                    <  2:266435/148‑1 (MQ=255)              cccATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCAt                                                                                                                                      >  1:143158/1‑151 (MQ=255)                  tCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATAccc                                                                                                                                  >  1:458086/1‑151 (MQ=255)                  tCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATAccc                                                                                                                                  >  2:345235/1‑151 (MQ=255)                    ttGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGc                                                                                                                                >  1:104882/1‑151 (MQ=255)                              ccGTCGCCGTTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCAGGCGGCACCCCCAGTCCGGGGATTACCTCTACAAGGGTGGTAAACCCTTGCTCTGCGGCCGGGCGGCCTTTTTAAAAAAGCGTCCCCATCCCCGCCTGTcg                                                                                                                          <  1:478533/147‑1 (MQ=255)                                gTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgtg                                                                                                                    >  1:570157/1‑151 (MQ=255)                                 tCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCACCCCCAGTGCGCGGATTAACTCTTAAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATACCCGCCTGTcgtgcgt                                                                                                                     <  1:416267/149‑1 (MQ=255)                                    ccGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGtt                                                                                                                >  1:361280/1‑151 (MQ=255)                                      gtttATGCCATGCCAGCGCCAGAGTTGCTTCCTCGCCATGCGGCACCCCCAGTCCGGGGATTAACTCTTAAAGGGTGGTTAAGCCTTCCTCTACGGCCAGGCAGCTTTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                                                                                               <  1:353796/148‑1 (MQ=255)                                        ttATCCCATGCCGGCGCCAGGGTTGCTTACTCGCCATGCGGCACCCCCAGTACGCGGTTTAACTCTAAAAAGGTGTTTAAGCCTTGCTCTCCGGCCGGGCAGCCGTTTTAAAAAGGCGTCCGCATACCCGCCTGTcgtg                                                                                                                        <  1:67250/139‑1 (MQ=255)                                           tgccatcccgGCGCCAGATTTTTTTCCTCGCCATGCGGCACCCCCAGTCCGCGGATTACCTCTTCAAAGTGGTTTAAGCTTTGTTCTACGCCCAGGCAGCTTTTTTAAAAAGGCTTCCACTTCCCCGCCTGTCGTGCGTGCGTTTCCAGc                                                                                                          <  1:518812/150‑1 (MQ=255)                                             ccatgccaGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTGCGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCg                                                                                                         <  2:546201/149‑1 (MQ=255)                                                tgccaGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTc                                                                                                     >  2:901977/1‑150 (MQ=255)                                                          agagTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAt                                                                                          >  2:157078/1‑151 (MQ=255)                                                                  ttccgcGCCAGGCGGCACCCCCAGCACGCGGTTTACCTCTTCAAGGGTGGTTAGGCCTTGCTCTCCGGCCAGGCGGCTGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCg                                                                                     <  2:267380/143‑1 (MQ=255)                                                                     cccGCCATGCGGCACCCCCAGTACGCGGATTAACTCTTAAAAGGTGGTAAAGCTTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaa                                                                                  <  2:594932/146‑1 (MQ=255)                                                                     cTCGCCATGCGGCACCCCCAGTACGGGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTAAAAAAGCGTAGGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa                                                                                 <  1:222953/149‑1 (MQ=255)                                                                      tCGCCATGCGGCACCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTCCGGCCAGGCGGCCTTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCgg                                                                                    <  1:970558/145‑1 (MQ=255)                                                                         ccATGCGCCACCCCCAGTCCGCGGATTACCTCTTAAAAGGTGTTTAAGCCTTGCTCTACGGCCAGGCAGCGTTTTTAAAAAGGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGc                                                                            <  2:915680/150‑1 (MQ=255)                                                                         cagtGCGGCCCCCCAAGTCCGCGGATTATCTCTTAAAGGGTGGTTAAGCTTTGCTCTCCGGCCAGGCCCCTGTTTTAAAAAAGCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGc                                                                            <  1:801858/147‑1 (MQ=255)                                                                             ggggCACCCCCAGTCCGGGGATTACCTCTTAAAGGGTGTTAAAGCTTTGCTCTACGGCCAGGCAGCCTTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGa                                                                         <  2:741674/147‑1 (MQ=255)                                                                             gCGGCCCCCCCAGTACGCGTATTACCTCTAAAAAGGTGTTTAAGCTTTCCTCTACGGCCAGGCCGCCGTTTTAAAAAAGCGTACCCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACg                                                                       <  1:412463/151‑1 (MQ=255)                                                                              cGGCACCCCCAGTCCGCGGATTATCTCTACAAGGGGGGTAAGGCCTTGCTCCACGCCGGGGCTGCCTTTTCCAAAAAGCGTGCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt                                                                           <  1:441727/146‑1 (MQ=255)                                                                               ggCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGa                                                                      >  2:569259/1‑150 (MQ=255)                                                                                  aTCCCCAGTACGCGGATTAACTCTTCAAAGGGGGGTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATg                                                                   >  1:538301/1‑150 (MQ=255)                                                                                    ccccAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGAcc                                                                >  1:345604/1‑151 (MQ=255)                                                                                     cccAGTACGCGGATTAACTCTGAAAGGGGGGTTAAGCTTTGCTCTGCGGCCGGGCTGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCg                                                               <  1:393394/151‑1 (MQ=255)                                                                                         gtccgcgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGc                                                             <  2:709483/146‑1 (MQ=255)                                                                                           aCGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGtt                                                          >  2:687537/1‑150 (MQ=255)                                                                                             gcgAATTACCTCTTAAAAGGTGGTTAAGCTTTGCTCTACGGCCGGGCCCCCTTTTTAAAAAAGGGTACGCATGCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCCGACGAATGACCGGCGTTATg                                                       <  2:706860/151‑1 (MQ=255)                                                                                               gtaTTAACTCTAAAAAGGGGGTTAAGCCTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATggg                                                     <  1:501193/149‑1 (MQ=255)                                                                                               ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATggg                                                     >  1:723099/1‑151 (MQ=255)                                                                                               ggATTAACTCTACAAGGGGGGTTAAGCCTTGCTCTACGCCCAGGCAGCCTTTTAAAAAAAGCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTAt                                                        <  1:713157/148‑1 (MQ=255)                                                                                                  ttAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGc                                                    >  2:475717/1‑149 (MQ=255)                                                                                                  taaaCTCTTAAAAGGTGGTAAAGCCTTGCTCTCCGCCCGGGCCGCGTTTTTAAAAAAGCTGACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGtt                                                          <  1:90917/141‑1 (MQ=255)                                                                                                   tAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGc                                                    >  1:839423/1‑148 (MQ=255)                                                                                                              aGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATa                                       >  1:417579/1‑150 (MQ=255)                                                                                                               ggtggtTAAGCTTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAg                                     <  1:546895/151‑1 (MQ=255)                                                                                                                gtggtTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAgg                                    >  2:138129/1‑151 (MQ=255)                                                                                                                  gtttAAGCCTTGCTCTACGGCCAGGCGGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAgg                                    <  1:700466/147‑1 (MQ=255)                                                                                                                  ggtTAAGCTTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGcc                                  <  2:40770/151‑1 (MQ=255)                                                                                                                   gtaaaGCCTTCCCCTCCGGGCGGGCGCCCGTTTCAAAAAAGTGAACGCAATCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATaa                                      <  1:797587/143‑1 (MQ=255)                                                                                                                       aGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACg                             >  1:728683/1‑151 (MQ=255)                                                                                                                        gCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACg                             >  2:380831/1‑150 (MQ=255)                                                                                                                        gCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACg                             >  2:468972/1‑150 (MQ=255)                                                                                                                          cTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGAcc                          >  1:240385/1‑151 (MQ=255)                                                                                                                            tGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGAca                          <  2:223511/149‑2 (MQ=255)                                                                                                                            tGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCa                         <  2:729716/150‑1 (MQ=255)                                                                                                                                caaCGGCCAGGCAGCCGTTTTCAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATaaa                     <  2:411039/148‑1 (MQ=255)                                                                                                                                  aCGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAAc                  >  1:19896/1‑151 (MQ=255)                                                                                                                                   cGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAAcc                 >  1:447691/1‑151 (MQ=255)                                                                                                                                     gCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGt               >  1:740784/1‑151 (MQ=255)                                                                                                                                      ccAGGCGCCCTTTTACAAAAAGCTGCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATaaa                     <  1:239935/144‑1 (MQ=255)                                                                                                                                        aGGCAGCCTTTTTAAAAAAGCGTACGAATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTgg             <  2:660323/150‑1 (MQ=255)                                                                                                                                          gCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGt            >  2:355408/1‑149 (MQ=255)                                                                                                                                               cGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATg     >  1:391245/1‑151 (MQ=255)                                                                                                                                                 ttttCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGt    >  2:878509/1‑150 (MQ=255)                                                                                                                                                  tttCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTac  >  2:675926/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTAC  >  NC_000913/115649‑115944

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