Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,795,591 G→A 21.8% S221S (TCG→TCA waaF → ADP‑heptose:LPS heptosyltransferase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,795,5910GA21.8% 401.1 / 98.8 197S221S (TCG→TCAwaaFADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base G (86/68);  new base A (20/23);  total (106/91)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 202 total aligned reads displayed.

GCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGATCAA  >  NC_000913/3795493‑3795690
                                                                                                  |                                                                                                   
gCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCgg                                                                                                    >  2:3893544/1‑100 (MQ=255)
 cGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGc                                                                                                   >  1:2720591/1‑100 (MQ=255)
 cGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGc                                                                                                   >  2:1720095/1‑100 (MQ=255)
  ggAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCg                                                                                                  <  2:2072265/100‑1 (MQ=255)
     gTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGaaa                                                                                               <  1:3841397/100‑1 (MQ=255)
      tttGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAg                                                                                              <  1:4073552/100‑1 (MQ=255)
       ttGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGa                                                                                             >  1:1956670/1‑100 (MQ=255)
        tGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGAt                                                                                            >  1:2415803/1‑100 (MQ=255)
         ggTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATc                                                                                           <  1:2698676/100‑1 (MQ=255)
         ggTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATc                                                                                           <  1:2478766/100‑1 (MQ=255)
          ggccGGCAAAACGTTGGGCACACGACTAAGATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGAt                                                                                            <  2:899152/96‑1 (MQ=255)
          gTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCa                                                                                          >  1:1480582/1‑100 (MQ=255)
          gTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCa                                                                                          >  1:312584/1‑100 (MQ=255)
          gTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCa                                                                                          <  2:466439/100‑1 (MQ=255)
          gTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCa                                                                                          <  2:3692932/100‑1 (MQ=255)
            gcgacAAAACGCGGGCCACACTACCACCATGCGGAGCTGGCAGAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGGTTGGCTCGGCGAACGATCATg                                                                                        >  2:3079891/5‑100 (MQ=255)
            ccGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATg                                                                                        <  2:2276117/100‑1 (MQ=255)
             cGGCAAAACGCTGGCGACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGATCGGCGAAAGATCATg                                                                                        >  2:4144972/1‑99 (MQ=255)
             cGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGTAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATg                                                                                        <  2:3848450/99‑1 (MQ=255)
             cGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGa                                                                                       <  1:2841409/100‑1 (MQ=255)
              ggCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCACAGCAGCCGATTGCTGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGaa                                                                                      >  2:3735136/1‑100 (MQ=255)
              ggCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGaa                                                                                      >  1:1890844/1‑100 (MQ=255)
               gCAAAGCGCTGGACTAATTACGACGATCCGGAGCGTGGAAAGGAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATg                                                                                        <  2:3962289/97‑1 (MQ=255)
               gCAAAACGCTGGTCACATTAACAGTATGCGGAGCTGGCAAAGCAGCGGACTGATGAAGGTTATCAGGTGTTTTTGTTTGGCTAAGCGAAAGATAATGaa                                                                                      <  2:728438/99‑1 (MQ=255)
               gCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGATGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGCAg                                                                                     >  1:2404211/1‑100 (MQ=255)
                cAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGaa                                                                                      >  1:1350897/1‑98 (MQ=255)
                 aaaaCGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGc                                                                                    >  1:3951297/1‑99 (MQ=255)
                     cGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCa                                                                               <  2:3299497/100‑1 (MQ=255)
                       cTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAAt                                                                             >  1:425037/1‑100 (MQ=255)
                         ggCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATga                                                                           <  2:3401931/100‑1 (MQ=255)
                         ggCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATga                                                                           >  1:3045891/1‑100 (MQ=255)
                         ggCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATga                                                                           >  2:243776/1‑100 (MQ=255)
                         ggCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATga                                                                           >  2:90412/1‑100 (MQ=255)
                         ggCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATga                                                                           <  2:3856853/100‑1 (MQ=255)
                         ggCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATga                                                                           <  2:539046/100‑1 (MQ=255)
                          gCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATgag                                                                          >  1:125679/1‑100 (MQ=255)
                          gCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATgag                                                                          <  2:3923740/100‑1 (MQ=255)
                           ccACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATgaga                                                                         >  1:2249925/1‑100 (MQ=255)
                            cgcacTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGAt                                                                        <  2:2526891/98‑1 (MQ=255)
                            cacaCTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGAt                                                                        >  2:3568526/1‑100 (MQ=255)
                            cacaCTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGAt                                                                        >  1:51165/1‑100 (MQ=255)
                            cacaCTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGAt                                                                        >  2:3204651/1‑100 (MQ=255)
                             acacTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGAtt                                                                       <  1:1452225/100‑1 (MQ=255)
                                cTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGcttctt                                                                    >  2:736340/1‑100 (MQ=255)
                                cTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCtt                                                                    <  1:976928/100‑1 (MQ=255)
                                 tACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTg                                                                   >  1:854171/1‑100 (MQ=255)
                                  aCCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGAGTCTTg                                                                   <  2:1558673/99‑1 (MQ=255)
                                    cACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCg                                                                >  1:2987167/1‑100 (MQ=255)
                                    cACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCg                                                                <  1:1856054/100‑1 (MQ=255)
                                    cACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCg                                                                <  1:2052005/100‑1 (MQ=255)
                                     ncTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGc                                                               >  2:495802/2‑100 (MQ=255)
                                        aTGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCttt                                                            >  1:939601/1‑100 (MQ=255)
                                         tGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTg                                                           <  1:2523022/100‑1 (MQ=255)
                                         tGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTg                                                           <  1:2904206/100‑1 (MQ=255)
                                           cGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGaa                                                         <  2:86611/100‑1 (MQ=255)
                                           cGGAGCTGGCAAAGCAGCTGATTGATCAAGGTTATCAGGGGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTg                                                           <  2:1300343/98‑1 (MQ=255)
                                            ggAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAAt                                                        >  1:535380/1‑100 (MQ=255)
                                                  ggCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGag                                                  >  2:1379290/1‑100 (MQ=255)
                                                   gCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGagc                                                 >  1:2370989/1‑100 (MQ=255)
                                                    cAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACTGag                                                  >  2:2127841/1‑98 (MQ=255)
                                                    cAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTCGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTCAATACCGacc                                                 >  2:3953283/1‑97 (MQ=255)
                                                      aaGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGagcagc                                              <  2:1177124/100‑1 (MQ=255)
                                                       agcagcTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAgcagca                                             >  2:2613543/1‑100 (MQ=255)
                                                        gcagcTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGcagcag                                            >  1:550525/1‑100 (MQ=255)
                                                        gcagcTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGcagcag                                            >  1:324844/1‑100 (MQ=255)
                                                         cagcTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAgg                                           >  2:1516529/1‑100 (MQ=255)
                                                               attgatGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGgt                                     >  2:429479/1‑100 (MQ=255)
                                                                  gatgaAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTc                                  <  1:1683697/100‑1 (MQ=255)
                                                                   atgaAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCg                                 >  2:3538766/1‑100 (MQ=255)
                                                                    tgaAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCgg                                >  1:2778436/1‑100 (MQ=255)
                                                                    tgaAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCgg                                >  1:3457467/1‑100 (MQ=255)
                                                                    tgaAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCgg                                <  1:1585817/100‑1 (MQ=255)
                                                                       aGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAAc                             <  1:3560702/100‑1 (MQ=255)
                                                                          ttATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTg                          <  2:424282/100‑1 (MQ=255)
                                                                          ttATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTg                          >  1:1846609/1‑100 (MQ=255)
                                                                            aTCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGc                        >  2:1850404/1‑100 (MQ=255)
                                                                             tCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCg                       >  1:3659278/1‑100 (MQ=255)
                                                                             tCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCg                       >  2:2561947/1‑100 (MQ=255)
                                                                             tCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCg                       >  1:1153732/1‑100 (MQ=255)
                                                                               cggtggtTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAAGGAGATTCTTGCCGCTATGAATACCGAGCAGCAGGCATGGTGCCGGAACCTGGCggg                     >  2:2209677/2‑100 (MQ=255)
                                                                                 gtggtTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCggggg                   <  2:3175379/100‑1 (MQ=255)
                                                                                 gtggtTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCggggg                   <  2:1786454/100‑1 (MQ=255)
                                                                                 gtggtTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTATGCATACCGGGCAGCAGGCATGGTGTCGGAACCTGGCggggg                   >  2:2925648/1‑100 (MQ=255)
                                                                                  tggtTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGa                  <  2:1912305/100‑1 (MQ=255)
                                                                                   ggtTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGaa                 <  1:3925412/100‑1 (MQ=255)
                                                                                   ggtTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGaa                 >  1:2662283/1‑100 (MQ=255)
                                                                                     ttCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAac               <  1:2117695/100‑1 (MQ=255)
                                                                                     ttCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAac               <  1:1418474/100‑1 (MQ=255)
                                                                                        tGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGCCTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAaca              >  2:3260771/1‑98 (MQ=255)
                                                                                        tGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAAcaca            >  2:2886398/1‑100 (MQ=255)
                                                                                         gTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAg           >  1:1291794/1‑100 (MQ=255)
                                                                                         gTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAg           >  2:1800722/1‑100 (MQ=255)
                                                                                            tGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCtt        <  1:397060/100‑1 (MQ=255)
                                                                                            tGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCtt        <  2:1252675/100‑1 (MQ=255)
                                                                                             ggCTCGGCGAAAGATCATGTAGCGGGCAATGAGCTTCTTGCCGCTTTGAATACCGAGCAGCCGGCATGGTGTCGGAACTTGGCGGGGGAAACACAGCTTg       >  2:2783175/1‑100 (MQ=255)
                                                                                             ggCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCACTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTg       <  2:68125/100‑1 (MQ=255)
                                                                                              gCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGa      <  1:1105845/100‑1 (MQ=255)
                                                                                               cTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGAt     >  2:4017487/1‑100 (MQ=255)
                                                                                               cTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGAt     >  2:3010351/1‑100 (MQ=255)
                                                                                               cTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGAt     >  1:1001503/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGATCAA  >  NC_000913/3795493‑3795690

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