Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 197,248 G→A 6.2% A235T (GCA→ACA)  rseP → inner membrane zinc RIP metalloprotease; RpoE activator, by degrading RseA; cleaved signal peptide endoprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913197,2480GA 6.2% 218.9 / 4.1 64A235T (GCA→ACA) rsePinner membrane zinc RIP metalloprotease; RpoE activator, by degrading RseA; cleaved signal peptide endoprotease
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base A (2/2);  total (37/27)
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.49e-01

AACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGA  >  NC_000913/197111‑197386
                                                                                                                                         |                                                                                                                                          
aaCGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGgca                                                                                                                                          <  2:1741296/140‑1 (MQ=255)
  cggcggGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGgcaag                                                                                                                                        <  1:900119/140‑1 (MQ=255)
       ggATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAgg                                                                                                                                   >  1:34134/1‑140 (MQ=255)
       ggATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAgg                                                                                                                                   >  1:1539423/1‑140 (MQ=255)
       ggATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAgg                                                                                                                                   >  1:1539121/1‑140 (MQ=255)
            tAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTAACTTCGCTGGGGATTCGTCCTCGTGGGCCTCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGt                                                                                                                              >  2:1019046/1‑140 (MQ=255)
               aGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTg                                                                                                                           >  1:295986/1‑140 (MQ=255)
                  tCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCaa                                                                                                                        >  2:947440/1‑140 (MQ=255)
                    gATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCTGCAAGCAAGGCAGGTTTGCAAGc                                                                                                                      >  2:273921/1‑140 (MQ=255)
                    gATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCTGCAAGCAAGGCAGGTTTGCAAGc                                                                                                                      >  2:1400058/1‑140 (MQ=255)
                    gATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGc                                                                                                                      >  2:324268/1‑140 (MQ=255)
                    gATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGc                                                                                                                      >  2:323861/1‑140 (MQ=255)
                     aTTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCa                                                                                                                     <  2:631345/140‑1 (MQ=255)
                     aTTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCGGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCa                                                                                                                     <  2:512231/140‑1 (MQ=255)
                              aCTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAgg                                                                                                            >  1:1590210/1‑140 (MQ=255)
                              aCTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAgg                                                                                                            >  1:1590791/1‑140 (MQ=255)
                                tGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGAt                                                                                                          >  1:24898/1‑140 (MQ=255)
                                    cGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACGTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGtt                                                                                                      <  2:1407641/140‑1 (MQ=255)
                                    cGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGtt                                                                                                      <  2:151311/140‑1 (MQ=255)
                                    cGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGtt                                                                                                      >  1:1477696/1‑140 (MQ=255)
                                    cGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAAATGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGtt                                                                                                      <  2:1431827/140‑1 (MQ=255)
                                      tttGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTaa                                                                                                    >  1:332546/1‑140 (MQ=255)
                                       ttGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTaaa                                                                                                   >  2:778551/1‑140 (MQ=255)
                                            ccTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGa                                                                                              >  1:1372960/1‑140 (MQ=255)
                                                agacagaagaCCCGGTATCTTCGCTGGGGATTCGTCCTCGTGCGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCACACTCGGCGGCAAGCAAGCCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGt                                                                                          <  2:1404504/136‑1 (MQ=255)
                                                aTAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGt                                                                                          <  2:397920/140‑1 (MQ=255)
                                                aTAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGt                                                                                          <  2:1223430/140‑1 (MQ=255)
                                                aTAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGt                                                                                          >  2:87208/1‑140 (MQ=255)
                                                  acagaagaTCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAATTTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCTAGCTAGGCTGTTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTct                                                                                        <  2:1816504/138‑2 (MQ=255)
                                                  aaaGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                                                                        <  2:12810/140‑1 (MQ=255)
                                                  aaaGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                                                                        <  1:1860714/140‑1 (MQ=255)
                                                      aagaTCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGccc                                                                                    >  2:76178/1‑140 (MQ=255)
                                                        gaTCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTATAGTCGATGGTCAGCCCtt                                                                                  <  2:322881/140‑1 (MQ=255)
                                                        caTCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCtt                                                                                  <  1:1169383/139‑1 (MQ=255)
                                                         aTCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTa                                                                                 >  2:833204/1‑140 (MQ=255)
                                                         aTCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTa                                                                                 >  2:788201/1‑140 (MQ=255)
                                                             ggTCTCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGc                                                                             <  2:1811535/140‑1 (MQ=255)
                                                             ggTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAATCTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGc                                                                             <  2:1814077/140‑1 (MQ=255)
                                                             ggTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGc                                                                             <  1:893293/140‑1 (MQ=255)
                                                             ggTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGc                                                                             <  1:1671540/140‑1 (MQ=255)
                                                             ggTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGc                                                                             <  1:218814/140‑1 (MQ=255)
                                                             ggTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGc                                                                             <  1:769932/140‑1 (MQ=255)
                                                              gTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCa                                                                            <  1:18362/140‑1 (MQ=255)
                                                                aTCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGt                                                                          >  2:1798238/1‑140 (MQ=255)
                                                                                 gTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAATGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATg                                                         >  2:190057/1‑140 (MQ=255)
                                                                                 gTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATg                                                         >  2:190032/1‑140 (MQ=255)
                                                                                 gTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATg                                                         >  2:904719/1‑140 (MQ=255)
                                                                                  tCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGc                                                        >  1:1768048/1‑140 (MQ=255)
                                                                                      cGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGt                                                    >  1:787654/1‑140 (MQ=255)
                                                                                      cGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGt                                                    >  2:75617/1‑140 (MQ=255)
                                                                                            ccGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGa                                              <  2:1679200/140‑1 (MQ=255)
                                                                                                      caccagcaCTGGACAATGTGCACCCAACCTCGGCGGCAAGCAAGGCAGGATTGCAAGCAGCCGACAGGATCGATAAAGGCGATGGTCAGCCCTTAACGCAGTGGGTGACCATTGTGATGCTAGTCCGGGATAACCCGGGt                                    <  2:1287785/133‑1 (MQ=255)
                                                                                                       aCCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTa                                   <  2:1753016/140‑1 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGTTGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAg                    >  1:1029406/1‑140 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAg                    >  2:903530/1‑140 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAg                    >  1:1031378/1‑140 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAg                    >  1:1664288/1‑140 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAg                    >  2:1459095/1‑140 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGACAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                                                                        <  1:1370096/72‑1 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGACAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                                                                        >  2:806676/1‑72 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGACAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                                                                        >  2:1370096/1‑72 (MQ=255)
                                                                                                                      tgtgCAGCCAAACTCGGCGACAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                                                                        <  1:806676/72‑1 (MQ=255)
                                                                                                                        tgCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGaa                  >  2:898456/1‑140 (MQ=255)
                                                                                                                              cAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGaa            <  2:1039668/140‑1 (MQ=255)
                                                                                                                                aaCTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAg          <  1:1475433/140‑1 (MQ=255)
                                                                                                                                 aCTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAgg         >  1:813222/1‑140 (MQ=255)
                                                                                                                                        ggcAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGTCATTTGTGATGCTTGTCCGGGCTAACCCGGGTAAATCCTTCGCGTTAGAAATCGAAAGGCAGGGGa  <  2:144274/140‑1 (MQ=255)
                                                                                                                                        ggcAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATATCCCGGGTTAATCCTTAGCGTAAGAAATCGAAAGGCAGGGGa  <  2:144283/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
AACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGA  >  NC_000913/197111‑197386

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