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

CTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGT  >  NC_000913/115647‑115942                                                                                                                                                        |                                                                                                                                                   ctctGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTc                                                                                                                                                    >  2:1203893/1‑151 (MQ=255)     gCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCaaaa                                                                                                                                                <  2:315241/151‑1 (MQ=255)       ccATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAg                                                                                                                                              <  1:1002224/151‑1 (MQ=255)         aTAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCg                                                                                                                                            <  1:672160/151‑1 (MQ=255)             cATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑Ta                                                                                                                                         <  1:1192762/149‑1 (MQ=255)             cATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACg                                                                                                                                       >  2:856171/1‑151 (MQ=255)                cccATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAATGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCa                                                                                                                                     >  1:1239153/1‑150 (MQ=255)                cccATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCAt                                                                                                                                    >  1:904099/1‑151 (MQ=255)                  cATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATAc                                                                                                                                  >  2:335371/1‑151 (MQ=255)                        gcgcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGcc                                                                                                                             >  2:96641/1‑150 (MQ=255)                        gcgcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCt                                                                                                                            >  1:1267667/1‑151 (MQ=255)                         cgcATTGCCGTCGCCGTTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCt                                                                                                                            <  1:681500/150‑1 (MQ=255)                         cgcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTg                                                                                                                           >  1:1119078/1‑151 (MQ=255)                                ccgccgccgTTTATGCCATGCCGGCGCAAGATTTTCTTACTCGCCATGCGCCACCCCCAGTACGCGGATTAA‑CTCTCAAAGGGTGTTAAAGCCTTCCTCTCCGGCCGGGCGGCCTTTTTAAAAAAGCG‑TACGCATACCCGCCTGTcgtgcg                                                                                                                    <  1:716897/151‑1 (MQ=255)                                ccGTCGCGGTTTATCCCATGCCGCCGCCAGATTTTCTTCCCCCCCAGGCGCCACCCCCAGTACGCGTTTTAA‑CTCTTAAAGGGTGTTTACCCCTTCCCCTCCGCCCGGGCGGCTGTTTTAAAAAAGCG‑TCCGCATACCCGCCTGTcgtgc                                                                                                                     <  1:631911/150‑1 (MQ=255)                                ccGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTcgtgcg                                                                                                                    >  1:632387/1‑151 (MQ=255)                                  gTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTcgtgcgtg                                                                                                                  >  2:606116/1‑151 (MQ=255)                                      ccgtttagcccatcccaGCGCGAGATTTCCTTCCTCGCCATGCGCCACCCCCATGCCGCGGATTAC‑CTCTACAAGGGTGGTTAACCCTTCCTCTCCGGCCAGGCGGCTGTTTAAAAAAGGCT‑TACGCATCCCCGCCTGTCGTgcgtgcgt                                                                                                               <  1:1110865/141‑1 (MQ=255)                                      ccGGTTATGCCATGCCAGCCCCAGATTTGCTTACTCGCCATGCGGCACCCCCAGTCCGCGGATTAA‑CTCTTAAAAGGTGGTTAAGCCTGGCTCTACGGCAGGGCAGCCGTTTTAAAAAAGCG‑TACGCATACCCGCCTGTCgtgcgtgc                                                                                                                 <  2:1317248/148‑1 (MQ=255)                                            atgccatgccAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGc                                                                                                        >  2:106685/1‑151 (MQ=255)                                                 atgccaGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAtt                                                                                                    >  1:151225/1‑150 (MQ=255)                                                    cccgcgcCAGGGGTTCTTACGCCCTAGGCGGCTCCCCGAGTGCGGGTATTAA‑CTCTGCAGGGGTGTTAAGGCCTGGCCCTCGGGCCGGGCAGCCTTTTCGAAAAAGCG‑TACGCATCGCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAAc                                                                                                <  1:553737/148‑1 (MQ=255)                                                        cgcCAAAGTTTTTTACTCCCCGTGCGGCATCCCCAGTCCGCGTATAAA‑CTCTACAAAGGTGGTTAAGCTTTGCTCTACGCCCAGGCGGCTTTTTAAAAAAAGCG‑TCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACg                                                                                               <  2:662495/148‑1 (MQ=255)                                                         gcCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCgg                                                                                           >  1:1020699/1‑151 (MQ=255)                                                          ccAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGt                                                                                          >  2:1298928/1‑151 (MQ=255)                                                          ccAGAGTTGCTTACTCCCCAGGGGGCACCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACCCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCg                                                                                            <  2:1179586/149‑1 (MQ=255)                                                           cAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGt                                                                                          >  1:436700/1‑150 (MQ=255)                                                           cAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGt                                                                                          >  2:580930/1‑150 (MQ=255)                                                                tttgttccTCGCAATGCGGAACCCCCAGTGCTCGAATAAC‑CTTTTAAACGGTGGTAAAGCCTTGCCCGAGTGCCGGGTGCCAATTTTAAAAAAGCG‑TGCGCATTTCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAt                                                                                        <  1:1254660/140‑1 (MQ=255)                                                                ttGCTAACTCCCCATGCGCCCTCCCCAGCCCGCGGATTAA‑CTCTTCAAAGGTGTTTATGCCTTGCTCTCCGGCCAGGCACCCTTTTACAAAAAGCG‑TCCGTATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTa                                                                                      <  2:1115797/149‑1 (MQ=255)                                                                    ttccccgccgggCGACACCCCCAGTCCG‑GGTTTAACCTCTTCAAAGGTGGTTAAGCTTTTCGCTACGGCAGGGCACCCTTTTAAAAAAAGCG‑TACCCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGc                                                                                    <  2:790690/136‑1 (MQ=255)                                                                    ttACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGa                                                                                 >  1:283962/1‑150 (MQ=255)                                                                       cTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAAt                                                                              >  1:292415/1‑150 (MQ=255)                                                                         cGCCATGCGGCCTCCCCAGTCCGCGGATTAC‑CTCTTAAAAGGTGGTAAAGCCTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCG‑TACCCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa                                                                               <  2:84286/147‑1 (MQ=255)                                                                            cATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTg                                                                        >  1:743192/1‑151 (MQ=255)                                                                             aTGCGGCACCCCCAGTCCGGGAATTAA‑CTCTTCAAAGGGGTTTAAGCCTTGCTCTCCGGCCAGGCGGCCTTTTTAAAAAAGCG‑TCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGa                                                                       <  2:42776/151‑1 (MQ=255)                                                                              tGCGTCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGAc                                                                      >  2:1028328/1‑151 (MQ=255)                                                                              tGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGa                                                                       >  2:820463/1‑150 (MQ=255)                                                                                      ccccAGTCCGGGGATAAC‑CTCTTCAAAGGTGTTAAACCCGTGCTCCCCGGCGCGGCAGCCGTTTAAAAAAAGCA‑TACTCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGa                                                                <  2:564255/149‑1 (MQ=255)                                                                                      ccccAGTCCGCGGATTAA‑CTCTTAAAGGGTGTTTAAGCCTTGCTCTCCGGCCAGGCAGCCTTTTTCAAAAAGCG‑TCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGAcc                                                              <  1:312701/151‑1 (MQ=255)                                                                                        cacgTACGGGAATTAC‑CCTTTAAAATGTGTTAAAGCCTTCCTCTACGGCCGGGCGGCTTTTTCCAAAAACCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGAcc                                                              <  1:144267/146‑1 (MQ=255)                                                                                          agcaCGGGGATTAT‑CTCTAAAAAGGTTTTAAAGCCTTCCTCTCCGGCCAGGCAGCC‑TTTTAAAAAAGGGATACGAACACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGc                                                           <  1:886330/147‑1 (MQ=255)                                                                                             aCGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGtt                                                        >  2:1073603/1‑150 (MQ=255)                                                                                                 ggATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTCCGGCCGGGCAGCCGTTTTAAAAAAGCG‑TACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGt                                                         <  1:361583/145‑1 (MQ=255)                                                                                                 ggATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATggg                                                   >  1:431863/1‑151 (MQ=255)                                                                                                  gATTAACCTCTTAAAGGGTGTTTAACCCTTGCTCTCCGCCCAGGCGGCCTTTTTAAAAAAGCG‑TACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATg                                                     <  1:598626/149‑1 (MQ=255)                                                                                                   aTTAC‑CTCTAAAAAGGTGGTAACGCCTTCCTCTGCGGCCGGGCAGCCGTTTAAAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATgg                                                    <  1:233783/148‑1 (MQ=255)                                                                                                       acctctTAAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCG‑TCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAg                                                <  2:221266/147‑1 (MQ=255)                                                                                                            taaaagggtgttAAAGCCTTCCTCTCCGCCCGGGCAGCCTTTTAAAAAAAGGT‑TACCCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACt                                            <  1:750298/142‑1 (MQ=255)                                                                                                              cAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCaa                                        >  2:503668/1‑150 (MQ=255)                                                                                                               aagggtgttTAGGCTTTCCTCTACGGCCGGGCAGCCTTTTTAAAAAGGCG‑TCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTc                                          <  2:688340/144‑1 (MQ=255)                                                                                                               aaggggggtAAAGCCTTGCTCTCCGGCCAGGCGCCGGTTTTAAAAAAGCG‑AACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAAt                                      <  1:556302/145‑1 (MQ=255)                                                                                                               aaaGGTGGTTAAGCCTTGCTCTCCGGCCAGGCAGCCGTTTTAAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTc                                          <  1:199294/147‑1 (MQ=255)                                                                                                               aaaGGTGGTTAAGCCTTGCTCTACGGCCGGGCGGCCGTTTTCAAAAAGCG‑TACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAAt                                      <  2:566889/151‑1 (MQ=255)                                                                                                               aaaGGTGGTTAAGCCTTGCTCTACGCCCAGGCGGCTTTTTTCAAAAAGCG‑TACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACtt                                           <  2:868053/146‑1 (MQ=255)                                                                                                                    tgtttAACCTTTGCTCTACGCCCGGGCAGCCGTTTCAAAAAAGCG‑AACGCATCTCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATaa                                    <  1:989651/145‑1 (MQ=255)                                                                                                                      gtTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCg                               >  2:565702/1‑151 (MQ=255)                                                                                                                         aaGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTAc                            >  1:912887/1‑151 (MQ=255)                                                                                                                                cctcaacgcgcAGGCAGCTTTTTAAAAAAACCG‑TCCGCAAGCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCa                       <  1:173747/139‑1 (MQ=255)                                                                                                                                      cGGCCGGGCGCCCTTTTAAAAAAATCG‑TCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATaaaaa                 <  1:815527/149‑1 (MQ=255)                                                                                                                                       gcccgggcggcCTTTTTAAAAAAGCT‑TACCTACCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATaaaaa                 <  1:670432/143‑1 (MQ=255)                                                                                                                                         ccgggCAGCCTTTTTAAAAAGGCG‑TCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTg            <  2:359702/148‑1 (MQ=255)                                                                                                                                            ggCAGCTTTTTAAAAAAAGGG‑ACCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTg            <  1:1160999/148‑1 (MQ=255)                                                                                                                                             gcgcccGTTTTAAAAAAGCG‑ACCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAAcc               <  1:348537/140‑1 (MQ=255)                                                                                                                                                gCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCa      >  1:1036918/1‑150 (MQ=255)                                                                                                                                                  cGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAAt    >  2:1003423/1‑150 (MQ=255)                                                                                                                                                    ttttCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGt  >  1:561645/1‑150 (MQ=255)                                                                                                                                                    ttttAAAAAGGCG‑TCCCGATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATaaa                   <  1:223546/133‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   CTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAA‑CTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCG‑TACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGT  >  NC_000913/115647‑115942

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