Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,591,705 G→T 7.3% A63A (GCC→GCA PP_1395 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,591,7050GT 7.3% 203.5 / 4.4 64A63A (GCC→GCAPP_1395AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base T (3/2);  total (36/28)
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.31e-01

CAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTC  >  NC_002947/1591558‑1591845
                                                                                                                                                   |                                                                                                                                            
caggtccaTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTggg                                                                                                                                             >  6:213334/1‑149 (MQ=255)
   gtccaTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCaa                                                                                                                                          <  5:178793/149‑1 (MQ=255)
   gtccaTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCaa                                                                                                                                          <  7:357355/149‑1 (MQ=255)
      caTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAgg                                                                                                                                        <  8:34762/148‑1 (MQ=255)
      caTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAgg                                                                                                                                        <  5:6266/148‑1 (MQ=255)
      caTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGa                                                                                                                                       <  7:236339/149‑1 (MQ=255)
          cgcAGGCCCAGGCTTGGGGTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGca                                                                                                                                   <  8:169864/149‑1 (MQ=255)
                 ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAggg                                                                                                                              >  2:26723/1‑147 (MQ=255)
                 ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAggg                                                                                                                              >  6:203034/1‑147 (MQ=255)
                        tGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAg                                                                                                                                         >  1:190257/1‑129 (MQ=255)
                        tGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAg                                                                                                                                         <  2:190257/129‑1 (MQ=255)
                           gTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTgcgc                                                                                                                    >  3:189990/1‑147 (MQ=255)
                               tcGGGGCTGGCATCTAGCCTCAAGGCGGTGAGTGGCAATGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                              <  8:340289/149‑1 (MQ=255)
                               tcGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                              >  6:281251/1‑149 (MQ=255)
                                 ggggCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                              <  3:232783/147‑1 (MQ=255)
                                  gggCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggtggg                                                                                                           <  1:26723/149‑1 (MQ=255)
                                     cTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTa                                                                                                         >  8:63667/1‑148 (MQ=255)
                                        gCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACa                                                                                                       <  5:203034/147‑1 (MQ=255)
                                                 cTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCaa                                                                                                                                          >  4:131528/1‑103 (MQ=255)
                                                 cTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCaa                                                                                                                                          <  3:131528/103‑1 (MQ=255)
                                                     cGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGAt                                                                                                                                      <  1:82256/103‑1 (MQ=255)
                                                     cGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGAt                                                                                                                                      >  2:82256/1‑103 (MQ=255)
                                                     cGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCg                                                                                        <  6:402526/149‑1 (MQ=255)
                                                       gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                      <  7:63667/149‑1 (MQ=255)
                                                       gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                      <  7:145248/149‑1 (MQ=255)
                                                           tGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTaa                                                                                  >  7:356218/1‑149 (MQ=255)
                                                                 gCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCAc                                                                             >  1:402237/1‑148 (MQ=255)
                                                                  cAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCa                                                                           <  4:24119/149‑1 (MQ=255)
                                                                         aCCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAg                                                                    >  7:308672/1‑149 (MQ=255)
                                                                           cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAgg                                                                         >  3:224822/1‑142 (MQ=255)
                                                                           cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAgg                                                                         <  4:224822/142‑1 (MQ=255)
                                                                           cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAgg                                                                   >  5:404030/1‑148 (MQ=255)
                                                                           cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGt                                                                  >  4:198308/1‑149 (MQ=255)
                                                                                   ccATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                          >  1:168515/1‑149 (MQ=255)
                                                                                   ccATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                          <  5:281251/149‑1 (MQ=255)
                                                                                                tGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGAc                                             <  8:356218/149‑1 (MQ=255)
                                                                                                tGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCTTGGCAAAGAAGCACAAGGCAGGCAGCGCCAGGGTTGGTACAGAAAGTACATTTTAGGTGTTACGCACCAGGGACAAGTCATAGAATGCCCAAGCTTGAc                                             >  2:136490/1‑148 (MQ=255)
                                                                                                    tCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCa                                           >  2:93873/1‑147 (MQ=255)
                                                                                                    tCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCa                                           >  6:277304/1‑147 (MQ=255)
                                                                                                        aGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGc                                     >  8:327011/1‑149 (MQ=255)
                                                                                                          gCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTc                                   >  2:229806/1‑149 (MQ=255)
                                                                                                               cTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgg                                                                                                                >  4:265967/1‑67 (MQ=255)
                                                                                                               cTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgg                                                                                                                <  3:265967/67‑1 (MQ=255)
                                                                                                                   tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                          >  1:278900/1‑117 (MQ=255)
                                                                                                                   tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                          <  2:278900/117‑1 (MQ=255)
                                                                                                                   tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATg                          >  2:32615/1‑149 (MQ=255)
                                                                                                                   tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATg                          >  3:95082/1‑149 (MQ=255)
                                                                                                                    cGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATgg                         >  1:358339/1‑149 (MQ=255)
                                                                                                                     gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTAc                                                                                                        >  2:328632/1‑69 (MQ=255)
                                                                                                                     gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTAc                                                                                                        <  1:328632/69‑1 (MQ=255)
                                                                                                                     gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGCCCTAAGGGCTCGCCCGCAAGGc                        >  8:384333/1‑149 (MQ=255)
                                                                                                                     gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGc                        <  5:213334/149‑1 (MQ=255)
                                                                                                                        aCCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCgg                      <  4:189990/148‑1 (MQ=255)
                                                                                                                          ccAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCACCCGCATGGCGGtt                    >  1:236927/1‑148 (MQ=255)
                                                                                                                              gAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCa               <  2:168515/149‑1 (MQ=255)
                                                                                                                                  aGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTc           >  6:204831/1‑149 (MQ=255)
                                                                                                                                    gcgcagcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCgg         >  6:335432/8‑149 (MQ=38)
                                                                                                                                    gCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCgg         >  4:1846/1‑149 (MQ=255)
                                                                                                                                     cgcagcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGc        <  5:335432/143‑1 (MQ=255)
                                                                                                                                      ttttGCTGCCCTGGGCAAGGATGCACAGGGCAGGCAGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGc        >  1:268745/1‑148 (MQ=255)
                                                                                                                                        agcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAg     >  6:172323/4‑149 (MQ=255)
                                                                                                                                          gctgcCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCt   >  7:29436/1‑149 (MQ=255)
                                                                                                                                          cctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTg                                                >  7:200816/2‑104 (MQ=38)
                                                                                                                                          cctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTg                                                <  8:200816/103‑1 (MQ=38)
                                                                                                                                            tgcCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTc  >  6:75856/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTC  >  NC_002947/1591558‑1591845

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: