Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,658,123 A→G 6.7% K60E (AAG→GAG)  sahR → methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,658,1230AG 6.7% 170.9 / 2.0 59K60E (AAG→GAG) sahRmethionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base G (2/2);  total (32/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.90e-01

ACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCCGA  >  NC_002947/5657978‑5658271
                                                                                                                                                 |                                                                                                                                                    
aCGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAgg                                                                                                                                                   >  5:23474/1‑149 (MQ=255)
 cgcgAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGt                                                                                                                                                  >  6:177322/1‑149 (MQ=255)
  gcgAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTCGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTg                                                                                                                                                 >  1:24971/1‑149 (MQ=255)
    gAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCt                                                                                                                                               <  8:19766/149‑1 (MQ=255)
         aTTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCt                                                                                                                                           >  1:222491/1‑148 (MQ=255)
          ttGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCt                                                                                                                                           >  6:206818/1‑147 (MQ=255)
          ttGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCt                                                                                                                                           >  6:258457/1‑147 (MQ=255)
            ggCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAg                                                                                                                                        <  7:143529/148‑1 (MQ=255)
            ggCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTC‑GGCCTGAGCCAGCACCTTGAGGTGATGGCTCatg                                                                                                                                       <  6:313571/148‑3 (MQ=38)
             gCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAgg                                                                                                                                       <  8:51856/148‑1 (MQ=255)
             gCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTC‑GGCCTGAGCCAGCACCTTGAGGTGATGGCTCATGcc                                                                                                                                     <  6:149456/149‑1 (MQ=38)
                gcgcTGTGCAAAGCGAGTGGCTATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCATGGTGCTGGCTCAGGCCGa                                                                                                                                   <  8:104285/146‑1 (MQ=255)
                        cAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGGTCAGGCCGACCTCGTgg                                                                                                                           <  2:24971/149‑1 (MQ=255)
                        cAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTg                                                                                                                            <  1:310141/148‑1 (MQ=255)
                             ccAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTCGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACGACGTGGCAAcc                                                                                                                      >  6:193531/1‑149 (MQ=255)
                                         cgctgcgcCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAgg                                                                                                           <  5:177322/148‑1 (MQ=255)
                                          gctgcgcCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGc                                                                                                          >  5:114730/1‑148 (MQ=255)
                                                   gAACGTTTTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCAtct                                                                                                >  4:318209/1‑149 (MQ=255)
                                                    aaCGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCAtctt                                                                                               <  7:236128/149‑1 (MQ=255)
                                                    aaCGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCAtctt                                                                                               <  6:114730/149‑1 (MQ=255)
                                                             cgcgcCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTAccgcc                                                                                       <  3:194534/148‑1 (MQ=255)
                                                                      gCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGtgccctgcc                                                                             >  8:106695/1‑149 (MQ=255)
                                                                        cAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCaa                                                                                                                        >  2:239832/1‑104 (MQ=255)
                                                                        cAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCaa                                                                                                                        <  1:239832/104‑1 (MQ=255)
                                                                         aaTGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGAGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGa                                                                          >  3:216890/1‑149 (MQ=255)
                                                                           tGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTTGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATg                                                                        <  6:23474/149‑1 (MQ=255)
                                                                                 gTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTgc                                                                  >  2:276172/1‑149 (MQ=255)
                                                                                    cGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGTGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTgcgc                                                                <  6:42797/148‑1 (MQ=255)
                                                                                      gCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTg                                                             <  5:89577/149‑1 (MQ=255)
                                                                                         tGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGt                                                          >  4:88038/1‑149 (MQ=255)
                                                                                                 cTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCt                                                  <  2:282299/149‑1 (MQ=255)
                                                                                                   ggCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTg                                                 >  2:176303/1‑148 (MQ=255)
                                                                                                    gcgcAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCa                                               >  6:38650/1‑149 (MQ=255)
                                                                                                      gcAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCAcg                                             >  2:122073/1‑149 (MQ=255)
                                                                                                      gcAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCAcg                                             >  2:103597/1‑149 (MQ=255)
                                                                                                             ttCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCgct                                      >  1:123280/1‑149 (MQ=255)
                                                                                                             ttCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCgct                                      >  6:243309/1‑149 (MQ=255)
                                                                                                                    tCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaa                               <  7:106695/149‑1 (MQ=255)
                                                                                                                     cGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATg                                                                        <  8:139138/107‑1 (MQ=255)
                                                                                                                     cGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATg                                                                        >  7:139138/1‑107 (MQ=255)
                                                                                                                       gCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCACCGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaag                              >  5:232023/1‑147 (MQ=255)
                                                                                                                       gCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaagaa                            <  1:122073/149‑1 (MQ=255)
                                                                                                                       gCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaagaa                            <  5:243309/149‑1 (MQ=255)
                                                                                                                       gCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaag                              >  8:83502/1‑147 (MQ=255)
                                                                                                                       gCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaag                              >  6:300422/1‑147 (MQ=255)
                                                                                                                       gCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaag                              >  4:15918/1‑147 (MQ=255)
                                                                                                                       gCCAGTCGGGCAGGAGCCATCACCTCAAGGGGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCgaagaa                            <  2:214751/149‑1 (MQ=255)
                                                                                                                         cAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGt                          <  7:83502/149‑1 (MQ=255)
                                                                                                                          aGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTc                         <  5:258457/149‑1 (MQ=255)
                                                                                                                            tCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTcga                       <  5:38650/149‑1 (MQ=255)
                                                                                                                            tCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTcga                       <  3:88038/149‑1 (MQ=255)
                                                                                                                                       ccAGCACCTTGAGGTGATGGCTCATGCCGACCTCGTGGCACCCCGCCGCGAAGGCACCGGCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCGGGCTGGGCGGCTGCACGCGGCGCGGCTCGATGAAGTCGACGAGCTGGccc            >  8:178361/1‑149 (MQ=255)
                                                                                                                                       ccAGCACCTTGAGGTGATGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGccc            >  1:297355/1‑149 (MQ=255)
                                                                                                                                         aTCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                            >  2:272119/1‑99 (MQ=255)
                                                                                                                                         aTCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTgg                                                            <  1:272119/99‑1 (MQ=255)
                                                                                                                                         aTCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTa          >  4:176952/1‑149 (MQ=255)
                                                                                                                                          tCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTAc         >  3:134907/1‑149 (MQ=255)
                                                                                                                                              cTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGcc     >  4:148038/1‑149 (MQ=255)
                                                                                                                                              cTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTCGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGcc     >  8:3365/1‑149 (MQ=255)
                                                                                                                                                 aaGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCCg   <  6:232023/148‑1 (MQ=255)
                                                                                                                                                 aaGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCCGa  <  2:222491/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCAGTCGGGCATGAGCCATCACCTCAAGGTGCTGGCTCAGGCCGACCTCGTGGCAACCCGCCGCGAAGGCAACGCCATCTTCTACCGCCGTGCCCTGCCCGATGGCCTGCGCCTGGGTGGCCGGCTGCACGCGGCGCTGCTCGAAGAAGTCGACGACCTGGCCCTACCGCCCGA  >  NC_002947/5657978‑5658271

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: ‑