Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 359,710 T→C 5.5% intergenic (+70/+50) gbdR → / ← PP_0299 AraC family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947359,7100TC 5.5% 225.2 / 2.2 73intergenic (+70/+50)gbdR/PP_0299AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/40);  new base C (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

CCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCACTGCCGGTGTACA  >  NC_002947/359562‑359854
                                                                                                                                                       |                                                                                                                                                
ccTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGAt                                                                                                                                                  <  4:227785/149‑1 (MQ=255)
ccTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGAt                                                                                                                                                  >  3:212076/1‑149 (MQ=255)
ccTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGAt                                                                                                                                                  >  3:79044/1‑149 (MQ=255)
    tCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCtgt                                                                                                                                              <  7:258531/149‑1 (MQ=255)
      gccgcACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCtgtg                                                                                                                                             <  5:214818/148‑1 (MQ=255)
      gccgcACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCtgtg                                                                                                                                             <  5:314575/148‑1 (MQ=255)
      gccgcACAGCGGGCCGATGGCGGCCCGGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTgg                                                                                                                                            <  7:343944/149‑1 (MQ=255)
        cgcACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTgggg                                                                                                                                          <  8:221523/149‑1 (MQ=255)
         gcACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTggggg                                                                                                                                         >  2:337765/1‑149 (MQ=255)
          cacaGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGgggc                                                                                                                                        <  6:69798/149‑1 (MQ=255)
            caGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGAACAAGCCCGCCCCCACAGATCACCGCTGGCTCAAGCTCAGCGGGTGACCTGTGGgagcg                                                                                                                                       <  5:365246/149‑1 (MQ=17)
            caGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCCCCCACAGATCACCGCTGGCTCAAGCTCAGCGGGTGACCTGTGGgagcgg                                                                                                                                      <  2:127805/149‑1 (MQ=21)
                gggCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCtt                                                                                                                                  >  3:62338/1‑149 (MQ=255)
                     gATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTcc                                                                                                                             >  6:189784/1‑149 (MQ=255)
                         gcggcCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGa                                                                                                                         <  3:282321/149‑1 (MQ=255)
                          cggcCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGaa                                                                                                                        <  7:11373/149‑1 (MQ=255)
                            gcCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAgg                                                                                                                      <  2:5342/149‑1 (MQ=255)
                             cccTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAggg                                                                                                                     <  2:203486/149‑1 (MQ=255)
                               cTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGc                                                                                                                   <  2:225677/149‑1 (MQ=255)
                                 gAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCg                                                                                                                 <  7:227624/149‑1 (MQ=255)
                                     cAGGCGCGCAATGAGTCGACATTTCCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTCACCAGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTAGGTCCGCGAAGGGGCCGGt                                                                                                               >  6:173538/1‑147 (MQ=39)
                                        gcgcgcAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTgg                                                                                                          <  6:96983/149‑1 (MQ=255)
                                         cgctcAATGAGTCGACATTTGCCTGTGTAAGGCTCTGACACGGGGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTggg                                                                                                         <  2:233843/149‑1 (MQ=255)
                                          gcgcAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGgggcggg                                                                                                                                     <  2:189839/120‑1 (MQ=255)
                                          gcgcAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGgggcggg                                                                                                                                     >  1:189839/1‑120 (MQ=255)
                                           cgcAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGt                                                                                                        >  3:305312/1‑148 (MQ=255)
                                                  aGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGAt                                                                                                <  4:103394/149‑1 (MQ=255)
                                                   gTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGTCGGTCCTGGGTAAATGGAtt                                                                                               <  5:189784/149‑1 (MQ=255)
                                                    tCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTa                                                                                              >  8:10064/1‑149 (MQ=255)
                                                     cGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCg                                                                                                                 >  5:87213/1‑129 (MQ=255)
                                                     cGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCg                                                                                                                 <  6:87213/129‑1 (MQ=255)
                                                         aTTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTAcgcgc                                                                                         <  5:1836/149‑1 (MQ=255)
                                                           ttGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGt                                                                                       >  3:126398/1‑149 (MQ=255)
                                                              ccAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTttg                                                                                    >  4:98884/1‑149 (MQ=255)
                                                               cAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCGCACAGGTCACCCGCTGAGCTTGAGCCTGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTttgt                                                                                   <  5:5777/149‑1 (MQ=255)
                                                               cAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTttgt                                                                                   <  4:374117/149‑1 (MQ=255)
                                                               cAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTttgt                                                                                   <  4:380156/149‑1 (MQ=255)
                                                                aGTGTACGGCTCTGACACCGCGCCGGCCTAGTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTttgtt                                                                                  <  7:321887/149‑1 (MQ=255)
                                                                aGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTttgtt                                                                                  <  3:19008/149‑1 (MQ=255)
                                                                    tAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGaa                                                                               <  4:305312/148‑1 (MQ=255)
                                                                       ggCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTgc                                                                           <  4:126398/149‑1 (MQ=255)
                                                                          tctGACACCGCGCCGGCCTCTTCGCGGGCACGCCCCCAACCCCAAGTCACCCCCGGGCCGTGAGCCACCGGGGATCAGTTGGGGCGGGGTTTTTCCGCGAAGGGGGCGGGACGGGGGAAATTGATTTAGCGCGTTTTTTGTACAGAGcc                                                                        >  2:193072/1‑149 (MQ=255)
                                                                            tGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCaa                                                                      >  6:95768/1‑149 (MQ=255)
                                                                            tGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGGAAATGGATTACGCGCGTTTGTTGAACTGCGCCaa                                                                      >  1:62651/1‑149 (MQ=255)
                                                                              acacCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCaa                                                                      <  4:317922/147‑1 (MQ=255)
                                                                              acacCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCaa                                                                      <  5:95768/147‑1 (MQ=255)
                                                                                acCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCt                                                                  >  3:147221/1‑149 (MQ=255)
                                                                                 ccGCGCCGGTCTCTTCGCGGGCACGCGCGCTCCCACAGGTAACCCGCTGAGGTTGAGCTCGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGATGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCATCGCTg                                                                 <  5:173538/149‑1 (MQ=255)
                                                                                 ccGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTg                                                                 >  4:324462/1‑149 (MQ=255)
                                                                                 ccGCGCCCGCCTCTTCGCGGGCACGCCCGCTCCCACTCGTCACCCGCTGTGCTTGCGCCCGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTg                                                                 <  2:62651/149‑1 (MQ=255)
                                                                                  cgcgCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTgg                                                                <  3:198811/149‑1 (MQ=255)
                                                                                   gcgcCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTgg                                                                                                          <  4:168269/106‑1 (MQ=255)
                                                                                   gcgcCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTgg                                                                                                          >  3:168269/1‑106 (MQ=255)
                                                                                     gcCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCaa                                                             <  4:79044/149‑1 (MQ=255)
                                                                                     gcCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCaa                                                             <  7:10064/149‑1 (MQ=255)
                                                                                          ccTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGtt                                                        >  8:204650/1‑149 (MQ=255)
                                                                                             ctTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTAcgc                                                                                           >  8:242342/1‑111 (MQ=255)
                                                                                                cgcgGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTGGAACTGCGCCAACGGTGGCCCCCGTTGCTGGt                                                  >  6:266394/1‑149 (MQ=255)
                                                                                                   gggCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGa                                               >  5:9606/1‑149 (MQ=255)
                                                                                                        aCGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGc                                           >  6:83485/1‑149 (MQ=255)
                                                                                                         cGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGcc                                          >  2:66014/1‑149 (MQ=255)
                                                                                                            ccGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTgg                                       >  1:148411/1‑149 (MQ=255)
                                                                                                            ccGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTgg                                       >  5:23489/1‑149 (MQ=255)
                                                                                                            ccGCCCCCACAGATCACCGCTGGCTCAAGCTCAGCGGGTGACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGc                                      >  8:149523/1‑149 (MQ=14)
                                                                                                             cGCCCCCACAGATCACCGCTGGCTCAAGCTCAGCGGGTGACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGCACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCg                                     >  1:177506/1‑149 (MQ=14)
                                                                                                                        gTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCgg                            >  8:308424/1‑148 (MQ=255)
                                                                                                                        gTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGc                           <  1:337765/149‑1 (MQ=255)
                                                                                                                         tCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCgg                            >  3:1371/1‑147 (MQ=255)
                                                                                                                           aCCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAg                         >  2:151912/1‑148 (MQ=255)
                                                                                                                                 ggAGCTTGCGCCAGCGGTGATCTGTGGGGGCGGGCGTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGAtcgtc                  <  7:336817/148‑1 (MQ=255)
                                                                                                                                     cTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCAGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCAc              >  8:302960/1‑149 (MQ=255)
                                                                                                                                               gCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCACTGCCGgtgt     <  3:324462/149‑1 (MQ=255)
                                                                                                                                                  ggTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCACTGCCGGTGTAc   <  6:9606/149‑1 (MQ=255)
                                                                                                                                                   gTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCACTGCCGGTGTACa  <  5:83485/149‑1 (MQ=255)
                                                                                                                                                      aTCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTg                                        >  2:103681/1‑108 (MQ=255)
                                                                                                                                                      aTCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTg                                        <  1:103681/108‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
CCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCACTGCCGGTGTACA  >  NC_002947/359562‑359854

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