Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 359,704:1 +G 5.6% intergenic (+64/+56) 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,7041.G 5.6% 303.1 / 9.0 72intergenic (+64/+56)gbdR/PP_0299AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (30/38);  new base G (2/2);  total (32/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 = 1.88e-01

CATGACCCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCACTGCCGGTGT  >  NC_002947/359556‑359851
                                                                                                                                                       |                                                                                                                                                   
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:314575/148‑1 (MQ=255)
            gccgcACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCtgtg                                                                                                                                          <  5:214818/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:380156/149‑1 (MQ=255)
                                                                     cAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTttgt                                                                                <  4:374117/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                                                                                                       >  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                                    >  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)
                                                                                                                                                       |                                                                                                                                                   
CATGACCCTGTCGCCGCACAGCGGGCCGATGGCGGCCCTGAGCCAGGCGCGCAATGAGTCGACATTTGCCAGTGTAAGGCTCTGACACCGCGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTCACCCGCTGAGCTTGAGCCAGCGGTGATCTGTGGGGGCGGGCTTGTTCCGCGAAGGGGCCGGTCCTGGGTAAATGGATTACGCGCGTTTGTTGAACTGCGCCAACGCTGGCAACAGTTGCTTGTCGATGGCCTGGCGTACTGCCGGCAGGATCGTCGCACTGCCGGTGT  >  NC_002947/359556‑359851

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: