Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,527,993 A→G 5.6% S402G (AGC→GGC)  imuB → DNA repair DNA linked protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,527,9930AG 5.6% 212.3 / 2.1 72S402G (AGC→GGC) imuBDNA repair DNA linked protein
Reads supporting (aligned to +/- strand):  ref base A (36/31);  new base G (2/2);  total (39/33)
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.21e-01

CCTGGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGG  >  NC_002947/3527847‑3528138
                                                                                                                                                  |                                                                                                                                                 
ccTGGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAgc                                                                                                                                                 >  8:4965/1‑149 (MQ=255)
 cTGGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAgcc                                                                                                                                                >  3:260400/1‑149 (MQ=255)
 cTGGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAgcc                                                                                                                                                >  3:178130/1‑149 (MQ=255)
     gcgcAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGccc                                                                                                                                            >  6:161565/1‑149 (MQ=255)
      cgcAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCg                                                                                                                                           <  5:44668/149‑1 (MQ=255)
        cAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCGGCTGCCCggt                                                                                                                                         <  8:159249/149‑2 (MQ=255)
        cAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCggg                                                                                                                                         <  8:173694/149‑1 (MQ=255)
             gcgcGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTgg                                                                                                                                     <  6:311603/148‑1 (MQ=255)
              cgcgAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCt                                                                                                                                   >  6:143885/1‑149 (MQ=255)
              cgcgAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCt                                                                                                                                   >  1:238004/1‑149 (MQ=255)
                 gAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgtt                                                                                                                                <  7:4965/149‑1 (MQ=255)
                     ggTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCg                                                                                                                            >  6:267743/1‑149 (MQ=255)
                      gTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCgg                                                                                                                           >  6:155504/1‑149 (MQ=255)
                           gcgcCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCAcc                                                                                                                       >  5:127758/1‑148 (MQ=255)
                             gcccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGAt                                                                                                                    >  6:277965/1‑149 (MQ=255)
                              cccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATc                                                                                                                   <  5:277965/149‑1 (MQ=255)
                               ccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATcg                                                                                                                  <  4:295063/149‑1 (MQ=255)
                                  ccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCgcgc                                                                                                               >  7:103870/1‑149 (MQ=255)
                                   cggggggAGGAGGCCGGCAAGGCGCTGGGCGCAGCCGCCGATCCCCGCCCCGACTGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGAGCGCGCt                                                                                                              <  8:176676/147‑1 (MQ=255)
                                   cTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCt                                                                                                              <  6:225809/149‑1 (MQ=255)
                                      ggggATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTg                                                                                                             <  2:238004/147‑1 (MQ=255)
                                      ggggATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGa                                                                                                           <  6:26325/149‑1 (MQ=255)
                                         gATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGa                                                                                                        >  4:223579/1‑149 (MQ=255)
                                         gATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGa                                                                                                        <  3:223579/149‑1 (MQ=255)
                                            gAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGAcgc                                                                                                     >  4:11035/1‑149 (MQ=255)
                                                 cGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACgccgcc                                                                                                 >  5:44961/1‑148 (MQ=255)
                                                  gTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTa                                                                                               >  7:178270/1‑149 (MQ=255)
                                                     aaGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTAcc                                                                                             <  6:44961/148‑1 (MQ=255)
                                                      aGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAg                                                                                           <  1:183797/149‑1 (MQ=255)
                                                           cTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTg                                                                                        >  7:179190/1‑147 (MQ=255)
                                                            ttGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAg                                                                                     <  2:163894/149‑1 (MQ=255)
                                                               gcgcAGCCGGCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAggg                                                                                   <  2:76897/148‑1 (MQ=255)
                                                               gcgcAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgttgt                                                                                                                              >  5:138122/1‑105 (MQ=255)
                                                               gcgcAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgttgt                                                                                                                              <  6:138122/105‑1 (MQ=255)
                                                               gcgcAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAggg                                                                                   <  5:143885/148‑1 (MQ=255)
                                                               gcgcAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGc                                                                                  >  4:102570/1‑149 (MQ=255)
                                                               gcgcAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGc                                                                                  >  5:117479/1‑149 (MQ=255)
                                                                        ccgATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCCGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGaa                                                                         >  1:45136/1‑149 (MQ=255)
                                                                         cgATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAAc                                                                        >  5:55710/1‑149 (MQ=255)
                                                                          gATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACg                                                                       >  7:72682/1‑149 (MQ=255)
                                                                                  ccccGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAAt                                                               <  5:155504/149‑1 (MQ=255)
                                                                                            gcCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtggtg                                                     <  6:320192/149‑1 (MQ=255)
                                                                                              cTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtggtgg                                                    >  1:171053/1‑148 (MQ=255)
                                                                                              cTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCATACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTggg                                                   >  6:5070/1‑149 (MQ=255)
                                                                                                 gCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGa                                                  <  4:178130/147‑1 (MQ=255)
                                                                                                    tttGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATggcg                                             <  3:11035/149‑1 (MQ=255)
                                                                                                       ggCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATggcggc                                           <  8:72682/148‑1 (MQ=255)
                                                                                                          ggAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCGCCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAtg                                       >  2:318347/1‑149 (MQ=255)
                                                                                                           gAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCTCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAtgt                                      <  4:260400/149‑1 (MQ=255)
                                                                                                             gCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCTGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGc                                    >  6:142292/1‑149 (MQ=255)
                                                                                                                gggTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTc                                 <  5:267743/149‑1 (MQ=255)
                                                                                                                  gTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCGCCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTcgc                               <  1:318347/149‑1 (MQ=255)
                                                                                                                      ccAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGac                            >  3:250572/1‑148 (MQ=255)
                                                                                                                      ccAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGac                            >  3:231572/1‑148 (MQ=255)
                                                                                                                        ggtgCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGactac                         >  5:104840/4‑149 (MQ=255)
                                                                                                                        aaGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGactac                         >  3:222410/1‑149 (MQ=255)
                                                                                                                          ggCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGActacta                       >  2:9289/1‑149 (MQ=255)
                                                                                                                            cAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTAcc                     <  6:127758/149‑1 (MQ=255)
                                                                                                                                aTGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCAt                 <  8:103870/149‑1 (MQ=255)
                                                                                                                                  gCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAtgt                                      <  5:35594/126‑1 (MQ=255)
                                                                                                                                  gCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAtgt                                      >  6:35594/1‑126 (MQ=255)
                                                                                                                                   ccGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTTCCGCAtag               >  6:300337/1‑146 (MQ=255)
                                                                                                                                    cGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGCCTACTACCGCATCGa              >  1:203036/1‑148 (MQ=255)
                                                                                                                                    cGGTTGCACCGGGCAGCCGCCCGGGGAGGCTGTTGTCGGCACCGATCGCGCTGGAAGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGATGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGaa             <  5:300337/149‑1 (MQ=255)
                                                                                                                                       ccaccccGGGCGGCTGCCCGGTGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAAcc          >  5:285487/6‑149 (MQ=255)
                                                                                                                                       ccaccccGGGCGGCTGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGAAGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAAGTGCGTCGCGACTACTACCGCATCGAAAcc          >  4:73273/6‑149 (MQ=255)
                                                                                                                                        caccccGGGCGGCTGCCCGGTGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAAccc         <  6:285487/145‑1 (MQ=255)
                                                                                                                                          cACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGt       <  2:171053/149‑1 (MQ=255)
                                                                                                                                           aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGt       >  1:206915/1‑148 (MQ=255)
                                                                                                                                           aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGt       >  6:207078/1‑148 (MQ=255)
                                                                                                                                           aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTCGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGATGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCAACGAACACCGt       >  8:252529/1‑148 (MQ=255)
                                                                                                                                               ggCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGG  <  4:222410‑M1/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCTGGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGG  >  NC_002947/3527847‑3528138

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