Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,501,401 T→C 6.4% E103G (GAA→GGA)  PP_1313 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,501,4010TC 6.4% 179.0 / 2.7 62E103G (GAA→GGA) PP_1313AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (33/25);  new base C (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

AGTTGCACCGCCGGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGGGAGGATGATCACATCGGCGT  >  NC_002947/1501259‑1501548
                                                                                                                                              |                                                                                                                                                   
tactGCGCCGGCGGTGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCgg                                                                                                                                                  >  7:181336/4‑146 (MQ=255)
tactGCGCCGGCGGTGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCgg                                                                                                                                                  <  8:181336/143‑1 (MQ=255)
 gTTGCACCGCCGGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCTGGCCGAATCCGgcctg                                                                                                                                              >  3:241627/1‑147 (MQ=255)
     cACCGCCGGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGgccagcca                                                                                                                                           <  6:104718/148‑1 (MQ=255)
      aCCGCCGGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGccagccag                                                                                                                                          >  6:157093/1‑148 (MQ=255)
      aCCGCCGGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGccagccag                                                                                                                                          >  2:148612/1‑148 (MQ=255)
       ccgccgGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGaa                                                                                                                                        <  8:255391/149‑1 (MQ=255)
        cgccgGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCTGGCCGAATCCGgcctgctcgatg                                                                                                                                       <  4:241627/149‑7 (MQ=255)
        cgccgGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGAATCCGgcctgctcgatg                                                                                                                                       <  4:175774/149‑7 (MQ=255)
                cAGTAGTTTTTGGCGGGATCGGGCATGTGCTTGTCCTGGTTCATGCGGATCTGGGGGAAGCGTTCGGCAAAGCTGTTGTAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCTGATTCGGCCAGCCAGAACTCACCACt                                                                                                                               <  6:53798/149‑1 (MQ=255)
                       tatTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCg                                                                                                                        <  1:148612/149‑1 (MQ=255)
                         ttGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGgc                                                                                                                      <  4:77127/149‑1 (MQ=255)
                               gCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAg                                                                                                                >  5:165212/1‑149 (MQ=255)
                                           tGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAACGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCCCCGCGCCCcg                                                                                                    >  8:145990/1‑149 (MQ=255)
                                                   cTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGgcc                                                                                                                                                <  1:316559/97‑1 (MQ=255)
                                                   cTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGgcc                                                                                                                                                >  2:316559/1‑97 (MQ=255)
                                                   cTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTg                                                                                              >  3:74494/1‑147 (MQ=255)
                                                   cTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTg                                                                                              <  4:74494/147‑1 (MQ=255)
                                                     ggTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACCCACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGt                                                                                          >  7:300987/1‑149 (MQ=255)
                                                      gTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACCCCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTggt                                                                                         >  3:243392/1‑147 (MQ=255)
                                                                 tCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGc                                                                                                                            >  3:164926/1‑103 (MQ=255)
                                                                 tCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGc                                                                                                                            <  4:164926/103‑1 (MQ=255)
                                                                 tCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCAc                                                                              >  8:25940/1‑149 (MQ=255)
                                                                    gCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACgg                                                                            >  4:84772/1‑148 (MQ=255)
                                                                           gCGTTCGGCAAAGCTGTTGAAGAAGCGCCGGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCt                                                                    >  4:255955/1‑149 (MQ=255)
                                                                             gTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCtgt                                                                  <  1:172204/149‑1 (MQ=255)
                                                                             gTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCtgt                                                                  >  2:172204/1‑149 (MQ=255)
                                                                             gTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCtgt                                                                  <  7:25940/149‑1 (MQ=255)
                                                                             gTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTAGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCtgt                                                                  >  3:275693/1‑149 (MQ=255)
                                                                                cGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTgg                                                                <  1:203900/148‑1 (MQ=255)
                                                                                 ggCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTggg                                                               >  2:324482/1‑148 (MQ=255)
                                                                                          gTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTg                                                       >  7:76268/1‑147 (MQ=255)
                                                                                           ttGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTgca                                                     <  1:324482/148‑1 (MQ=255)
                                                                                             gaagaaGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTgcagc                                                   >  6:29362/1‑148 (MQ=255)
                                                                                              aagaagCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCaa                                                 <  7:306222/149‑1 (MQ=255)
                                                                                               agaagCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAg                                                >  7:135500/1‑149 (MQ=255)
                                                                                               agaagCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAg                                                >  6:297778/1‑149 (MQ=255)
                                                                                                gaagCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGt                                               >  1:109812/1‑149 (MQ=255)
                                                                                                 aagCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGtt                                              >  3:122586/1‑149 (MQ=255)
                                                                                                 aagCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGtt                                              >  2:215779/1‑149 (MQ=255)
                                                                                                   gcgcCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGt                                            >  3:299818/1‑149 (MQ=255)
                                                                                                     gcCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTc                                           >  2:156922/1‑148 (MQ=255)
                                                                                                      ccAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTc                                                                                        >  7:76072/1‑102 (MQ=255)
                                                                                                      ccAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTc                                                                                        <  8:76072/102‑1 (MQ=255)
                                                                                                       cAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGa                                         >  2:221063/1‑148 (MQ=255)
                                                                                                         gTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAat                                      <  1:221063/149‑1 (MQ=255)
                                                                                                            ggtggtCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAatcat                                   <  1:156922/149‑1 (MQ=255)
                                                                                                                   gCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTa                            <  1:215779/149‑1 (MQ=255)
                                                                                                                   gCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTa                            <  5:157093/149‑1 (MQ=255)
                                                                                                                            ccATCGAGCAGGCCGGATTCGGCCAGCCAGAACCCACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGgga                    <  8:300987/148‑1 (MQ=255)
                                                                                                                             cATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGggag                   >  5:3048/1‑148 (MQ=255)
                                                                                                                             cATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGggag                   >  6:57959/1‑148 (MQ=255)
                                                                                                                               tCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGcc                        <  8:117106/141‑1 (MQ=255)
                                                                                                                               tCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGcc                        >  7:117106/1‑141 (MQ=255)
                                                                                                                               tCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGGGAGgat                <  4:275693/149‑1 (MQ=255)
                                                                                                                                   gctggCCGAATCCGGCCTGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGaaa                                       >  2:116283/4‑122 (MQ=255)
                                                                                                                                   gCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGGGAGGATGATc            <  4:299818/149‑1 (MQ=255)
                                                                                                                                       gCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTc                                           >  5:213628/1‑114 (MQ=255)
                                                                                                                                       gCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTc                                           <  6:213628/114‑1 (MQ=255)
                                                                                                                                        ccGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGGGAGGATGATCACATc       >  3:256786/1‑149 (MQ=255)
                                                                                                                                         cGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGGGAGGATGATCACATCg      >  2:233509/1‑149 (MQ=255)
                                                                                                                                            aTTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGGGAGGATGATCACATCGGcg   <  4:256786/149‑1 (MQ=255)
                                                                                                                                             ttCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGGGAGGATGATCACATCGGcgt  <  3:255955/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
AGTTGCACCGCCGGCGCAGTAGATATTGTCGGCATCGGTCAGGTGCTTGTCCTGGTTCAGGCGGATCTGCGGGAAGCGTTCGGCAAAGCTGTTGAAGAAGCGCCAGTAGGTGGTCGCTTCCTTGCCATCGAGCAGGCCGGATTCGGCCAGCCAGAACACACCACTGGCCTCGGCGCACAGCACCGCGCCACGCGCGTGCTGTTCACGCAGCCACGGCAGCACCTGTGGGTAACGCTGCAGCAAGTTGTCGAAATCATCCCAGTAGGCCGGGAGGATGATCACATCGGCGT  >  NC_002947/1501259‑1501548

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: ‑