Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,276,845 A→T 7.2% intergenic (‑149/‑195) cioA ← / → PP_4652 cyanide insensitive ubiquinol oxidase subunit I/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,276,8450AT 7.2% 277.5 / 7.7 92intergenic (‑149/‑195)cioA/PP_4652cyanide insensitive ubiquinol oxidase subunit I/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/54);  new base T (3/4);  total (34/58)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCTTGTTACA  >  NC_002947/5276697‑5276987
                                                                                                                                                      |                                                                                                                                              
ggTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGa                                                                                                                                                >  8:154711/1‑149 (MQ=255)
ggTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGa                                                                                                                                                >  2:217286/1‑149 (MQ=255)
  tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGAGTAGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATATGGAACGATTGAc                                                                                                                                               >  1:325582/1‑149 (MQ=255)
  tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAcc                                                                                                                                              <  3:106717/149‑1 (MQ=255)
  tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAcc                                                                                                                                              <  2:68176/149‑1 (MQ=255)
    cTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGTGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAccc                                                                                                                                             <  1:84202/149‑1 (MQ=255)
    cTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAccc                                                                                                                                             <  4:165307/148‑1 (MQ=255)
         tctTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATc                                                                                                                                        >  6:97435/1‑148 (MQ=255)
         tctTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATc                                                                                                                                        >  3:16932/1‑148 (MQ=255)
             cAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaatg                                                                                                                                    >  8:326309/1‑148 (MQ=255)
             cAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaat                                                                                                                                     <  3:267517/147‑1 (MQ=255)
               gATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGGTCAATCGTTCCAGATCaatgaat                                                                                                                                 <  7:249087/149‑1 (MQ=39)
               gATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaatgaat                                                                                                                                 >  6:175268/1‑149 (MQ=255)
               gATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaatgaat                                                                                                                                 >  1:352507/1‑149 (MQ=255)
                 tAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGc                                                                                                                               <  2:22873/149‑1 (MQ=255)
                  aaTCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGTAACGATTGACCCAGATCAATGAATGCt                                                                                                                              <  4:224924/149‑1 (MQ=255)
                  aaTCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCt                                                                                                                              <  5:393789/149‑1 (MQ=255)
                         cGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTTCGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGc                                                                                                                       <  5:331542/149‑1 (MQ=255)
                          gTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCa                                                                                                                      <  2:267771/149‑1 (MQ=255)
                                  ttGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTccc                                                                                                               >  3:56726/1‑148 (MQ=255)
                                    gtggaaccAGTCGCCCACTGACCCCACTGGGTCGAGTACGGCGTGTTGGAAATGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACTATTGACCCACATCAATGAATGCTAGGAAGCATAGTCCCGaa                                                                                                            >  3:246251/7‑149 (MQ=255)
                                       gCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATg                                                                                                          <  6:151965/148‑1 (MQ=255)
                                         cccAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGcc                                                                                                       >  3:316694/1‑149 (MQ=255)
                                         cccAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGcc                                                                                                       >  7:100587/1‑149 (MQ=255)
                                         cccAGGCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGcc                                                                                                       <  6:110180/149‑1 (MQ=255)
                                           cAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCg                                                                                                     <  5:97435/149‑1 (MQ=255)
                                            aGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCgg                                                                                                    >  2:3722/1‑149 (MQ=255)
                                             gTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGt                                                                                                   <  1:217286/149‑1 (MQ=255)
                                              tCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTg                                                                                                  <  7:159127/149‑1 (MQ=255)
                                                   cACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGCCGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCg                                                                                             <  5:265548/149‑1 (MQ=255)
                                                     cTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCgtg                                                                                           <  8:317601/149‑1 (MQ=255)
                                                      tGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCgtgt                                                                                          <  3:71184/149‑1 (MQ=255)
                                                        cccccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtgg                                                                                        <  6:312890/149‑1 (MQ=255)
                                                        cccccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAAGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtgg                                                                                        >  8:84533/1‑149 (MQ=255)
                                                         ccccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggt                                                                                       <  5:177829/149‑1 (MQ=255)
                                                          cccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCGGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      >  5:81935/1‑149 (MQ=255)
                                                          cccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      <  4:196011/149‑1 (MQ=255)
                                                          cccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      <  2:325582/149‑1 (MQ=255)
                                                          cccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      <  4:289234/149‑1 (MQ=255)
                                                           ccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      >  8:306214/1‑148 (MQ=255)
                                                           ccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGtggtgg                                                                                     <  3:398738/149‑1 (MQ=255)
                                                                 ggTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGcc                                                                               <  2:352507/149‑1 (MQ=255)
                                                                 ggTCGAGTACGGCGTGTTGGAATTGTGCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGcc                                                                               <  7:154711/149‑1 (MQ=255)
                                                                    tCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCcgc                                                                             <  4:316694/149‑1 (MQ=255)
                                                                    tCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCcgc                                                                             >  7:331858/1‑149 (MQ=255)
                                                                      gattACGGCGTGTTGGAATTGTTCTGTTCGATCTCAGGGTTTACGGCCAGACGGCGAATTCATTGATCTGGAACGATTGACCCATATCAATGAATGCTAGGAAGCATAGTCCCGAATGGGCCGGTGGGCCGTGTGGTTTATTGCCTCGt                                                                           <  4:246251/146‑1 (MQ=255)
                                                                         tACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGAc                                                                        <  7:306214/149‑1 (MQ=255)
                                                                           cGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGAcc                                                                       <  4:56726/148‑1 (MQ=255)
                                                                             gCGTGTTGGAATTGATCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCGCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGt                                                                    <  4:112741/149‑1 (MQ=255)
                                                                             gCGTGTTGGAATTGATCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGt                                                                    >  3:112741/1‑149 (MQ=255)
                                                                              cGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCGGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGtt                                                                   <  6:81935/149‑1 (MQ=255)
                                                                                     gAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCa                                                             <  4:387176/148‑1 (MQ=255)
                                                                                     gAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCa                                                             <  8:100587/148‑1 (MQ=255)
                                                                                       aTTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCAGAAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCAccc                                                          >  1:333978/1‑149 (MQ=255)
                                                                                          gttctgttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCctct                                                       >  8:180702/1‑149 (MQ=255)
                                                                                           ttctgttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCctct                                                       <  4:16932/148‑1 (MQ=255)
                                                                                           ttctgttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCtt                                                      >  4:111988/1‑149 (MQ=255)
                                                                                             ctgttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCttt                                                     >  6:28851/1‑148 (MQ=255)
                                                                                             ctgttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCttt                                                     <  5:246156/148‑1 (MQ=255)
                                                                                             ctgttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCttt                                                     <  5:175268/148‑1 (MQ=255)
                                                                                               gttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGa                                                   <  8:331858/148‑1 (MQ=255)
                                                                                                ttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATTAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGAt                                                  <  5:164705/148‑1 (MQ=255)
                                                                                                ttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGAt                                                  <  3:111988/148‑1 (MQ=255)
                                                                                                ttcGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGAt                                                  >  8:233231/1‑148 (MQ=255)
                                                                                                       gCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTg                                          >  4:57602/1‑149 (MQ=255)
                                                                                                        cAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTg                                          >  4:126533/1‑148 (MQ=255)
                                                                                                           ggATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaa                                      <  7:326309/149‑1 (MQ=255)
                                                                                                                tCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTg                                  >  7:361453/1‑149 (MQ=255)
                                                                                                                  cGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTa                                                <  7:73495/133‑1 (MQ=255)
                                                                                                                  cGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTa                                                >  8:73495/1‑133 (MQ=255)
                                                                                                                         aCGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGa                         <  3:126533/149‑1 (MQ=255)
                                                                                                                         aCGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGa                         <  1:3722/149‑1 (MQ=255)
                                                                                                                              cAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCaaa                    <  7:233231/149‑1 (MQ=255)
                                                                                                                               caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     >  5:126279/2‑131 (MQ=38)
                                                                                                                               caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     <  6:132063/130‑1 (MQ=38)
                                                                                                                               caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     >  5:132063/2‑131 (MQ=38)
                                                                                                                               caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     <  6:126279/130‑1 (MQ=38)
                                                                                                                               aaTTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGt                                                                           <  6:381341/93‑1 (MQ=255)
                                                                                                                               aaTTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGt                                                                           >  5:381341/1‑93 (MQ=255)
                                                                                                                                aTTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACa                  <  7:180702/149‑1 (MQ=255)
                                                                                                                                    aTTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTg                          <  5:182885/137‑1 (MQ=39)
                                                                                                                                    aTTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTg                          >  6:182885/1‑137 (MQ=39)
                                                                                                                                    aTTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGAc                                                                        >  4:127886/1‑91 (MQ=255)
                                                                                                                                    aTTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGAc                                                                        <  3:127886/91‑1 (MQ=255)
                                                                                                                                         tCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTg                                          <  2:313286/116‑1 (MQ=255)
                                                                                                                                         tCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTg                                          >  1:313286/1‑116 (MQ=255)
                                                                                                                                            ggAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGa                              >  6:378447/1‑125 (MQ=255)
                                                                                                                                            ggAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGa                              <  5:378447/125‑1 (MQ=255)
                                                                                                                                             gAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCtt        <  2:333978/146‑1 (MQ=255)
                                                                                                                                             gAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCtt        <  8:361453/146‑1 (MQ=255)
                                                                                                                                             gAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCTTg       <  7:84533/147‑1 (MQ=255)
                                                                                                                                                 gATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCTTGTTACa  >  5:374085/1‑148 (MQ=255)
                                                                                                                                                  aTTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     <  8:254660/112‑1 (MQ=255)
                                                                                                                                                  aTTGACCCAGATCAATGAATGCGAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     >  7:254660/1‑112 (MQ=255)
                                                                                                                                                      aCCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTaa              >  8:199026/1‑131 (MQ=255)
                                                                                                                                                      aCCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTaa              <  7:199026/131‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
GGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCTTGTTACA  >  NC_002947/5276697‑5276987

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: