Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 98,451 A→T 6.4% S17C (AGT→TGT)  ftsW → putative lipid II flippase; integral membrane protein; FtsZ ring stabilizer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091398,4510AT 6.4% 184.0 / 3.1 72S17C (AGT→TGT) ftsWputative lipid II flippase; integral membrane protein; FtsZ ring stabilizer
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base T (3/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

ATATGGTTCTGCTCTCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATGGTGACCTCGG  >  NC_000913/98303‑98598
                                                                                                                                                    |                                                                                                                                                   
atatGGTTCTGCTCTCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCa                                                                                                                                                     <  1:215411/149‑1 (MQ=255)
  atGGTTCTGCTCTCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGt                                                                                                                                                   >  1:623342/1‑149 (MQ=255)
      ttCTGCTCTCCCCAGCCTGTGCCAGCCTTGCTCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTCTCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATcc                                                                                                                                               >  2:538687/1‑149 (MQ=255)
       tCTGCTCTCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCACTCGCCAGAAAATGCCGCGCCTGCCAGGATTCAGTATGCt                                                                                                                                              >  2:1089444/1‑149 (MQ=255)
        cTGCTCTCCCCCGCCTCTGCCAGCCTTGATCCGTTCAACAACTTTGAACAACGAGGCAAAGAGTTTGCCCGGCTGGCGAAGGAGCTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTg                                                                                                                                             <  2:127907/149‑1 (MQ=255)
          gCTCTCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGt                                                                                                                                           >  2:232581/1‑149 (MQ=255)
          gCTCTCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGt                                                                                                                                           <  1:400210/149‑1 (MQ=255)
          gCTCTCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATACTGAATCCTggc                                                                                                                                           <  1:368369/149‑2 (MQ=255)
           ctctCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTc                                                                                                                                          >  2:591473/1‑149 (MQ=255)
            tctcCCCAGCGTGTGCGAGCCTTGATCAGTTCAAGAACATTGCGCAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGGTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGGATCGTGGTCt                                                                                                                                         >  2:1130939/1‑149 (MQ=255)
            tctcCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCt                                                                                                                                         >  1:944017/1‑149 (MQ=255)
              tcCCCAGCCTGTGCCACCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAAAGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCAGGTCTgg                                                                                                                                       >  2:124609/1‑149 (MQ=255)
              accccAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGGCTGGCGAAGGAGTTAGGTTGATTCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTAGCCTGGTCTgg                                                                                                                                       <  1:786464/148‑1 (MQ=255)
                cccAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGGTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGAt                                                                                                                                     <  1:472126/149‑1 (MQ=255)
                cccAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTATGAACAACGAGGCAATGAGTATGCCCGTCAGGCGAAGGAGTTAGGTTGATGCGTTTAACTCTCCCTCGCCTGAAAATGCCGCGCATGCCAGCATTCAGTATCCGGGACAGGAt                                                                                                                                     >  1:1148283/1‑149 (MQ=255)
                   aGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCATCGCCTGAAAATGCCGCGACTGCCAGGATTCAGTATCCTGGTCTGGAtctc                                                                                                                                  >  1:602454/1‑149 (MQ=255)
                     ccTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGTAAATGAGTTTTACCGTCTGGCTAAGGAGTTAGGTAGATGCGATTATCTCTCCCTCGCCTGAAAATTCCGAGCCTGCCATTATTCAGAATCCTAGTCTGGATCTCCa                                                                                                                                >  2:661743/1‑149 (MQ=255)
                        gtgCCAGCCTTGATCAGTTCAAGAACTTTGAATAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACgg                                                                                                                             <  2:312457/149‑1 (MQ=255)
                             aGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTATGGATCTCCACGGCGCTa                                                                                                                        >  1:789713/1‑149 (MQ=255)
                                  tgctccgTTCAACCACCTTGAACAACCAGGCAATGCCTTTGCCCGGCGGGCGAAGGAGGTACGGTGCTGCGGTTATCCCCCCCCCGCCCGAAAATGCCCCCCCGGCCACGCTTCCGGATCCCTGGCGGGCTCTCCACGGCGCTAAAggg                                                                                                                   <  2:338843/143‑1 (MQ=255)
                                   gATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGc                                                                                                                  >  2:549674/1‑149 (MQ=255)
                                        gTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTTGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGt                                                                                                             <  1:1127845/149‑1 (MQ=255)
                                           cAAGAACTTTGAACAACGACGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGAt                                                                                                          <  2:335981/149‑1 (MQ=255)
                                            aaGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATg                                                                                                         <  2:625394/149‑1 (MQ=255)
                                             aGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTCGGTAGATGCGTTTATCTCTCCCTCGCCTGAAATTGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGTATCTCCACGGCGCTAAAGGGCTGGGTGATgg                                                                                                        <  2:784497/149‑1 (MQ=255)
                                              gAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGCAGGAGTTAGGTTGATGCGGTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGAATCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATggg                                                                                                       <  2:240510/149‑1 (MQ=255)
                                                     gAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGGTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCgg                                                                                                <  1:829734/149‑1 (MQ=255)
                                                      aacaacGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGCTGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGCTGGGCTCGCggg                                                                                               >  1:402802/1‑149 (MQ=255)
                                                      aacaacGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGTTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCggg                                                                                               >  1:562279/1‑149 (MQ=255)
                                                        caacGAGGCAATGAGTTTGCCCGTCTCGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCCCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGaa                                                                                             >  1:617121/1‑149 (MQ=255)
                                                            gAGGCAATGCGTATGCCCGTCTGGCGTCGGAGTTTGCTTGATGGGATAATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTTTCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAg                                                                                         <  2:476519/149‑1 (MQ=255)
                                                             aGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGa                                                                                        >  2:236644/1‑149 (MQ=255)
                                                                           cccGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGGGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGa                                                                          <  1:841335/149‑1 (MQ=255)
                                                                           cccGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCAGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTgc                                                                          >  2:688344/1‑148 (MQ=255)
                                                                              gTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGCCAGCCTGATCa                                                                       >  1:201902/1‑149 (MQ=255)
                                                                                cTGGCGAAGGAGTTAGGTTGAGGCGTTTATCTCTCCCTCGCCTGAAACTGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTACACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATg                                                                     <  2:239879/149‑1 (MQ=255)
                                                                                 tGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAATATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGt                                                                    <  1:232581/149‑1 (MQ=255)
                                                                                  ggCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATCCCGCGCCTGCCATGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTa                                                                   <  1:124609/149‑1 (MQ=255)
                                                                                    cGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACg                                                                 >  2:1148840/1‑149 (MQ=255)
                                                                                     gAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGa                                                                >  1:850186/1‑149 (MQ=255)
                                                                                     gAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGa                                                                >  1:846132/1‑149 (MQ=255)
                                                                                           ggtAGGTTGATGCTGTTATCACTCCCTCGCCTGAAAATGCCCCGCCCGCGAGGATTCAGTATCCTGGTCTCGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCAc                                                          <  2:268564/147‑1 (MQ=255)
                                                                                               ggatGATGCGTTAATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTAACCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTa                                                      <  1:347991/146‑1 (MQ=255)
                                                                                               ggTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTa                                                      <  1:33518/149‑1 (MQ=255)
                                                                                                     tGCGTTTATCTCTCCCTCTCCTGAAAATGCCGCGGCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGATAATGATACCGACAGCCTGATCATGTACGAGCGCACCTTACTGTgg                                                <  1:661743/149‑1 (MQ=255)
                                                                                                        gTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGCTACCGACAGCCTGATCATGGACGATCGCACCTTACTGTGGCTg                                             >  2:314015/1‑149 (MQ=255)
                                                                                                              ctctcCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCGGGGTGATGGGCGCGCGGGAAAAAGATACCGACAGCCTGATCCTGTACGAGCGCACCTTACTGTGGCTGACCTTc                                       >  1:848518/1‑149 (MQ=255)
                                                                                                               tctcCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCg                                      >  2:243658/1‑149 (MQ=255)
                                                                                                                ctcCCTCGCCTGTAAATGCGGCGCCTGCCAGGTTTCAGGAGCCTGTTCTTGATCTCCACGGCGCTAAAGGGCGGGGTGATGGGCTAGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGAACCTTACTGTGGCTGACCTTCgg                                     <  2:58157/149‑1 (MQ=255)
                                                                                                                ctcCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCCCGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGCTACCGACAGCCTGATCATGTACGATCGCACCTGACTGTGGCTGACCTTCgg                                     >  1:281925/1‑149 (MQ=255)
                                                                                                                ctcCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCgg                                     <  1:549674/149‑1 (MQ=255)
                                                                                                                ctcCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATACTGGTCTGGATCTCCACGGCGGTTAAGGGCTGGGTGATGGGCTCGCGGGAAAATGATACCGACAGCCTGATCTTGTACGATCGCACCTTACTGTGGCTGCCCTTCgg                                     >  1:1087176/1‑149 (MQ=255)
                                                                                                                ctcCCTCGCCTGAAAATACCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCgg                                     >  1:456347/1‑149 (MQ=255)
                                                                                                                cacccTCGCCTGAAACTGCCGCGCCTGCCAGGATTCAGTCTCCTGGTCGGGATCTCCACCGCGCTAAAGGTCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCgg                                     <  2:614902/147‑1 (MQ=255)
                                                                                                                     tCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTcg                                >  2:594280/1‑149 (MQ=255)
                                                                                                                        ccTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGCCCTTCGGCCTCgcgg                             >  2:351668/1‑149 (MQ=255)
                                                                                                                         cTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGATAAAGATACCGACAGCCTGATCATGGACGATCGCACCTTAGTGTGGCTGACCTTCGGCCTCgcggc                            >  1:822368/1‑149 (MQ=255)
                                                                                                                            aaaaTGCCGCGCCTGCCAGGATTCAGTATCGTGGTCTGGATCTCCACGGCGCTATAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTCCTGTGGCTGCCCTTCGGCAGCGCGGCGag                         >  1:365359/1‑148 (MQ=255)
                                                                                                                            aaaaTGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGAt                         <  1:1148840/149‑1 (MQ=255)
                                                                                                                                tGCCGCGCCTGCCAGGATTCTGTAGCCTGGGCGGGAGCTCCACGGCGCTAAAGGGCGGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGc                     <  1:1051134/149‑1 (MQ=255)
                                                                                                                                   cgcgCCTGCCAGGATTCAGTATCCTGGTCTGGAGCTCCACGGCGCTAAACGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCttt                  <  1:591473/149‑1 (MQ=255)
                                                                                                                                      ccagaccAGGATACTGAATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCCCACCATACTCTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTAta               >  2:995316/6‑148 (MQ=255)
                                                                                                                                       cagaccAGGATACTGAATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCa              >  1:94718/5‑149 (MQ=255)
                                                                                                                                         tGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATg            <  2:955350/149‑1 (MQ=255)
                                                                                                                                          accAGGATACTGAATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATgg           >  1:978703/2‑149 (MQ=255)
                                                                                                                                           ccAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATGGt          <  1:1043304/149‑1 (MQ=255)
                                                                                                                                           ccAGGATGCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATGGt          <  2:623342/149‑1 (MQ=255)
                                                                                                                                            cAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATGGTg         <  1:594280/149‑1 (MQ=255)
                                                                                                                                             aGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTa                 <  2:1104276/140‑1 (MQ=255)
                                                                                                                                             aGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTa                 >  1:1104276/1‑140 (MQ=255)
                                                                                                                                              ggATTCAGTATCCTGGTCTTGATCTCCACGGCGCTAAACGGCTGGGTGATGGTCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATGGTGAc       <  2:861552/149‑1 (MQ=255)
                                                                                                                                               gATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATGGTGAcc      >  1:578365/1‑149 (MQ=255)
                                                                                                                                                   cAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATGGTGACCTCgg  <  2:944017/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATATGGTTCTGCTCTCCCCAGCCTGTGCCAGCCTTGATCAGTTCAAGAACTTTGAACAACGAGGCAATGAGTTTGCCCGTCTGGCGAAGGAGTTAGGTTGATGCGTTTATCTCTCCCTCGCCTGAAAATGCCGCGCCTGCCAGGATTCAGTATCCTGGTCTGGATCTCCACGGCGCTAAAGGGCTGGGTGATGGGCTCGCGGGAAAAAGATACCGACAGCCTGATCATGTACGATCGCACCTTACTGTGGCTGACCTTCGGCCTCGCGGCGATTGGCTTTATCATGGTGACCTCGG  >  NC_000913/98303‑98598

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