Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 908,703 T→A 6.3% Q169L (CAG→CTG)  cra ← catabolite DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947908,7030TA 6.3% 181.9 / 2.6 64Q169L (CAG→CTG) cracatabolite DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base A (2/2);  total (32/32)
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.63e-01

GCCCTGATAGCGGCGCACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAA  >  NC_002947/908558‑908845
                                                                                                                                                 |                                                                                                                                              
gCCCTGATAGCGGCGCACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTgg                                                                                                                                              <  1:366877/148‑1 (MQ=255)
  ccTGATAGCGGCGCACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTggcgg                                                                                                                                           >  8:158456/1‑149 (MQ=255)
     gATAGCGGCGCACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTg                                                                                                                                        >  2:114307/1‑149 (MQ=255)
      aTAGCGGCGCACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTgg                                                                                                                                       >  6:154692/1‑149 (MQ=255)
         gcggcgCACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCGCGGGGCGTGCGCCTATCAGGGCGATGCGGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCCGCCAGCTGGcc                                                                                                                                     <  3:32007/148‑2 (MQ=255)
             cgcACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATcgcg                                                                                                                                <  8:316771/149‑1 (MQ=255)
             cgcACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATcgcg                                                                                                                                <  1:114307/149‑1 (MQ=255)
             cgcACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATcgcg                                                                                                                                <  2:171447/149‑1 (MQ=255)
                  ttCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGgtcgt                                                                                                                           >  1:15510/1‑149 (MQ=255)
                       cGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCAGGCGGCTGGCATCGCGGTCGTCGCt                                                                                                                       >  7:11373/1‑148 (MQ=255)
                       cGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCAGGCGGCTGGCATCGCGGTCGTCGCt                                                                                                                       >  7:16793/1‑148 (MQ=255)
                        ggTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGa                                                                                                                     >  6:331008/1‑149 (MQ=255)
                          tataGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATc                                                                                                                   <  5:304256/149‑1 (MQ=255)
                               ccTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCgg                                                                                                                               >  6:197870/1‑132 (MQ=255)
                               ccTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCgg                                                                                                                               <  5:197870/132‑1 (MQ=255)
                               ccTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCAcc                                                                                                                >  6:67092/1‑147 (MQ=255)
                               ccTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCAcc                                                                                                                >  3:362567/1‑147 (MQ=255)
                                 ttGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCAAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAg                                                                                                             >  4:19363/1‑148 (MQ=255)
                                         cTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAgtgt                                                                                                    <  2:15510/149‑1 (MQ=255)
                                          ttCGTCGAATCAGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCAATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGtgtg                                                                                                   >  2:247022/1‑149 (MQ=255)
                                              tcgAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCggg                                                                                               >  2:307854/1‑149 (MQ=255)
                                              tcgAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCggg                                                                                               <  1:144444/149‑1 (MQ=255)
                                              tcgAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCggg                                                                                               >  3:326206/1‑149 (MQ=255)
                                                gAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGt                                                                                             <  4:201010/149‑1 (MQ=255)
                                                        cgGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGAcgg                                                                                     >  5:106920/1‑149 (MQ=255)
                                                            aCGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGGGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGatcgga                                                                                   >  1:121817/1‑142 (MQ=255)
                                                            aCGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCgg                                                                                                >  5:318224/1‑134 (MQ=255)
                                                            aCGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCgg                                                                                                <  6:318224/134‑1 (MQ=255)
                                                            aCGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGGCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTc                                                                                            <  2:121817/138‑1 (MQ=255)
                                                                gcCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAg                                                                                                          <  7:245284/120‑1 (MQ=255)
                                                                gcCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAg                                                                                                          >  8:245284/1‑120 (MQ=255)
                                                                gcCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGAt                                                                             <  4:170726/149‑1 (MQ=255)
                                                                 ccTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATg                                                                            <  3:287567/149‑1 (MQ=255)
                                                                     gCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCTGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCaa                                                                        >  4:202960/1‑149 (MQ=255)
                                                                       tAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAAt                                                                       >  4:165180/1‑148 (MQ=255)
                                                                           cGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCg                                                                  >  2:318077/1‑149 (MQ=255)
                                                                             aCAGCTGGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGc                                                                >  6:252133/1‑149 (MQ=255)
                                                                             aCAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCt                                                                                                                                                <  8:220329/69‑1 (MQ=255)
                                                                             aCAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCt                                                                                                                                                >  7:220329/1‑69 (MQ=255)
                                                                             aCAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCgg                                                                 >  5:230601/1‑148 (MQ=255)
                                                                                gCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAgg                                                             <  5:154692/149‑1 (MQ=255)
                                                                                gCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAgg                                                             <  3:159201/149‑1 (MQ=255)
                                                                                gCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCGGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAgg                                                             <  4:338934/149‑1 (MQ=255)
                                                                                 cTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGc                                                            <  4:362567/149‑1 (MQ=255)
                                                                                     gggCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCttt                                                        <  1:157025/149‑1 (MQ=255)
                                                                                      ggCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACTGCGGTCGATGGCAATCACCGGCAGGCGTTTg                                                       <  1:247022/149‑1 (MQ=255)
                                                                                      ggCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTg                                                       <  7:158456/149‑1 (MQ=255)
                                                                                      ggCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCAGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGATCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTg                                                       <  8:16793/149‑1 (MQ=255)
                                                                                      ggCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCAGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGATCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTg                                                       <  8:11373/149‑1 (MQ=255)
                                                                                            cgcCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTg                                                  <  5:67092/148‑1 (MQ=255)
                                                                                                 aTCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGCGCAGAAGTGTGCGGGGTCCAGCCGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTc                                            >  8:158/1‑149 (MQ=255)
                                                                                                    aGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTcgcg                                         >  7:341735/1‑149 (MQ=255)
                                                                                                     gggCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGTTGAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCgg                                        <  1:307854/149‑1 (MQ=255)
                                                                                                     gggCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCgg                                        >  5:174244/1‑149 (MQ=255)
                                                                                                        cGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAg                                     <  1:318077/149‑1 (MQ=255)
                                                                                                                cGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCa                               <  3:90666/147‑1 (MQ=255)
                                                                                                                            cTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCggcag                 <  3:19363/149‑1 (MQ=255)
                                                                                                                              gAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCggcagg                >  8:318315/1‑148 (MQ=255)
                                                                                                                                gCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCaa             >  2:277501/1‑149 (MQ=255)
                                                                                                                                    gCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTgg         >  2:123078/1‑149 (MQ=255)
                                                                                                                                    gCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTgg         <  3:202960/149‑1 (MQ=255)
                                                                                                                                    gCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTgg         <  8:341735/149‑1 (MQ=255)
                                                                                                                                    gCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTgg         >  8:79509/1‑149 (MQ=255)
                                                                                                                                      gcgatgcCAGCCGCCAGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGc        >  7:238703/6‑148 (MQ=255)
                                                                                                                                        gcggcCAGCTGGCGGCAGGCATCGCCGTCGTCGCTGCTCAAAGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCAACGGGAGGCCTGTGTCCTGCAGCTAGCGCTAGCTGTCATCTTCTGGCGGCaggcaacaggcaac     >  6:177530/1‑149 (MQ=255)
                                                                                                                                           gcCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGaa  <  4:326206/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GCCCTGATAGCGGCGCACTTCGCCGGTATAGCCTTGCAGGGCTTCGTCGAATCCGCCGGCACGCGCCTGGCTAACCGACAGCTCGGGGCGTGCGCCTATCAGGGCGATGCTGCGTGGGGCCGAGCTGAGCAGGCTGGCGGCCAGCTGGCGGCTGGCATCGCGGTCGTCGCTGATCACCGAGCAGAAGTGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAA  >  NC_002947/908558‑908845

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