Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,279 C→G 9.4% K569N (AAG→AAC)  PP_4872 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,539,2790CG 9.4% 255.9 / 14.2 85K569N (AAG→AAC) PP_4872hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/40);  new base G (3/5);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATT  >  NC_002947/5539131‑5539427
                                                                                                                                                    |                                                                                                                                                     
aCGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAc                                                                                                                                                       >  4:235332/1‑149 (MQ=255)
     tAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtctc                                                                                                                                                  >  6:197464/1‑149 (MQ=255)
       gAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCg                                                                                                                                                 >  3:252696/1‑148 (MQ=255)
        aaCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGaa                                                                                                                                               >  8:1369/1‑149 (MQ=255)
         aCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaa                                                                                                                                              <  7:306690/149‑2 (MQ=255)
          cGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaag                                                                                                                                             <  1:334575/149‑3 (MQ=255)
           gCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaag                                                                                                                                             <  1:173590/148‑3 (MQ=255)
             cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGTGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgc                                                                                                                                          <  6:63461/149‑1 (MQ=255)
              tGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgcg                                                                                                                                         >  7:61571/1‑149 (MQ=255)
                tACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGa                                                                                                                                                >  2:133076/1‑140 (MQ=255)
                tACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGa                                                                                                                                                <  1:133076/140‑1 (MQ=255)
                  ccTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCtt                                                                                                                                      <  4:112205/148‑1 (MQ=255)
                  ccTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCtt                                                                                                                                      >  4:60651/1‑148 (MQ=255)
                   cTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACAGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTg                                                                                                                                     >  6:9937/1‑148 (MQ=255)
                       tCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCg                                                                                                                                <  8:339238/149‑1 (MQ=255)
                        cGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCgg                                                                                                                               <  4:252696/149‑1 (MQ=255)
                         ggCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGt                                                                                                                              <  8:103919/149‑1 (MQ=255)
                           cgcgGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCgg                                                                                                                               >  3:144926/1‑146 (MQ=255)
                           cgcgGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCgg                                                                                                                               <  4:144926/146‑1 (MQ=255)
                            gcgGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGt                                                                                                                           <  2:142740/149‑1 (MQ=255)
                              ggAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAg                                                                                                                         >  6:339080/1‑149 (MQ=255)
                              ggAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAg                                                                                                                         >  5:19206/1‑149 (MQ=255)
                               gAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaag                                                                                                                                             >  8:172212/1‑126 (MQ=255)
                               gAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaag                                                                                                                                             <  7:172212/128‑3 (MQ=255)
                                    cGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGaa                                                                                                                                               <  5:324086/121‑1 (MQ=255)
                                    cGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGaa                                                                                                                                               >  6:324086/1‑121 (MQ=255)
                                        gCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGTTCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGCACTTCTCCAACCTGCCCTTGAAGCGCTCGTAGGTGCCGATAc                                                                                                               >  6:389569/1‑149 (MQ=255)
                                         cccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCg                                                                                                                                >  1:270901/1‑131 (MQ=255)
                                         cccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCg                                                                                                                                <  2:270901/131‑1 (MQ=255)
                                          ccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATAccc                                                                                                             >  7:229676/1‑149 (MQ=255)
                                           cTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCg                                                                                                            >  4:207390/1‑149 (MQ=255)
                                            ttCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCgg                                                                                                           >  2:47262/1‑149 (MQ=255)
                                              cAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCgg                                                                                                           <  4:272797/147‑1 (MQ=255)
                                                  cacaGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTtcat                                                                                                     <  4:209922/149‑1 (MQ=255)
                                                   acaGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTtcat                                                                                                     >  1:47284/1‑148 (MQ=255)
                                                    caGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTtcatca                                                                                                   <  3:235332/149‑1 (MQ=255)
                                                     aGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTcatcat                                                                                                  <  3:218320/149‑1 (MQ=255)
                                                       ggCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCt                                                                                                <  4:116801/149‑1 (MQ=255)
                                                        gCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCtt                                                                                               <  1:391065/149‑1 (MQ=255)
                                                           aGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTc                                                                                            <  5:339080/149‑1 (MQ=255)
                                                             gTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTc                                                                                            >  7:71441/1‑147 (MQ=255)
                                                             gTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTc                                                                                            >  6:71266/1‑147 (MQ=255)
                                                                 atcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑ggg                                                                                      <  3:178920/148‑1 (MQ=255)
                                                                    atGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGc                                                                                   >  3:206833/1‑148 (MQ=255)
                                                                     tGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCa                                                                                 <  8:71441/149‑1 (MQ=255)
                                                                      gTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCAc                                                                                >  5:368666/1‑149 (MQ=255)
                                                                        cTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCa                                                                                 >  5:378297/1‑146 (MQ=255)
                                                                        cTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCa                                                                                 <  6:378297/146‑1 (MQ=255)
                                                                           cgtcgGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCt                                                                             <  8:229676/147‑1 (MQ=255)
                                                                           cgtcgGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTc                                                                            <  8:28286/148‑1 (MQ=255)
                                                                             tcgGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCtt                                                                                                        <  6:204478/119‑1 (MQ=255)
                                                                                gtggtACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGtgtg                                                                      <  6:19206/149‑1 (MQ=255)
                                                                                     aCCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCg                                                                 <  3:60651/149‑1 (MQ=255)
                                                                                        ggggTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCgggg                                                              <  4:206833/149‑1 (MQ=255)
                                                                                         gggTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCgggg                                                              >  1:409340/1‑148 (MQ=255)
                                                                                         gggTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCgggg                                                              >  7:316868/1‑148 (MQ=255)
                                                                                          ggTGCGCGGCAATGAAGTACGGGATCAGGGACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGGCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCa                                                            <  6:245038/149‑1 (MQ=255)
                                                                                            tGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAAt                                                          >  1:211093/1‑149 (MQ=255)
                                                                                              cgcgGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGGCCACGGGGGCCACGCTCGGTGTGCTCCGGGGCAATc                                                         <  1:47262/149‑1 (MQ=255)
                                                                                                 ggCAATGAAGTACGGGATCAGGTACTGCTCCTTACCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCtt                                                       >  5:269462/1‑147 (MQ=255)
                                                                                                              gggATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAAc                                        <  8:316868/149‑1 (MQ=255)
                                                                                                              gggATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAAc                                        <  4:416964/149‑1 (MQ=255)
                                                                                                                gATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgt                                      >  8:321950/1‑149 (MQ=255)
                                                                                                                gATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgt                                      >  4:74052/1‑149 (MQ=255)
                                                                                                                 aTCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtg                                     <  5:71266/149‑1 (MQ=255)
                                                                                                                  tCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtgg                                    <  2:47284/149‑1 (MQ=255)
                                                                                                                    aGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtggtg                                  <  2:211093/149‑1 (MQ=255)
                                                                                                                    aGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtggtg                                  <  3:207390/149‑1 (MQ=255)
                                                                                                                     ggTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGGAGCCGCCAACGtggggg                                 <  5:389569/149‑1 (MQ=255)
                                                                                                                     ggTACTGCTCCTTGCCAGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACACGGCTTCAGCATCTTGTCAA‑GTGGGCCACGCGCGGTGAGCTCCTGGGAAATCTTCAGGTAGCCGCCAACGaggtgg                                 >  8:384946/1‑149 (MQ=255)
                                                                                                                            cTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGAc                            >  5:70714/1‑147 (MQ=255)
                                                                                                                             tCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAg                         <  8:61571/149‑1 (MQ=255)
                                                                                                                               cTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAg                         <  6:368666/147‑1 (MQ=255)
                                                                                                                               cTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACACTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGAGCAGtt                       >  3:77160/1‑149 (MQ=255)
                                                                                                                                    ccGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCtt                   >  5:126515/1‑148 (MQ=255)
                                                                                                                                     cGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGa                 <  2:409340/149‑1 (MQ=255)
                                                                                                                                      cgcatgttCGAGAAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGAc                >  5:162864/7‑149 (MQ=255)
                                                                                                                                           tttCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACt           >  5:433133/1‑149 (MQ=255)
                                                                                                                                              cGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCg        >  5:24218/1‑149 (MQ=255)
                                                                                                                                               gagaAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGa                 >  7:299991/1‑139 (MQ=255)
                                                                                                                                               gagaAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGa                 <  8:299991/139‑1 (MQ=255)
                                                                                                                                               gagaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCATTGACGTACTc          <  6:243918/146‑1 (MQ=255)
                                                                                                                                               gagaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCATTGACGTACTc          >  5:243918/1‑146 (MQ=255)
                                                                                                                                                agaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGc      >  2:77067/1‑149 (MQ=255)
                                                                                                                                                agaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGc      <  3:74052/149‑1 (MQ=255)
                                                                                                                                                 gaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGGAATCTTCGGGTAGCGGCCAACGTGGTGGGTGACCAGTTCCTTGAGGTACTCCGGc      >  6:66867/1‑148 (MQ=255)
                                                                                                                                                    cTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGAtt  >  6:442952/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATT  >  NC_002947/5539131‑5539427

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