Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,5960GA40.0% 85.5 / 43.5 71F13F (TTC→TTT) ymgDuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (37/5);  new base A (0/28);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 8.49e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACA  >  NC_000913/1222446‑1222740
                                                                                                                                                        |                                                                                                                                                  
ggCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCg                                                                                                                                                    <  1:1119107/151‑1 (MQ=255)
  cAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGaa                                                                                                                                                  <  1:499226/151‑1 (MQ=255)
   aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGaa                                                                                                                                                  >  1:638654/1‑150 (MQ=255)
      gCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaac                                                                                                                                               >  1:354003/1‑150 (MQ=255)
       cTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaaca                                                                                                                                              >  1:604177/1‑150 (MQ=255)
           cattacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑g                                                                                                                                       <  2:735462/151‑1 (MQ=255)
              tacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGTCCAGGTCATTGGAGCGAAAACAACCA‑‑GGcc                                                                                                                                    >  2:343271/1‑151 (MQ=255)
               acccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGcc                                                                                                                                    <  1:467864/150‑1 (MQ=255)
                   attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGcc                                                                                                                                >  2:728378/1‑150 (MQ=255)
                   attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCa                                                                                                                               <  2:290471/151‑1 (MQ=255)
                   attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCa                                                                                                                               <  2:1012320/151‑1 (MQ=255)
                       ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTa                                                                                                                            >  2:675022/1‑150 (MQ=255)
                             caaTTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCaa                                                                                                                      >  2:626595/1‑150 (MQ=255)
                               aTTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAAtt                                                                                                                   >  1:65216/1‑151 (MQ=255)
                                ttccggaaaCCTGTGCGCTGTCAGTTTTTC‑CAATTTCTTTGGTTTCTGCGAAAGTAATATGGTTTTCTGCCGGAAACCCTTTTTAAA‑TTTTGAAGTCCCGCGCCCAGGTAATTGGTGCAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAAtt                                                                                                                   <  1:374096/143‑1 (MQ=255)
                                 ttCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAAtttt                                                                                                                 >  2:434735/1‑151 (MQ=255)
                                   cTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAAtttttt                                                                                                               >  2:949116/1‑151 (MQ=255)
                                       aaCCTGTGCGTTGTCATTTTTTG‑CCATTTCGTTGGTTTCTGCGCAAGAAATATAGGTTTCTGACGGCAACCCTTTTTAAA‑TATTTTAGTCCCGCCCCAAGGTCTTTGGTGAAAAAACAAACG‑‑GGCCCGCCAGTAAAGCAAAttttt                                                                                                                <  1:525121/146‑1 (MQ=255)
                                       aaCCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                           >  2:1114182/1‑151 (MQ=255)
                                         ccggtgCGCTGTCATTTTTTG‑CTTTTTCGTTGGTTTCTGCGAAAGTAATATAGTTTTCTGACGACAAACCTTTTTAAA‑TTTTGTAGTCCGGCGCCAGGGTCTTTGGTGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTtcatca                                                                                                          <  2:1226515/147‑1 (MQ=255)
                                              gctcTGTCTTTTTTTG‑CCTTTTTGTTGTTTTCTGCACAAGTAATATAGGTTTCGGACGGAAAACTTTTTTAAA‑TATTGTAGTCCTGCCCCCAGTTCATTGGTGAAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                      <  2:113275/146‑1 (MQ=255)
                                              gcgtTGTAATTTTTTC‑CCTTTTCTTTGGCTTCTGCGAAAGAAATATGGTTTTCGGACGGCAACCCTTTTTAAA‑TTTTGTAGTCCGGCCCCCAGTTCATTGGTGAAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCAtt                                                                                                     <  1:655545/150‑1 (MQ=255)
                                                 ggtgtcgTTTTTGCCCATTTCGTTGCTTTCTGCGAAATTAATATAGGTATCAGACGAAAACCCTTTTAAAA‑TATTGAAGCCCGGCGCCATGGTCATTGTAACAAAAAAAAAAA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCAtt                                                                                                     <  2:495878/142‑1 (MQ=255)
                                                  tGTCTGTTTTTG‑CCTTTTCGTTGGTTTCTGCGAAAGTAATATAGGTTTCGGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                  <  1:128875/149‑1 (MQ=255)
                                                  tGTCATTTTTTG‑CCTTTTCGTTGGTTTCTGCGCAATTAATATAGGTTTCTGCCGGCAAACCTTTTTTAA‑TTTTGTAGTCCTGCCCCAAGGTCATTGGTGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTcc                                                                                                   <  1:446402/148‑1 (MQ=255)
                                                    tCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAtt                                                                                              >  2:954853/1‑151 (MQ=255)
                                                       ttttttG‑CCGTTTGTTTGGTTTTTGCGCAGTAAATTTTTTTTTCTGGCTGCAACCTTTTTTTAA‑TTTTTCGGTCCTGCACCC‑GGTTCTTGG‑GCTATAAAAAACATGGGGCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                            <  2:1230395/149‑1 (MQ=255)
                                                        tttttG‑CTATTTCGTTGGTTTCTGCGAAATTAATATAGGTTTCTGACGGCAAACCTTTTTTAA‑TATTGGAGTCCGCCCCCAAGTTTTTTGGTGAAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCa                                                                                           <  2:478712/150‑1 (MQ=255)
                                                         ttttG‑CTTTTTCGTTGGCTTCTGCAAAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAA‑TATTGTAGTCCTGCGCCAAGGTCTTGGGTGAAAAAACAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAt                                                                                          <  2:147395/150‑1 (MQ=255)
                                                               ccATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                    >  2:697199/1‑151 (MQ=255)
                                                               ccATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                    >  2:44400/1‑151 (MQ=255)
                                                                cATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                   >  2:581841/1‑151 (MQ=255)
                                                                  tttCGTTGGTTTCTGCGCAAGTAATATAGTTTTCGGACGGAAAACTGTTTTAAA‑TTTTGTAGCCCTGCGCCCAGTTCATTGGTGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCa                                                                                  <  1:241501/150‑1 (MQ=255)
                                                                  tttCGTTGGCTTCTGCGCAAGTAATATAGGTTTCGGACGGCAACCCTTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGGGAAAAAACAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcca                                                                                   <  1:432685/149‑2 (MQ=255)
                                                                       ttGGTTTCTGCGAAAGTAATATGGTTTTCTGACGGCAAACCTTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCTTTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                             <  1:1168255/150‑1 (MQ=255)
                                                                       ttGGCTTCCGCGAAAGTAATATAGTTTTGTGACGGCAACCCTTTTTAAA‑TATTGAAGCCCGGCCCCCAGGTCATTGGTGAAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                              <  1:697208/149‑1 (MQ=255)
                                                                        tGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                            >  2:551247/1‑150 (MQ=255)
                                                                         gttttCTGCGCAAGTAATATGGGTTGCTGCCGGAAAACTGTTTTAAA‑TTTTGTAGTCCTGCGCCCAGTTTTTTGGTGCAAAAACAAACG‑‑GCCCCGCCAGTAAAGCAAATTTTTTCAACATCATTCCTCATTCAAGTTTACCCAGAATTg                                                                            <  2:1082188/146‑1 (MQ=255)
                                                                         gttttCTGCACAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAA‑TATTGTAGTCCTGCCCCCAGTTCTTGGGTGAAAAAACAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                            <  1:1010474/146‑1 (MQ=255)
                                                                         ggCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                           >  2:712158/1‑150 (MQ=255)
                                                                         ggCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                           >  1:774836/1‑150 (MQ=255)
                                                                                 cgcAAGAAATATGGTTTTTTGACGGCAACCTTTTTTAAA‑TATTTTGGTCCGCCCCCCAGTTCATGGGTGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGt                                                                  <  1:343285/151‑1 (MQ=255)
                                                                                    aaGTATTAGAGTTTTCGGCGGGCAACCTTTTTTTAA‑TTTTGTAGTCCTGCCCCAAGTTCATTGGAGAAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAc                                                                <  2:639716/150‑1 (MQ=255)
                                                                                    aaGTAATATGGGTTTCGGACGGCAACCCTTTTTTAA‑TATTGTAGCCCTGCGCCAAGTTCTTTGGTGCAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTa                                                                 <  1:85605/149‑1 (MQ=255)
                                                                                    aaGTAATATAGGTTTCGGACGGCAACCCTTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTa                                                                 <  1:812605/149‑1 (MQ=255)
                                                                                               gTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACgg                                                    >  1:21645/1‑151 (MQ=255)
                                                                                                tttCTGCCGTCAACCCTTTTTTAA‑TTTAGTAGCCCCGCCCCCAGGCTTTTGGTGAAAAAACAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACg                                                     <  1:1099375/149‑1 (MQ=255)
                                                                                                 ttCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCTGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCCCCCGTACCGCCAATCACGGc                                                   >  2:290092/1‑150 (MQ=255)
                                                                                                     gACGGCAACCTGTTTTAAA‑TTTTGTAGTCCTGCGCCCAGTTCTTGGGTGAAAAAACAAACG‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCtt                                              <  1:1074459/151‑1 (MQ=255)
                                                                                                      aCGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                             >  1:963961/1‑151 (MQ=255)
                                                                                                         ggaaaCCTTTTTTAATTTTTGTAGCCCTGCCCCCAGTTCTTGGGGGAAAAAAAAACCG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTa                                            <  1:84334/148‑1 (MQ=255)
                                                                                                          cAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAAt                                          >  1:953735/1‑150 (MQ=255)
                                                                                                           aacctttttttAA‑TATTGTACCCCGGCCCCCAGGTCTTTGGTGAAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCAc                                                 <  1:1186200/136‑1 (MQ=255)
                                                                                                             aCCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGcc                                      >  1:678665/1‑151 (MQ=255)
                                                                                                              ccGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGccc                                     >  1:1019000/1‑151 (MQ=255)
                                                                                                                gTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCt                                   >  2:95991/1‑151 (MQ=255)
                                                                                                                 tttttAA‑TTTTGTAGCCCGGCCCCCAGGTCATTGGTGCAAAAAAAAACG‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCt                                   <  1:338378/150‑1 (MQ=255)
                                                                                                                 tttttAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCtt                                  >  2:482875/1‑151 (MQ=255)
                                                                                                                  tttaAA‑TATTGTGGTCTTCCCCCCGGGTTATTGGTCAAAAAAAAAACG‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCtt                                  <  1:1233165/150‑1 (MQ=255)
                                                                                                                   tttAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                >  1:1218155/1‑151 (MQ=255)
                                                                                                                     tAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                >  2:6971/1‑149 (MQ=255)
                                                                                                                      aa‑TATTGTACCCCGCCCCCCAGCTCTTTGGCGAAAAAAAAAAGG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAg                              <  1:728385/150‑1 (MQ=255)
                                                                                                                                  ccTACGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCg                    >  2:269143/1‑149 (MQ=255)
                                                                                                                                   cTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                 >  1:657740/1‑151 (MQ=255)
                                                                                                                                     gcgcCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                 >  1:117707/1‑149 (MQ=255)
                                                                                                                                       gcCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa              >  2:646668/1‑150 (MQ=255)
                                                                                                                                         ccAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACg           >  2:968537/1‑151 (MQ=255)
                                                                                                                                              tCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCATGGCACCTTTAATCGCCCCTTCGAGGCCATTCCCGGTCAGTACGACa        >  2:447896/1‑149 (MQ=255)
                                                                                                                                                aTTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg     >  1:895408/1‑150 (MQ=255)
                                                                                                                                                 ttGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc   >  1:640057/1‑151 (MQ=255)
                                                                                                                                                  tGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc   >  1:292142/1‑150 (MQ=255)
                                                                                                                                                   ggTGCAAAAACAACCA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACa  <  1:832921/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTG‑CCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAA‑TATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACA  >  NC_000913/1222446‑1222740

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