Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 359,700:1 +T 5.3% intergenic (+60/+60) 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,7001.T 5.3% 320.0 / 8.5 75intergenic (+60/+60)gbdR/PP_0299AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (33/38);  new base T (2/2);  total (35/40)
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.59e-01

CATGACCCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCAC  >  NC_002947/359556‑359842
                                                                                                                                                  |                                                                                                                                               
cATGACCCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTCAGCCAGCCGCGCAATGAGACGACATTTGCCAGTCTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTACCACAGGACACCCGCAGAGCTTGAGCCAGc                                                                                                                                             >  1:233843/1‑149 (MQ=255)
 aTGACCCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCg                                                                                                                                           <  4:222433/149‑1 (MQ=255)
 aTGACCCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCg                                                                                                                                           >  3:227785/1‑149 (MQ=255)
 aTGACCCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCg                                                                                                                                           >  3:103394/1‑149 (MQ=255)
      ccTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGAt                                                                                                                                      <  4:227785/149‑1 (MQ=255)
      ccTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGAt                                                                                                                                      >  3:79044/1‑149 (MQ=255)
      ccTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGAt                                                                                                                                      >  3:212076/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                                                                                                                         >  1:189839/1‑120 (MQ=255)
                                                gcgcAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGgggcggg                                                                                                                         <  2:189839/120‑1 (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                                                                                                     <  6:87213/129‑1 (MQ=255)
                                                           cGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCg                                                                                                     >  5:87213/1‑129 (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                                                          <  5:95768/147‑1 (MQ=255)
                                                                                    acacCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCaa                                                          <  4:317922/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                                                                                              >  3:168269/1‑106 (MQ=255)
                                                                                         gcgcCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTgg                                                                                              <  4:168269/106‑1 (MQ=255)
                                                                                           gcCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCaa                                                 <  7:10064/149‑1 (MQ=255)
                                                                                           gcCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCaa                                                 <  4:79044/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                           >  5:23489/1‑149 (MQ=255)
                                                                                                                  ccGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTgg                           >  1:148411/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)
                                                                                                                                                  |                                                                                                                                               
CATGACCCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCAC  >  NC_002947/359556‑359842

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: