Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,143,312 T→G T261T (ACT→ACG)  frwC → putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,143,3120TG98.9% 632.0 / ‑2.8 176T261T (ACT→ACG) frwCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (84/90);  total (84/92)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 176 total aligned reads displayed.

GATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCG  >  NC_000913/4143216‑4143409
                                                                                                 |                                                                                                 
gATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGaa                                                                                                >  2:3852892/1‑100 (MQ=255)
 aTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAg                                                                                               >  1:3836438/1‑100 (MQ=255)
 aTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAg                                                                                               <  2:3233948/100‑1 (MQ=255)
  tGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAgg                                                                                              >  2:708694/1‑100 (MQ=255)
     ggCTGGCGACGCTGATTGGTTGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATTGGCTGTGTT‑GGCGGTACGGAAGGGGc                                                                                           >  2:141059/1‑100 (MQ=255)
     ggCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGc                                                                                           <  2:5479576/100‑1 (MQ=255)
     ggCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGc                                                                                           <  2:3410671/100‑1 (MQ=255)
      gCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCg                                                                                          <  1:2586692/100‑1 (MQ=255)
      gCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCg                                                                                          >  2:1187346/1‑100 (MQ=255)
        tGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGAt                                                                                        <  1:1326922/100‑1 (MQ=255)
        tGGCGACGCGGATTGGTCGTAAAAATTTCTCCGAAGAAGAGCGCGAAACCGGCAAAGCGGCAATGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGAt                                                                                        <  2:2167990/100‑1 (MQ=255)
         ggCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGAtt                                                                                       >  2:5385898/1‑100 (MQ=255)
         ggCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGAtt                                                                                       >  1:1055596/1‑100 (MQ=255)
           cGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTcc                                                                                     >  1:4478979/1‑100 (MQ=255)
            gACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCt                                                                                    >  1:2166185/1‑100 (MQ=255)
             aCGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCtt                                                                                   <  2:4144169/100‑1 (MQ=255)
             aCGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCtt                                                                                   <  2:3103976/100‑1 (MQ=255)
              cGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCttt                                                                                  <  1:3657598/100‑1 (MQ=255)
                  gATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTcgcc                                                                              <  1:79757/100‑1 (MQ=255)
                     tGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCt                                                                           >  2:2784604/1‑100 (MQ=255)
                     tGGTCGTAAAAATTACTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGAGGCGATTCCTTTCGCCGCt                                                                           <  2:5033361/100‑1 (MQ=255)
                      ggTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTg                                                                          >  2:874212/1‑100 (MQ=255)
                      ggTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTg                                                                          <  1:2971991/100‑1 (MQ=255)
                       gTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTAATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGc                                                                         >  2:3739351/1‑100 (MQ=255)
                        tCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGcc                                                                        <  2:2050832/100‑1 (MQ=255)
                        tCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGcc                                                                        <  1:2027018/100‑1 (MQ=255)
                        tCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGcc                                                                        <  1:1441644/100‑1 (MQ=255)
                         cGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCg                                                                       <  2:3993278/100‑1 (MQ=255)
                          gTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGa                                                                      >  2:4941030/1‑100 (MQ=255)
                           tAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGAt                                                                     >  2:468859/1‑100 (MQ=255)
                             aaaaTTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATcc                                                                   >  2:1236255/1‑100 (MQ=255)
                              aaaTTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGTCGATTCCTTTCGCCGCTGCCGATccn                                                                  <  2:4730452/100‑2 (MQ=255)
                                  ttCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGc                                                              >  1:5483235/1‑100 (MQ=255)
                                    ctcCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCgt                                                            >  2:426330/1‑100 (MQ=255)
                                    ctcCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCgt                                                            >  1:1921533/1‑100 (MQ=255)
                                     tcCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCgtg                                                           <  1:2749052/100‑1 (MQ=255)
                                      ccGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCgtgt                                                          >  1:3397723/1‑100 (MQ=255)
                                       cGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTCCGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTc                                                         <  1:2614278/100‑1 (MQ=255)
                                        gCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCa                                                        <  2:5211439/100‑1 (MQ=255)
                                          agaagaGCGCGAAACCGGCAAAGCGGCACTGGTGATGGACTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTGCGCCGCTGCCGATCCGCTGCGTGTCAtt                                                      <  2:5299179/100‑1 (MQ=255)
                                            aagaGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttcc                                                    <  1:717539/100‑1 (MQ=255)
                                             agagCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttcct                                                   <  2:1116943/100‑1 (MQ=255)
                                                gcgcGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttccttcc                                                <  1:3799561/100‑1 (MQ=255)
                                                gcgcGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttccttcc                                                >  1:116487/1‑100 (MQ=255)
                                                   cgAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCcatc                                             <  1:5308560/100‑1 (MQ=255)
                                                    gAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCcatca                                            <  1:344455/100‑1 (MQ=255)
                                                     aaaCCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTGATTCCTTCcatcat                                           >  1:3012174/1‑100 (MQ=255)
                                                      aaCCGGCAAAGCTGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATg                                          <  2:1996725/100‑1 (MQ=255)
                                                      aaCCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATg                                          >  2:1711583/1‑100 (MQ=255)
                                                        ccGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGt                                        >  1:1117079/1‑100 (MQ=255)
                                                         cGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTc                                       <  2:4592955/100‑1 (MQ=255)
                                                          ggCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCg                                      >  2:2665483/1‑100 (MQ=255)
                                                          ggCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCg                                      >  2:2288012/1‑100 (MQ=255)
                                                          ggCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCg                                      >  2:2288201/1‑100 (MQ=255)
                                                           gCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCgg                                     <  1:5027131/100‑1 (MQ=255)
                                                           gCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCgg                                     <  1:3962640/100‑1 (MQ=255)
                                                            cAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGt                                    >  2:70445/1‑100 (MQ=255)
                                                             aaaGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGtt                                   <  1:2311735/100‑1 (MQ=255)
                                                             aaaGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGtt                                   <  1:1203421/100‑1 (MQ=255)
                                                              aaGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTc                                  <  1:1793006/100‑1 (MQ=255)
                                                                gCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAg                                >  1:2156729/1‑100 (MQ=255)
                                                                gCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAg                                <  2:2479229/100‑1 (MQ=255)
                                                                  ggCACTGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGtt                              <  2:3328151/100‑1 (MQ=255)
                                                                       tGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTgtgg                         <  2:256361/100‑1 (MQ=255)
                                                                       tGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTgtgg                         <  2:2134163/100‑1 (MQ=255)
                                                                       tGGTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTgtgg                         <  1:590567/100‑1 (MQ=255)
                                                                        ggTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTgtggt                        <  2:1933275/100‑1 (MQ=255)
                                                                        ggTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTgtggt                        <  1:1731625/100‑1 (MQ=255)
                                                                        ggTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTgtggt                        <  1:3795303/100‑1 (MQ=255)
                                                                         gTGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCCGTTTgtggtg                       >  1:4773462/1‑100 (MQ=255)
                                                                          tggtGGGCT‑TGTTGGGCGTTACGGAAGGGGCGGTTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTAGTGc                      <  2:4671806/97‑1 (MQ=255)
                                                                          tcttGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGc                      <  1:2946584/97‑1 (MQ=255)
                                                                          tGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGc                      >  1:3117768/1‑100 (MQ=255)
                                                                          tGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGc                      >  2:1351391/1‑100 (MQ=255)
                                                                          tGATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGc                      <  1:4192866/100‑1 (MQ=255)
                                                                           gATGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCa                     <  1:1054240/100‑1 (MQ=255)
                                                                            aTGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAg                    <  2:1245542/100‑1 (MQ=255)
                                                                            aTGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAg                    <  2:3299866/100‑1 (MQ=255)
                                                                            aTGGGCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAg                    <  1:4127937/100‑1 (MQ=255)
                                                                              gggCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTa                  >  2:5055362/1‑100 (MQ=255)
                                                                              gggCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTa                  <  2:4435431/100‑1 (MQ=255)
                                                                                gCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAAc                <  1:3437273/100‑1 (MQ=255)
                                                                                gCTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAAc                <  1:5352011/100‑1 (MQ=255)
                                                                                gCTGTGTT‑GGCATTACGGAAAGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAAc                <  1:4321745/100‑1 (MQ=255)
                                                                                 cTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACt               >  2:91791/1‑100 (MQ=255)
                                                                                 cTGTGTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACt               <  2:3405740/100‑1 (MQ=255)
                                                                                  tgtgTT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTg              >  2:82379/1‑100 (MQ=255)
                                                                                   gtgtT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTgc             >  2:1209865/1‑100 (MQ=255)
                                                                                   gtgtT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTgc             >  1:2426499/1‑100 (MQ=255)
                                                                                    tgtT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCGTGGTCGGTTCAGTTTGTGGTGCAGTAACTgcg            >  1:5106597/1‑100 (MQ=255)
                                                                                    tgtT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTgcg            >  1:4102567/1‑100 (MQ=255)
                                                                                     gtT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTACCTgcgg           >  2:1449115/1‑100 (MQ=255)
                                                                                     gtT‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTgcgg           >  1:1562868/1‑100 (MQ=255)
                                                                                     ggt‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCAGGGTCGGTTCAGTTTGTGGTGCAGTAACTgcgg           >  2:388991/3‑100 (MQ=255)
                                                                                      tt‑GGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGACGCGTTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTgcggc          <  2:5248146/100‑1 (MQ=255)
                                                                                          gCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCt       <  2:4334753/100‑1 (MQ=255)
                                                                                           cGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTg      <  2:1486762/100‑1 (MQ=255)
                                                                                            gTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTgg     >  2:5282239/1‑100 (MQ=255)
                                                                                              tACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTc   >  1:373882/1‑100 (MQ=255)
                                                                                               aCTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCg  <  1:3316815/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
GATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTT‑GGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCG  >  NC_000913/4143216‑4143409

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