Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,974,392 G→A 6.2% R27R (CGC→CGT PP_3503 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,974,3920GA 6.2% 188.3 / 2.3 64R27R (CGC→CGTPP_3503sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

CAATCCATTCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCGTTTAAC  >  NC_002947/3974244‑3974538
                                                                                                                                                      |                                                                                                                                                  
cAATCCATTCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTg                                                                                                                                                    >  8:169478/1‑149 (MQ=255)
  aTCCATTCGGTGTTGCTGCGCGTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTgcg                                                                                                                                                  <  2:323565/149‑1 (MQ=255)
       ttCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGc                                                                                                                                             <  2:179303/149‑1 (MQ=255)
          ggTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAg                                                                                                                                          <  5:327378/149‑1 (MQ=255)
           gtgtTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAgg                                                                                                                                         >  4:158518/1‑149 (MQ=255)
              ttGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTa                                                                                                                                       <  1:239956/149‑1 (MQ=255)
                     gcTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCgg                                                                                                                               <  6:183524/149‑1 (MQ=255)
                      cTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGt                                                                                                                              >  7:155232/1‑149 (MQ=255)
                          aTCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCACGCGCAATAGCCCGACATCGCACGTGTGGTCCCGCGCAGCCCCGAGCATACAGCTGTCCACATCCCCGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGc                                                                                                                           >  2:297989/1‑148 (MQ=255)
                          aTCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGc                                                                                                                           >  2:293541/1‑148 (MQ=255)
                               gTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAAtcgt                                                                                                                     >  7:27012/1‑149 (MQ=255)
                                  cTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAAtcgtcg                                                                                                                   >  6:268869/1‑148 (MQ=255)
                                  cTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAAtcgtcg                                                                                                                   >  4:79123/1‑148 (MQ=255)
                                  cTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAAtcgtcg                                                                                                                   >  7:258679/1‑148 (MQ=255)
                                      aCCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTg                                                                                                                             <  5:194236/134‑1 (MQ=255)
                                      aCCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTg                                                                                                                             >  6:194236/1‑134 (MQ=255)
                                      aCCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAggg                                                                                                              <  8:258679/149‑1 (MQ=255)
                                        cgcgTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGt                                                                                                            >  7:75486/1‑149 (MQ=255)
                                           gTCCGGGTGCCTCAGGTGCGAAGTCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTcga                                                                                                         <  1:126975/149‑1 (MQ=255)
                                                 gTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGAt                                                                                                     >  6:100791/1‑147 (MQ=255)
                                                 gTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATaa                                                                                                   <  7:42481/149‑1 (MQ=255)
                                                   gCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAg                                                                                                  >  6:9335/1‑148 (MQ=255)
                                                      tCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATaa                                                                                                   <  7:70231/144‑1 (MQ=255)
                                                      tCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATaa                                                                                                   >  8:70231/1‑144 (MQ=255)
                                                      tCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCaaa                                                                                              >  7:193299/1‑149 (MQ=255)
                                                      tCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCaaa                                                                                              >  1:202466/1‑149 (MQ=255)
                                                        aGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACg                                                                                            >  6:212191/1‑149 (MQ=255)
                                                           tGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCg                                                                                         <  5:9335/149‑1 (MQ=255)
                                                            gCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGt                                                                                        <  5:100791/149‑1 (MQ=255)
                                                             cGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGtt                                                                                       >  7:145639/1‑149 (MQ=255)
                                                             cGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGtt                                                                                       >  1:19143/1‑149 (MQ=255)
                                                                  ccTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCCGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGGCGACGATAAGCAAACGGCGTTGAGCg                                                                                  <  1:297989/149‑1 (MQ=255)
                                                                  ccTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCg                                                                                  <  1:293541/149‑1 (MQ=255)
                                                                         gcgcAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGgct                                                                           >  5:893/1‑149 (MQ=255)
                                                                              aTAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCtgctgc                                                                      <  8:155232/149‑1 (MQ=255)
                                                                                                tgtgATCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGCGCGGCAGGCTGCTGCATGACCGATCCTTGgcgc                                                    >  4:48658/1‑149 (MQ=255)
                                                                                                 gtgGTCCAGCGCAGCCCCGAGCATACAGCTGTCCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGgc                                                      >  1:203440/1‑147 (MQ=255)
                                                                                                 gtgGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGgc                                                      <  2:203440/146‑1 (MQ=255)
                                                                                                     tCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGtt                                               >  2:192383/1‑149 (MQ=255)
                                                                                                     tCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGtt                                               >  2:73330/1‑149 (MQ=255)
                                                                                                        aGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGttttt                                            <  3:79123/149‑1 (MQ=255)
                                                                                                         gcgcAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTg                                           >  8:262567/1‑149 (MQ=255)
                                                                                                         gcgcAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTg                                           >  1:119345/1‑149 (MQ=255)
                                                                                                                 cccGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttttt                                    <  1:192383/149‑1 (MQ=255)
                                                                                                                   cGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACCAGCGGTGGCAACCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGa                                  <  7:169478/149‑1 (MQ=255)
                                                                                                                    gAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 <  6:893/149‑1 (MQ=255)
                                                                                                                           cAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTaaaa                           <  8:193299/148‑1 (MQ=255)
                                                                                                                           cAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTaaaa                           <  5:188580/148‑1 (MQ=255)
                                                                                                                              cTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGc                       <  8:75486/149‑1 (MQ=255)
                                                                                                                                gTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 <  6:2347/137‑1 (MQ=255)
                                                                                                                                gTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 >  5:2347/1‑137 (MQ=255)
                                                                                                                                    cacaTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGtt                  >  4:313499/1‑149 (MQ=255)
                                                                                                                                    cacaTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGtt                  >  1:47838/1‑149 (MQ=255)
                                                                                                                                      caTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTtgt                <  2:202466/149‑1 (MQ=255)
                                                                                                                                       aTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTtgtg               >  6:324322/1‑149 (MQ=255)
                                                                                                                                       aTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTtgtg               <  5:268869/149‑1 (MQ=255)
                                                                                                                                         cctggCCTGCGCAACGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 >  7:275410/5‑129 (MQ=255)
                                                                                                                                         cctggCCTGCGCAACGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 <  8:275410/125‑1 (MQ=255)
                                                                                                                                          ctggCCTGCGCAACGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttt                                      >  7:113447/4‑123 (MQ=255)
                                                                                                                                          ctggCCTGCGCAACGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttt                                      <  8:113447/120‑1 (MQ=255)
                                                                                                                                            aGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGcc          <  7:262567/149‑1 (MQ=255)
                                                                                                                                              ccTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttttt                                    <  1:213482/121‑1 (MQ=255)
                                                                                                                                              ccTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttttt                                    >  2:213482/1‑121 (MQ=255)
                                                                                                                                              ccTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCg        <  5:324322/149‑1 (MQ=255)
                                                                                                                                                   gTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCGTTTaa   >  7:79836/1‑149 (MQ=255)
                                                                                                                                                    ttGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCGTTTAAc  <  1:73330/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CAATCCATTCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCGTTTAAC  >  NC_002947/3974244‑3974538

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: