Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,276,840 G→A 7.4% intergenic (‑144/‑200) 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,8400GA 7.4% 276.6 / 7.2 94intergenic (‑144/‑200)cioA/PP_4652cyanide insensitive ubiquinol oxidase subunit I/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (31/55);  new base A (3/4);  total (35/59)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

ACATGGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCTTGTTACA  >  NC_002947/5276693‑5276987
                                                                                                                                                     |                                                                                                                                                   
aCATGGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCACACAGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGa                                                                                                                                                    >  6:164705/1‑149 (MQ=255)
 cATGGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGAt                                                                                                                                                   <  4:261470/149‑1 (MQ=255)
  aTGGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGAtt                                                                                                                                                  <  6:264359/149‑1 (MQ=255)
    ggTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGa                                                                                                                                                >  2:217286/1‑149 (MQ=255)
    ggTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGa                                                                                                                                                >  8:154711/1‑149 (MQ=255)
      tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGAGTAGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATATGGAACGATTGAc                                                                                                                                               >  1:325582/1‑149 (MQ=255)
      tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAcc                                                                                                                                              <  2:68176/149‑1 (MQ=255)
      tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAcc                                                                                                                                              <  3:106717/149‑1 (MQ=255)
        cTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGTGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAccc                                                                                                                                             <  1:84202/149‑1 (MQ=255)
        cTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAccc                                                                                                                                             <  4:165307/148‑1 (MQ=255)
             tctTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATc                                                                                                                                        >  3:16932/1‑148 (MQ=255)
             tctTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATc                                                                                                                                        >  6:97435/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                                                                                                                                 >  1:352507/1‑149 (MQ=255)
                   gATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaatgaat                                                                                                                                 >  6:175268/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:289234/149‑1 (MQ=255)
                                                              cccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      <  4:196011/149‑1 (MQ=255)
                                                              cccACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      <  2:325582/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                                                             <  8:100587/148‑1 (MQ=255)
                                                                                         gAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCa                                                             <  4:387176/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                                     <  6:132063/130‑1 (MQ=38)
                                                                                                                                   caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     >  5:126279/2‑131 (MQ=38)
                                                                                                                                   caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     >  5:132063/2‑131 (MQ=38)
                                                                                                                                   caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     <  6:126279/130‑1 (MQ=38)
                                                                                                                                   aaTTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGt                                                                           >  5:381341/1‑93 (MQ=255)
                                                                                                                                   aaTTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGt                                                                           <  6:381341/93‑1 (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                              <  5:378447/125‑1 (MQ=255)
                                                                                                                                                ggAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGa                              >  6:378447/1‑125 (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)
                                                                                                                                                     |                                                                                                                                                   
ACATGGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGGTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCTTGTTACA  >  NC_002947/5276693‑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: