Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,966,414 T→G 100% D674A (GAT→GCT)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,966,4140TG97.2% 412.7 / 0.8 122D674A (GAT→GCT) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (52/66);  total (54/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 123 total aligned reads displayed.

ATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGAC  >  NC_000913/2966318‑2966510
                                                                                                |                                                                                                
aTGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAgcgc                                                                                               >  1:1239872/1‑100 (MQ=255)
   cgccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGcgcccg                                                                                            >  2:133452/1‑100 (MQ=255)
   cgccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGcgcccg                                                                                            <  2:2047902/100‑1 (MQ=255)
   cgccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGcgcccg                                                                                            >  1:500843/1‑100 (MQ=255)
     ccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGcgcccgcc                                                                                          <  1:1827007/100‑1 (MQ=255)
         agcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCAttt                                                                                      <  1:2553733/100‑1 (MQ=255)
          gcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTc                                                                                     <  2:125576/100‑1 (MQ=255)
          gcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTc                                                                                     <  1:630798/100‑1 (MQ=255)
           cagGTATTTTGCCCGCGCTACAGAACGCCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCCCACGGGAGcccccgcccctca                                                                                     >  2:2009261/1‑94 (MQ=255)
            agGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCAc                                                                                   <  1:1506009/100‑1 (MQ=255)
            agGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCAc                                                                                   >  2:2483440/1‑100 (MQ=255)
             ggTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCAcc                                                                                  <  1:2258408/100‑1 (MQ=255)
             ggTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCAcc                                                                                  <  1:1933493/100‑1 (MQ=255)
             ggTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCAcc                                                                                  <  1:2554327/100‑1 (MQ=255)
              gTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGTGAGCGCCCGCCATTTCACCg                                                                                 <  2:1980603/100‑1 (MQ=255)
               tATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGc                                                                                <  1:2550254/100‑1 (MQ=255)
                 ttttGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGcac                                                                              <  1:1351913/100‑1 (MQ=255)
                 ttttGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGcac                                                                              <  1:1613233/100‑1 (MQ=255)
                     gCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAAc                                                                          >  1:67331/1‑100 (MQ=255)
                     gCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAAc                                                                          >  1:1329382/1‑100 (MQ=255)
                      cccGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCACTGAGGATTGCCACGCACATAGGAGCGCCCGCCATATCACCGCACAAAc                                                                          >  2:245492/1‑99 (MQ=255)
                      cccGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACg                                                                         >  1:1277597/1‑100 (MQ=255)
                        cgcgCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGga                                                                       <  1:870851/100‑1 (MQ=255)
                         gcgcTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGgag                                                                      >  1:2065793/1‑100 (MQ=255)
                            cTACAGAGCGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                   <  1:133454/100‑1 (MQ=255)
                            cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                   <  1:113220/100‑1 (MQ=255)
                            cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                   <  2:838101/100‑1 (MQ=255)
                            cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                   <  1:2249281/100‑1 (MQ=255)
                              aCAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCACCCGCCATTTCACCGCACAAACGGAGATCGa                                                                 >  2:2361015/1‑100 (MQ=255)
                                 gAACGTCCGTTCATAGACAGATGGCGTTACCCAAGCCCAATGGGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGTTCgct                                                                >  2:73093/1‑96 (MQ=255)
                                  aaCGTACGTTCATAGACAGATGGCGACACCCAAGCCCAATGAGGATTGCCCCGCACATGGAAGCGCCCGCCATTTCACCGCACAAACGGAGATCGAtt                                                               <  2:2043502/98‑1 (MQ=255)
                                   aCGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCa                                                            <  2:480746/100‑1 (MQ=255)
                                    cGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCa                                                            <  2:1858137/99‑1 (MQ=255)
                                    cGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCAt                                                           <  2:2150339/100‑1 (MQ=255)
                                    cGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCAt                                                           <  2:1928389/100‑1 (MQ=255)
                                     gTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATg                                                          <  2:2047689/100‑1 (MQ=255)
                                     gTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATg                                                          <  2:1204774/100‑1 (MQ=255)
                                        cGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTAt                                                       >  1:1857499/1‑100 (MQ=255)
                                         gctCATAGACAGAGGGCGATACCCAACCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTa                                                        <  2:2298833/96‑1 (MQ=255)
                                          ttCATAGACAGATGGCGATACCCAAGCGCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTAttt                                                     >  2:556407/1‑100 (MQ=255)
                                            cATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTcc                                                   >  2:1886080/1‑100 (MQ=255)
                                            cATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTcc                                                   <  1:1980389/100‑1 (MQ=255)
                                              tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                                 >  1:1064043/1‑100 (MQ=255)
                                                gACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCtt                                               <  2:1098814/100‑1 (MQ=255)
                                                  cAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTcc                                             >  1:1392923/1‑100 (MQ=255)
                                                   aGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTccc                                            <  1:991369/100‑1 (MQ=255)
                                                      tGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCggg                                         <  1:2345643/100‑1 (MQ=255)
                                                      tGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCggg                                         >  2:1172556/1‑100 (MQ=255)
                                                      tGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCggg                                         <  2:1979920/100‑1 (MQ=255)
                                                         cGATACCCAAGCCCAATGAGGATTGCCACGCACATTGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGa                                      <  1:2491995/100‑1 (MQ=255)
                                                          gATACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGat                                     >  2:1054556/1‑100 (MQ=255)
                                                            tACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCGCAAACGGAGGTCGATACCCGGTATTTCCGCTTCCCGGGCg                                       >  2:1685684/1‑96 (MQ=255)
                                                            tACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatca                                   <  1:1300310/100‑1 (MQ=255)
                                                            tACCCAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatc                                    <  1:185646/99‑1 (MQ=255)
                                                              cccAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTCCCCGGGCGatcatc                                 >  2:1648640/1‑100 (MQ=255)
                                                              cccAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                                 >  2:1594617/1‑100 (MQ=255)
                                                              cccAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                                 >  2:1094761/1‑100 (MQ=255)
                                                              cccAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGAGCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                                 >  2:1089218/1‑100 (MQ=255)
                                                               ccAAGCCGATTGACGGTTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCCCATACGGAGATCGATTCCATTTATTTCCGCTGCCCGGGCGATCATCg                                <  2:916685/100‑1 (MQ=255)
                                                                cAAGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGc                               >  2:2506068/1‑100 (MQ=255)
                                                                 aaGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGcc                              <  1:732731/100‑1 (MQ=255)
                                                                  aGCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCa                             <  1:1830571/100‑1 (MQ=255)
                                                                   gCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCag                            <  1:939154/100‑1 (MQ=255)
                                                                   gCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCag                            <  1:1315857/100‑1 (MQ=255)
                                                                   gCCCAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCag                            <  2:526063/100‑1 (MQ=255)
                                                                    cccAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCaga                           >  2:156823/1‑100 (MQ=255)
                                                                    cccAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCaga                           <  2:870538/100‑1 (MQ=255)
                                                                     ccAATGAGGATTGCCACGCACATGGGAGCGCCGGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCagag                          <  2:28476/100‑1 (MQ=255)
                                                                      cAATGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGc                         <  1:2097237/100‑1 (MQ=255)
                                                                       aaTGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCt                        >  1:1137116/1‑100 (MQ=255)
                                                                        aTGAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACCAACGGAGATCGATTCCATGTATTTCCGCTGCCCGGGCGATGATCGcccgggaag                       <  2:1664556/100‑8 (MQ=255)
                                                                          gAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGt                     >  1:1191642/1‑100 (MQ=255)
                                                                          gAGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGt                     <  1:2344518/100‑1 (MQ=255)
                                                                           aGGATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAATCGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTa                    <  1:2175972/100‑1 (MQ=255)
                                                                            ggATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTaa                   >  1:441362/1‑100 (MQ=255)
                                                                            ggATTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTaa                   <  2:678374/100‑1 (MQ=255)
                                                                              aTTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACa                 >  1:1490403/1‑100 (MQ=255)
                                                                              aTTGCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACCGAGATCGATTCCATGTATTTCCGCTTCCCGCGCGATCATCGCCAGAGCTCg                      >  2:2253331/1‑95 (MQ=255)
                                                                                 gCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg              >  2:568916/1‑100 (MQ=255)
                                                                                 gCCACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg              >  1:1023892/1‑100 (MQ=255)
                                                                                  ccACGCACATGGGAGCTCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGc             <  2:191885/100‑1 (MQ=255)
                                                                                  ccACGCACATGGGAGCACCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg              >  1:220377/1‑99 (MQ=255)
                                                                                   cACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa            <  1:2483454/100‑1 (MQ=255)
                                                                                   cACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa            >  1:399219/1‑100 (MQ=255)
                                                                                   cACGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTCCa            >  2:746741/1‑100 (MQ=255)
                                                                                    aCGCACATGGGAGCGCCCGCCATTTCACCGCACACACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGc             <  2:851511/98‑1 (MQ=255)
                                                                                    aCGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa            <  2:2265333/99‑1 (MQ=255)
                                                                                    aCGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAg           >  1:1753123/1‑100 (MQ=255)
                                                                                    aCGCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAg           <  1:2281119/100‑1 (MQ=255)
                                                                                      gCACATGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGa         <  1:2521965/100‑1 (MQ=255)
                                                                                      gCACATGGGAGCACCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGa         <  1:2361030/100‑1 (MQ=255)
                                                                                       cacaTGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGAt        <  2:2008403/100‑1 (MQ=255)
                                                                                        acaTGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCTATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATg       <  1:959103/100‑1 (MQ=255)
                                                                                         caTGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTCACATTGCAGGATg       >  1:940652/1‑99 (MQ=255)
                                                                                         caTGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATg       >  2:2186187/1‑99 (MQ=255)
                                                                                         caTGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGa      >  2:2032443/1‑100 (MQ=255)
                                                                                         caTGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGa      <  1:529517/100‑1 (MQ=255)
                                                                                           tGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAg    >  2:801814/1‑100 (MQ=255)
                                                                                           tGGGAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAg    >  1:767150/1‑100 (MQ=255)
                                                                                             ggAGCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGAc  >  1:1053813/1‑100 (MQ=255)
                                                                                                |                                                                                                
ATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGAC  >  NC_000913/2966318‑2966510

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: