Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,359,353 G→C 5.3% A306G (GCC→GGC)  nusA ← transcription termination/antitermination factor L

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,359,3530GC 5.3% 316.9 / 6.1 95A306G (GCC→GGC) nusAtranscription termination/antitermination factor L
Reads supporting (aligned to +/- strand):  ref base G (43/47);  new base C (3/2);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 105 total aligned reads displayed.

ACCGGTTTCGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCTGGACGCGCGAACC  >  NC_002947/5359205‑5359498
                                                                                                                                                    |                                                                                                                                                 
aCCGGTTTCGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCcatgg                                                                                                                                                   >  1:368362/1‑149 (MQ=255)
aCCGGTTTCGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCcatgg                                                                                                                                                   >  3:56686/1‑149 (MQ=255)
aCCGGTTTCGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCcatgg                                                                                                                                                   >  7:149130/1‑149 (MQ=255)
    gTTTCGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCcatggcatg                                                                                                                                               >  4:265743/1‑149 (MQ=255)
     tttCGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCcatggcatgg                                                                                                                                              <  4:333476/149‑1 (MQ=255)
      ttCGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATggg                                                                                                                                             >  1:214570/1‑149 (MQ=255)
        cGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGc                                                                                                                                            <  7:316402/148‑1 (MQ=255)
           ccTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCt                                                                                                                                        >  2:222187/1‑149 (MQ=255)
            cTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCtt                                                                                                                                       <  6:155422/149‑1 (MQ=255)
             tGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGCCCATGCCATGGACATc                                                                                                                                         <  8:237859/146‑1 (MQ=255)
               ttGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTtc                                                                                                                                      <  4:207833/147‑1 (MQ=255)
                  cTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCaa                                                                                                                                 <  1:261951/149‑1 (MQ=255)
                   ttGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATAAAc                                                                                                                                >  7:137910/1‑149 (MQ=255)
                    tGGCCGGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAAc                                                                                                                                >  7:369605/1‑148 (MQ=255)
                      gCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACga                                                                                                                              <  1:222187/148‑1 (MQ=255)
                          ggATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATc                                                                                                                         >  7:152769/1‑149 (MQ=255)
                          ggATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATc                                                                                                                         >  2:140415/1‑149 (MQ=255)
                              gTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATcgccg                                                                                                                     >  1:257679/1‑149 (MQ=255)
                               tCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATTATcgccgc                                                                                                                    <  2:368362/149‑1 (MQ=255)
                                ccTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATcgccggc                                                                                                                   >  5:10310/1‑147 (MQ=255)
                                    ctcGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACtt                                                                                                               >  5:388121/1‑149 (MQ=255)
                                    ctcGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACtt                                                                                                               >  2:363250/1‑149 (MQ=255)
                                     tcGGTCCTCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACtt                                                                                                               >  6:43698/1‑148 (MQ=255)
                                        gTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGc                                                                                                           <  3:285340/149‑1 (MQ=255)
                                         tcatcaCGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGcc                                                                                                          <  8:149130/149‑1 (MQ=255)
                                           atcaCGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCgg                                                                                                        >  3:298156/1‑149 (MQ=255)
                                            tcaCGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGc                                                                                                       >  7:223064/1‑149 (MQ=255)
                                            tcaCGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGc                                                                                                       >  1:3118/1‑149 (MQ=255)
                                              aCGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGa                                                                                                     <  1:40606/149‑1 (MQ=255)
                                                gTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACa                                                                                                   >  6:314215/1‑149 (MQ=255)
                                                gTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACa                                                                                                   >  7:23762/1‑149 (MQ=255)
                                                  tCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGCCCATGCCATGGACATcg                                                                                                                                        <  2:355532/110‑2 (MQ=255)
                                                  tCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGCCCATGCCATGGACATcg                                                                                                                                        >  1:355532/1‑109 (MQ=255)
                                                    aGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGc                                                                                               <  4:386351/149‑1 (MQ=255)
                                                     gggTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCcatggcat                                                                                                                                                >  3:117416/1‑99 (MQ=255)
                                                     gggTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCcatggcat                                                                                                                                                <  4:117416/99‑1 (MQ=255)
                                                     gggTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCg                                                                                              >  3:160538/1‑149 (MQ=255)
                                                        tCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTtg                                                                                           <  8:352552/149‑1 (MQ=255)
                                                        accAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACAGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGtt                                                                                            >  7:7392/2‑148 (MQ=255)
                                                          cAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGTCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCTCGACTTCAGGCGGCGACATGGCGTtga                                                                                          <  6:185127/148‑1 (MQ=255)
                                                          cAGCCGATCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTtgat                                                                                         <  1:140415/149‑1 (MQ=255)
                                                            gCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTtgatga                                                                                       <  8:369605/149‑1 (MQ=255)
                                                              cGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGAc                                                                                      >  6:233131/1‑148 (MQ=255)
                                                               ggTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACg                                                                                     >  7:202411/1‑148 (MQ=255)
                                                                gTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCCCGACTTCAGCCGGCGACATGGCGTTGATGACGaa                                                                                   >  3:120990/1‑149 (MQ=255)
                                                                  cAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACt                                                                                 >  1:241466/1‑149 (MQ=255)
                                                                    actgactgGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTCTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTgc                                                                               <  6:161677/149‑1 (MQ=255)
                                                                       cattgGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGcc                                                                            >  1:151759/4‑149 (MQ=255)
                                                                        actgGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACTAACTgc                                                                               <  8:85812/145‑1 (MQ=255)
                                                                        actgGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTgcgc                                                                             <  8:152769/147‑1 (MQ=255)
                                                                        actgGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTgc                                                                               >  7:85812/1‑145 (MQ=255)
                                                                          tgGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCggg                                                                         >  5:250006/1‑149 (MQ=255)
                                                                          tgGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCggg                                                                         >  5:288663/1‑149 (MQ=255)
                                                                           ggCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGt                                                                        >  1:132244/1‑149 (MQ=255)
                                                                            gCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCggg                                                                         <  2:257679/147‑1 (MQ=255)
                                                                             ccAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTa                                                                      >  1:321887/1‑149 (MQ=255)
                                                                             ccAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTa                                                                      <  1:363250/149‑1 (MQ=255)
                                                                              cAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTtgat                                                                                         <  7:152321/129‑1 (MQ=255)
                                                                              cAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTtgat                                                                                         >  8:152321/1‑129 (MQ=255)
                                                                               aGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATc                                                                    <  4:56686/149‑1 (MQ=255)
                                                                                  cGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTc                                                                 <  6:288663/149‑1 (MQ=255)
                                                                                   gAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTcc                                                                <  7:181282/149‑1 (MQ=255)
                                                                                   gAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTcc                                                                >  6:328286/1‑149 (MQ=255)
                                                                                          cTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGa                                                         >  2:24133/1‑149 (MQ=255)
                                                                                             aCCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGAc                                                        <  8:180935/147‑1 (MQ=255)
                                                                                             aCCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGAc                                                        <  4:329041/147‑1 (MQ=255)
                                                                                             aCCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACg                                                       <  4:120990/148‑1 (MQ=255)
                                                                                             aCCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACg                                                       <  4:20212/148‑1 (MQ=255)
                                                                                              ccGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGat                                                     <  2:132244/149‑1 (MQ=255)
                                                                                               cGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGata                                                    <  8:223064/149‑1 (MQ=255)
                                                                                               cGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGata                                                    <  2:241466/149‑1 (MQ=255)
                                                                                               cGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGGCCCAAAGGACGata                                                    <  3:244223/149‑1 (MQ=255)
                                                                                                gCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCTTTGGGACGataa                                                   <  4:218776/149‑2 (MQ=255)
                                                                                                   aCGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCca                                                >  6:90522/1‑149 (MQ=255)
                                                                                                    cGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCcac                                               <  2:321887/149‑1 (MQ=255)
                                                                                                    cGGCAAATCGCCTGGCCCAGCTAGTCCTCGTAGACCGCCATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCCAGACTTCAGCAGCCGACATGGAGTTGATGACGAACTCCGACGTGTTATAGTCACAAAGCACGATATCca                                                >  4:115804/1‑148 (MQ=255)
                                                                                                      gCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCacac                                             <  5:233131/149‑1 (MQ=255)
                                                                                                        cAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACg                                            >  4:222871/1‑148 (MQ=255)
                                                                                                          aTGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTc                                         <  5:43698/149‑1 (MQ=255)
                                                                                                          aTGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTc                                         <  5:314215/149‑1 (MQ=255)
                                                                                                             gCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTcac                                       <  8:23762/148‑1 (MQ=255)
                                                                                                              ccTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTcacca                                     >  8:342873/1‑149 (MQ=255)
                                                                                                              ccTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGAGATCCAGACGCTcacca                                     >  8:126119/1‑149 (MQ=255)
                                                                                                                  ggCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCaccacc                                   >  1:378091/1‑147 (MQ=255)
                                                                                                                  ggCCAGGGTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCaccacc                                   <  4:298156/147‑1 (MQ=255)
                                                                                                                     cAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCt                              <  8:34741/149‑1 (MQ=255)
                                                                                                                        gTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTccc                           <  2:378091/149‑1 (MQ=255)
                                                                                                                         ttGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTcccc                          >  1:148352/1‑149 (MQ=255)
                                                                                                                                 ggcgacggcgaTGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACg                   <  3:265743/148‑1 (MQ=255)
                                                                                                                                   cgacggcgaTGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGcc                <  5:328286/149‑1 (MQ=255)
                                                                                                                                   cgacggcgaTGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGcc                <  2:148352/149‑1 (MQ=255)
                                                                                                                                    gacggcgaTGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCca                                                >  8:81514/1‑116 (MQ=255)
                                                                                                                                    gacggcgaTGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCca                                                <  7:81514/116‑1 (MQ=255)
                                                                                                                                    gacggcgaTGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCt               >  6:350695/1‑149 (MQ=255)
                                                                                                                                    gacggcgaTGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCt               <  6:63399/149‑1 (MQ=255)
                                                                                                                                     atgaagaTGCCCATGCCATGGACATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCTg              >  5:340425/6‑149 (MQ=255)
                                                                                                                                              tCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCTGGAcgcgcg      <  2:3118/148‑1 (MQ=255)
                                                                                                                                                catggcatggGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCTGGACGCGCGAAc   >  6:232007/1‑149 (MQ=255)
                                                                                                                                                 atggcatggGCATCTTCATCAACGATGATCGCCGCTACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCTGGACGCGCGAAcc  <  7:386663/149‑1 (MQ=255)
                                                                                                                                                 atggcatggGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCTGGACGCGCGAAc   >  1:367811/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACCGGTTTCGGCCTGTTGCTTGGCCTGGATGTCCTTCTCGGTCATCACGTTCAGGGTCCAGCCGGTCAACTGACTGGCCAGGCGAACGTTCTGACCGCCACGGCCAATGGCCTGGGCCAGGTTGTCCTCGGCGACGGCGATGTCCATGGCATGGGCATCTTCATCAACGATGATCGCCGCGACTTCAGCCGGCGACATGGCGTTGATGACGAACTGCGCCGGGTTATCGTCCCAAAGGACGATATCCACACGCTCACCACCCAGCTCCCCGGATACGGCCTGGACGCGCGAACC  >  NC_002947/5359205‑5359498

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