Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,724 A→C 15.3% N105T (AAC→ACC)  gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,7240AC15.3% 257.6 / 42.2 111N105T (AAC→ACC) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (56/38);  new base C (7/10);  total (63/48)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

Only 100 of 115 total aligned reads displayed.

AAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGAT  >  NC_000913/869625‑869823
                                                                                                   |                                                                                                   
aaaaCGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCaa                                                                                                     >  1:643874/1‑100 (MQ=255)
aaaaCGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCaa                                                                                                     >  1:505335/1‑100 (MQ=255)
      tGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcgccgc                                                                                               <  1:2625729/100‑1 (MQ=255)
        cTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGACATTAAATTCCAGGATGGCACCGATTTCAACGCcgcg                                                                                              >  1:2676843/1‑99 (MQ=255)
         tGGCGGAGAGTTATACTGTTTCCGATGACGGCATTACTTACACCGGGAAATTGCGGGAAGGCATTAAATTCCATGATGGCACCGATTTCAACGCCgcggc                                                                                            <  2:1491093/100‑1 (MQ=255)
         tGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCgcggc                                                                                            >  2:2189241/1‑100 (MQ=255)
          ggCGGAGAGTTGTACCGTTTCCGATAAAGGCATTACTTACACGGTGAAATAGCGGGAAGGCATTAAAGTCCAGGAAGGAACCGATTTCAAc                                                                                                    >  2:1580616/1‑91 (MQ=255)
          ggCGGAGAGTTATAACGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAAGCATTAAATTCCAGGATGGCACCGATGTCAACAACgcggcg                                                                                           >  2:525631/1‑100 (MQ=255)
           gCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCgcggcgt                                                                                          <  1:1541788/100‑2 (MQ=255)
           gCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGcggcgg                                                                                          <  1:1969344/100‑1 (MQ=255)
            cGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGt                                                                                         >  1:1233545/1‑100 (MQ=255)
            cGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGt                                                                                         >  1:2990992/1‑100 (MQ=255)
            cGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGt                                                                                         >  1:2472414/1‑100 (MQ=255)
             ggAGAGTTATACCGTTTCCGATTACGGCATTACTTATACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTg                                                                                        <  1:2563995/100‑1 (MQ=255)
                  gTTATGCCGTTTCCGATGACGGCATTACTGACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGTTTTCAACGCCGCGGCGGTGAAAg                                                                                    <  2:995214/99‑1 (MQ=255)
                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGaaatc                                                                                   <  1:1462852/100‑3 (MQ=255)
                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGaaatc                                                                                   <  1:1241001/100‑3 (MQ=255)
                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGaaatc                                                                                   <  2:2483356/100‑3 (MQ=255)
                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGc                                                                                   <  1:2218612/100‑1 (MQ=255)
                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGc                                                                                   <  1:2780418/100‑1 (MQ=255)
                    tataCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGAAATCgg                                                                                 <  2:1748989/100‑2 (MQ=39)
                    tataCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGAAATCgg                                                                                 <  2:610794/100‑2 (MQ=39)
                     ataCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGaa                                                                                >  1:1781126/1‑100 (MQ=255)
                      taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAATGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGaaatcggtg                                                                               <  1:1475999/100‑7 (MQ=38)
                      taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGCTTTCACCGCCGCGGCGTTGAAATCggt                                                                                <  2:1940552/99‑3 (MQ=38)
                      taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGaaatcggtg                                                                               <  2:879826/100‑7 (MQ=38)
                      taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAAt                                                                               <  2:641185/100‑1 (MQ=255)
                       aCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATc                                                                              <  2:2103807/100‑1 (MQ=255)
                        ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATc                                                                              >  2:200050/1‑99 (MQ=255)
                        ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCt                                                                             >  1:1791830/1‑100 (MQ=255)
                        ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCt                                                                             >  1:1527601/1‑100 (MQ=255)
                             tCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGAcc                                                                        <  2:2361192/100‑1 (MQ=255)
                             tCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGAcc                                                                        <  2:2251276/100‑1 (MQ=255)
                             tCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGAcc                                                                        <  2:687465/100‑1 (MQ=255)
                               cGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCgg                                                                      >  2:1327655/1‑100 (MQ=255)
                               cGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGAGGCGGTGAAAGCGAATCTAGACCg                                                                       >  2:870804/1‑99 (MQ=255)
                                   gACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGcc                                                                   <  2:2174167/99‑1 (MQ=255)
                                     cGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGc                                                                <  2:1997787/100‑1 (MQ=255)
                                      ggCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCg                                                               <  1:2803611/100‑1 (MQ=255)
                                      ggCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCg                                                               >  2:115994/1‑100 (MQ=255)
                                      ggCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATATGGACCGGGCCa                                                                  >  1:2948643/1‑97 (MQ=255)
                                       gCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGa                                                              >  1:651335/1‑100 (MQ=255)
                                       gCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGCGAAAGCGAATCTGGACCGGGCCAGCGa                                                              >  1:1530400/1‑100 (MQ=255)
                                        cATTACTTACGCCGTGAAGTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGAt                                                             <  2:473120/100‑1 (MQ=255)
                                        cATTACTTACACCGTTAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACACCGCGGCGGTGAAAGAGAATCTAGACCGGGCCAGCGAt                                                             >  1:1314256/1‑100 (MQ=255)
                                        cATTACTTACACAGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGc                                                                <  1:728668/97‑1 (MQ=255)
                                           tacttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCg                                                          >  1:641547/1‑100 (MQ=255)
                                            acttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCgg                                                         >  1:2616273/1‑100 (MQ=255)
                                             cttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGTATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGc                                                        <  1:527555/100‑1 (MQ=255)
                                             cttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGc                                                        >  1:748973/1‑100 (MQ=255)
                                              ttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCg                                                       <  1:3044205/100‑1 (MQ=255)
                                               tacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGa                                                      >  1:763932/1‑100 (MQ=255)
                                               tacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGa                                                      >  2:2856041/1‑100 (MQ=255)
                                               tacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGAGGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCg                                                       <  2:1218337/99‑1 (MQ=255)
                                                 cacCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAat                                                    >  2:2706950/1‑100 (MQ=255)
                                                   ccGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAatca                                                  >  1:1596483/1‑100 (MQ=255)
                                                   ccGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGAACCGGCGAatca                                                  >  2:1898336/1‑100 (MQ=255)
                                                     gTGAAATTGCGGGAAGGCATTAAATTCCTGGTTGGTACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAatcatc                                                <  1:1197506/100‑1 (MQ=255)
                                                      tGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCt                                               >  1:2672152/1‑100 (MQ=255)
                                                        aaaTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTa                                             <  1:1873488/100‑1 (MQ=255)
                                                         aaTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGAACCGGGCCAGCGATCCGGCGAATCATCTTaa                                            >  1:406749/1‑100 (MQ=255)
                                                         aaTTGCGGGAAGGCATTAAATTCCAGGATAGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTaa                                            >  2:746614/1‑100 (MQ=255)
                                                          aTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTaaa                                           >  2:1161313/1‑100 (MQ=255)
                                                            tgggggAATGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACg                                         <  1:2770913/97‑1 (MQ=255)
                                                            tGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACAGGGCCAGCGATCCGGCGAATCATCTTAAACg                                         >  1:2701300/1‑100 (MQ=255)
                                                              cGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCt                                       >  2:574488/1‑100 (MQ=255)
                                                                   aGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAAc                                  >  1:1277034/1‑100 (MQ=255)
                                                                    ggCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGTCCAGCGATTCGGCGAATCATCTTAAACGCTATAAcc                                 <  2:1483102/100‑1 (MQ=255)
                                                                    ggCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAAcc                                 <  1:2706982/100‑1 (MQ=255)
                                                                      cATTAAATTCCAGGATGACACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTg                               >  2:1516006/1‑100 (MQ=255)
                                                                       aTTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGt                              >  2:261228/1‑100 (MQ=255)
                                                                        ttAAATTCCATGTTGACACCTATTTCAACGCCGCGGCGGTGAATGCTAATCTGAACCGGGCCATCGATCAGTCGAATCATCTTAAACGCTATAACGTGta                             <  2:1113253/100‑1 (MQ=255)
                                                                          aaaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGCAAGCGAATCTGGACCGGGCCAGCGAACCGACGAATCATCTTAAACGCTAAAACCtatata                           >  1:518856/1‑100 (MQ=255)
                                                                          aaaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGtata                           >  1:518850/1‑100 (MQ=255)
                                                                          aaaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGAACCGGCGAATCATCTTAAACGCAATAACCTGtata                           >  1:518855/1‑100 (MQ=255)
                                                                          aaaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAAAGAACCGGCGAATCATCTTAAACGATATAACCTGtata                           >  1:518840/1‑100 (MQ=255)
                                                                            aTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAg                         <  2:763928/100‑1 (MQ=255)
                                                                            aTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAg                         >  1:172951/1‑100 (MQ=255)
                                                                                cAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGACTCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAata                     >  1:1406540/1‑100 (MQ=255)
                                                                                    ccagattcgctttCACCGCCGCGGCGTTGAAATCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGc                 >  2:1115214/11‑100 (MQ=25)
                                                                                    aTGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGc                 <  1:2978201/100‑1 (MQ=255)
                                                                                     tGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCt                >  1:416908/1‑100 (MQ=255)
                                                                                       gCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCAAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaa              >  1:1350748/1‑100 (MQ=255)
                                                                                       gCACCGATTTCAACGCCGCGGCGGTGAAAGAGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaa              >  2:1911652/1‑100 (MQ=255)
                                                                                        cACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaa             >  2:1320014/1‑100 (MQ=255)
                                                                                         tctgctttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaa            <  1:2008729/95‑1 (MQ=255)
                                                                                         aCCGATTTCAACTCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaa            <  2:739161/100‑1 (MQ=255)
                                                                                           cgctttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACg          >  1:2404184/4‑100 (MQ=255)
                                                                                           cgctttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACACTATAACCTGTATAAGAATATTGCTAAAACg          >  1:683856/4‑100 (MQ=255)
                                                                                            gctttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaaaga         >  2:2008706/3‑97 (MQ=255)
                                                                                            gctttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACgg         >  2:2456621/3‑100 (MQ=255)
                                                                                            gctttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACgg         >  2:1305781/3‑100 (MQ=255)
                                                                                            gATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACgg         >  1:504328/1‑100 (MQ=255)
                                                                                             aTTTCAACGCCGCGGCGGTGAATGCGAATTTGGACTGGGTCAGCGATTCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGa        <  1:1533557/100‑1 (MQ=255)
                                                                                              tttCAACGCCGCGGCGGTAAAAGCGAATCTGGGCCTGGCCAGCGATCCGACGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaa             <  2:2569662/94‑1 (MQ=255)
                                                                                               ttCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAg      >  2:2245146/1‑100 (MQ=255)
                                                                                                 cAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcg    >  1:2571764/1‑100 (MQ=255)
                                                                                                 cAACGCCGCGGCGGTGAAAGCGAATCTGGAACGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGc     >  2:1913691/1‑99 (MQ=255)
                                                                                                  aaCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcga   >  2:511349/1‑100 (MQ=255)
                                                                                                   aCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgat  >  2:2209123/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGAT  >  NC_000913/869625‑869823

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: