Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,393,679	0	A	T	38.1%	87.9 / 23.3	63	S860R (AGT→AGA)	yhdP	DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (13/26);  new base T (24/0);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 2.61e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAACGTTATGTGTTGTGGGAAACTCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCAC  >  NC_000913/3393534‑3393827                                                                                                                                                  |                                                                                                                                                        gtaACGTTATGTGTTGTGGGAAACTCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGGGTTTTTCTg                                                                                                                                                      >  1:569928/1‑149 (MQ=255)   aaCGTTATGTGTTGTGGGAAACTCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGGTTTTCTGt                                                                                                                                                     >  2:759651/1‑148 (MQ=255)   aaCGTTATGTGGTGTGGGAAACTCGCTGCCCCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCt                                                                                                                                                   <  2:413119/150‑1 (MQ=255)        tATGTGTTGTGGGAAACTCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGGTTTTCTGTCTg‑ccgcc                                                                                                                                            >  1:500015/1‑151 (MQ=255)         aTGTGTTGTGGGAAACTCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTg‑ccgcc                                                                                                                                            <  2:604113/150‑1 (MQ=255)               tgtgGGAAACTCGCTGCACCACCGACCCTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAAt                                                                                                                                       <  2:36929/149‑1 (MQ=255)               tgtgGGAAACTCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGGTTTTCTTTCTG‑CCCCCCcaa                                                                                                                                        >  1:805999/1‑148 (MQ=255)                 cacttaaCTCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGTTTTACTGTCTC‑CCCCCCaaaaaac                                                                                                                                    >  2:753841/6‑146 (MQ=255)                         tCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGCGCGTTTTTCTGTCTT‑CCGCCCCAATa                                                                                                                                      >  2:840059/1‑140 (MQ=255)                             tGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGTGGTTTTCTGTGCG‑CCGCCCCAATAGCACgcaccaacggc                                                                                                                       >  2:224830/1‑141 (MQ=255)                                accaccGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGGTTTACTGTCTG‑CCGCCCCAATAGCACGAACGAGGGGCAGc                                                                                                                    >  1:885768/1‑151 (MQ=255)                                 ccaccGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGc                                                                                                                    <  2:437775/150‑1 (MQ=255)                                 ccaccGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGAGGTTTTTAGGCTGCCCGCCCAAA‑AGGACCagcgagcggccg                                                                                                                     >  1:710239/1‑145 (MQ=255)                                      gACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCCCCCCAAAAGCACGATCGAGCGTCAGCtttt                                                                                                                >  1:252704/1‑149 (MQ=255)                                       acacTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTT‑CCGCCCCAAAAGCACCAACGCGCGTCAGCtttt                                                                                                                >  2:744219/1‑148 (MQ=255)                                            ctcCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGGTTTTCTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTTACCGAg                                                                                                         >  1:771532/1‑150 (MQ=255)                                              ccGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGGTTTTCTGTCTG‑CCCCCCCAATAACACGAACGAGCGGCAGCCTTTTAcc                                                                                                            >  2:667955/1‑145 (MQ=255)                                              ccGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGc                                                                                                        <  1:367690/149‑1 (MQ=255)                                               cGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTTCTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑aa                                                                                                     >  1:865317/1‑150 (MQ=255)                                                gccgccccTTGCTGGAACAGTGCCGGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGATCGTTTTACTGTCTG‑CCGCCCAAGTAGCACGATCGAGCGTCAGCTTTTGACCGAgg                                                                                                        <  2:66673/141‑2 (MQ=255)                                                 cAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCCAATAGCACGATCGAGCGTCAGCTTTTTACCCAGC‑aaac                                                                                                   >  1:404495/1‑148 (MQ=255)                                                      ccTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGGTTTTCTGTCTG‑CCGCCCAAATAGCACGACCGAGCGTCAGCTTTTTACCGAGC‑AACCAGCgg                                                                                              >  2:696949/1‑150 (MQ=255)                                                         ttCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTGGACCGAGC‑AACCAGCGGCt                                                                                            <  2:313662/149‑1 (MQ=255)                                                           cTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGGTTTTCTGTCTC‑CCGCCCCAATAGCACGATCCAGCGGCAACTTTTGACCGAGC‑AACCAGCGGCTAtt                                                                                         >  1:157388/1‑150 (MQ=255)                                                             ggAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGa                                                                                       <  1:581894/150‑1 (MQ=255)                                                              gAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGaaa                                                                                     <  1:552601/151‑1 (MQ=255)                                                               aaCAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAAt                                                                                    <  1:454556/151‑1 (MQ=255)                                                                 cAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTTCTTTCTG‑CCCCCCAAATAACACGATCCAGCGTCAGCTTTTGACCCAGC‑AACCAGCGGGCATTTaaa                                                                                     >  1:484667/1‑148 (MQ=255)                                                                  aGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTTCTGTCTG‑CCGCCCAAATAGCACGATCGAGCGGCAGCTTTTTACCCAGC‑AACCAGCGGCTATTTaaa                                                                                     >  1:78920/1‑147 (MQ=255)                                                                    ggCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGGTTTTCTGTCTG‑CCGCCCCAAAAGCACGATCGGGCGGCAGCCTTTTACCGCGC‑AACCAGCGGCCATTGAAATGAtt                                                                                >  2:794737/1‑150 (MQ=255)                                                                       cagccaCTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATcc                                                                            <  1:243386/151‑1 (MQ=255)                                                                       cagccaCTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATcc                                                                            <  2:830577/151‑1 (MQ=255)                                                                        agccaCTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGATTTTATTTTTT‑CCCCCCCAATAGCACGATGGAGGGTCTGCGTTTTGTCGAGC‑GATCTGCGGCCAttcacaagaatttcc                                                                            >  2:851616/1‑137 (MQ=255)                                                                            aCTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCAcc                                                                       <  1:207537/151‑1 (MQ=255)                                                                              tCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGGGGTTTTCTGCTTC‑ACGCCCAAATAGCACCAACGAGCGTCAGCTTTTTATCGAGG‑AACCCACGGCTATTGAAATGGTTAtttccacc                                                                       >  2:847886/1‑142 (MQ=255)                                                                                ggCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTTCTGTCTG‑CCGCCCCAATAGCACGATCGAGCGTCAGCTTTTTACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGcc                                                                   >  1:830538/1‑151 (MQ=255)                                                                                    ccATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTTCTGTCTG‑CCGCCCCAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAAAGATTATCCCCACCAGCCC‑GTc                                                              >  2:193505/1‑151 (MQ=255)                                                                                         cATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGtt                                                         <  1:904685/151‑1 (MQ=255)                                                                                          aTCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGt                                                          <  1:478191/149‑1 (MQ=255)                                                                                             ggcggcATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGGTTTTCTGTCTT‑CCGCCCCAAAAGCCAGAAAGGGGGGCAGGCTTTGACCGAGC‑AAACCACGGGTAATGGAATAATTCTTCGCCCCAACCTGGCCCGGTTatt                                                      >  2:609282/1‑149 (MQ=255)                                                                                             ggcggcATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGGTTTACAGTCTC‑CCCCCCAAATAAGCCGCGCGGGAGGCAGCTTTTGTCCGAGCAAACCAG‑GGCAATTGAAAAGATTATCCGCCCCAGcccctcccggtaa                                                        >  1:699431/1‑138 (MQ=255)                                                                                                ggcATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCaa                                                  <  1:403221/151‑1 (MQ=255)                                                                                                      ttGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGCGCGCTTTTATGTCTG‑CCGCCCCAATAGCCCGCACCAGCGGCCGCTTTTGACCGCGC‑ACACCGCGGCCATTTAAACGATTCACCCCACCCGCCT‑GTCCGGTTAATTCCAAGCt                                              >  2:209914/1‑149 (MQ=255)                                                                                                       tGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTTACCGAGC‑AACCAGCGGCCATTGAAATGATTATCCGCACCAGCCC‑GTCCGGGTAATTCAAAGCTGtt                                           >  1:666855/1‑151 (MQ=255)                                                                                                          gTTCAACGCCACTTTGTTCCGGCGACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCGCCAGCCT‑GTCCGGTTAATTCAAAGCTGTTga                                         <  2:57616/150‑1 (MQ=255)                                                                                                            tCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTTCTTTCTG‑CCGCCCAAATAACACGATCGAGCGTCAGCTTTTTACCGAGC‑AACCAGCGGCTATTTAAAAGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGtt                                           >  1:511232/1‑146 (MQ=255)                                                                                                             cAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTT‑CCGCCCCAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCCGCGGCTATTGAAATGATTATCCGCCCCAGCCT‑CTCCGGTTAATTCAAAGCTGGTGGGAtt                                     >  2:43281/1‑151 (MQ=255)                                                                                                             cAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTgaga                                       <  2:207221/149‑1 (MQ=255)                                                                                                              aaCGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTT‑CCGCCCCAATAACACGGTCGAGCGTCAGCTTTTTACCGAGC‑AACCAGCAGCTTTTGAAATGATTATCCGCCCCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGGGAtt                                     >  1:919871/1‑150 (MQ=255)                                                                                                                    aCTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCa                              <  2:634699/151‑1 (MQ=255)                                                                                                                    aCTTTGTTCCGGCAACGGCGGGAGCGTTTTAATGTCTG‑CCCCCCAAATAGCAACATCGAGCGGCCGCTTTTTACCGCGC‑AACCAGCGGCTATTGAAAAGATTTTCCGCACCCGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCCTc                               >  1:581705/1‑150 (MQ=255)                                                                                                                     cTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCa                              <  2:110686/150‑1 (MQ=255)                                                                                                                      tttGTTCCGGCAACGGCGGGAGCGTTTTTCTGTCTT‑CCGCCCCAATAGCACGAACGAGCGGCAGCTTTTTACCGAGC‑AACCAGCGGCCATTGAAATGATTATCCGCACCAGCCC‑GTCCGGGTAATTTAAAACTGTTGAGATTGCCATCa                              >  2:629089/1‑149 (MQ=255)                                                                                                                        tGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAAcc                           <  2:859587/150‑1 (MQ=255)                                                                                                                        tGTTCCGGCAACGGCGGGAGCGGTTTTTTGTCTG‑CCCCCCCAATAGCACGATCGGGGGTCTGCTTTTTACCGAGC‑AAACCCCGGCTATTGGAAAGATTAACCGCCCCAACCC‑GTCCCGGTAATTTAAAACTGTTGAGATTGCCATCaaac                           >  1:353495/1‑148 (MQ=255)                                                                                                                             cGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACg                     <  2:384426/151‑1 (MQ=255)                                                                                                                               gCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGtt                   <  1:95226/151‑1 (MQ=255)                                                                                                                                    ggcggGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCg              <  2:348481/151‑1 (MQ=255)                                                                                                                                       gggAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTa           <  1:841566/151‑1 (MQ=255)                                                                                                                                        ggAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCt            >  1:656243/1‑149 (MQ=255)                                                                                                                                           gCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTgg       <  1:760084/151‑1 (MQ=255)                                                                                                                                               tttACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTc    >  2:35961/1‑150 (MQ=255)                                                                                                                                                ttACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCAc  >  2:741245/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                        GTAACGTTATGTGTTGTGGGAAACTCGCTGCACCACCGACACTCTCCGCAGCGCCTTTCTGGAACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTG‑CCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGC‑AACCAGCGGCTATTGAAATGATTATCCGCACCAGCCT‑GTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCAC  >  NC_000913/3393534‑3393827

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