Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,276,843 T→C 7.4% intergenic (‑147/‑197) 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,8430TC 7.4% 275.5 / 8.1 94intergenic (‑147/‑197)cioA/PP_4652cyanide insensitive ubiquinol oxidase subunit I/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/54);  new base C (3/4);  total (36/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCTTGTTACA  >  NC_002947/5276695‑5276987
                                                                                                                                                      |                                                                                                                                                
aTGGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGAtt                                                                                                                                                  <  6:264359/149‑1 (MQ=255)
  ggTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGa                                                                                                                                                >  2:217286/1‑149 (MQ=255)
  ggTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGa                                                                                                                                                >  8:154711/1‑149 (MQ=255)
    tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGGAGTAGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATATGGAACGATTGAc                                                                                                                                               >  1:325582/1‑149 (MQ=255)
    tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAcc                                                                                                                                              <  3:106717/149‑1 (MQ=255)
    tGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAcc                                                                                                                                              <  2:68176/149‑1 (MQ=255)
      cTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGGTGATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAccc                                                                                                                                             <  1:84202/149‑1 (MQ=255)
      cTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGAccc                                                                                                                                             <  4:165307/148‑1 (MQ=255)
           tctTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATc                                                                                                                                        >  6:97435/1‑148 (MQ=255)
           tctTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATc                                                                                                                                        >  3:16932/1‑148 (MQ=255)
               cAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaatg                                                                                                                                    >  8:326309/1‑148 (MQ=255)
               cAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaat                                                                                                                                     <  3:267517/147‑1 (MQ=255)
                 gATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGGTCAATCGTTCCAGATCaatgaat                                                                                                                                 <  7:249087/149‑1 (MQ=39)
                 gATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaatgaat                                                                                                                                 >  6:175268/1‑149 (MQ=255)
                 gATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCaatgaat                                                                                                                                 >  1:352507/1‑149 (MQ=255)
                   tAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGc                                                                                                                               <  2:22873/149‑1 (MQ=255)
                    aaTCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGTAACGATTGACCCAGATCAATGAATGCt                                                                                                                              <  4:224924/149‑1 (MQ=255)
                    aaTCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCt                                                                                                                              <  5:393789/149‑1 (MQ=255)
                           cGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTTCGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGc                                                                                                                       <  5:331542/149‑1 (MQ=255)
                            gTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCa                                                                                                                      <  2:267771/149‑1 (MQ=255)
                                    ttGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTccc                                                                                                               >  3:56726/1‑148 (MQ=255)
                                      gtggaaccAGTCGCCCACTGACCCCACTGG‑GTCGAGTACGGCGTGTTGGAAATGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACTATTGACCCACATCAATGAATGCTAGGAAGCATAGTCCCGaa                                                                                                            >  3:246251/7‑149 (MQ=255)
                                         gCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATg                                                                                                          <  6:151965/148‑1 (MQ=255)
                                           cccAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGcc                                                                                                       >  3:316694/1‑149 (MQ=255)
                                           cccAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGcc                                                                                                       >  7:100587/1‑149 (MQ=255)
                                           cccAGGCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGcc                                                                                                       <  6:110180/149‑1 (MQ=255)
                                             cAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCg                                                                                                     <  5:97435/149‑1 (MQ=255)
                                              aGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCgg                                                                                                    >  2:3722/1‑149 (MQ=255)
                                               gTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGt                                                                                                   <  1:217286/149‑1 (MQ=255)
                                                tCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTg                                                                                                  <  7:159127/149‑1 (MQ=255)
                                                     cACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGCCGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCg                                                                                             <  5:265548/149‑1 (MQ=255)
                                                       cTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCgtg                                                                                           <  8:317601/149‑1 (MQ=255)
                                                        tGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCgtgt                                                                                          <  3:71184/149‑1 (MQ=255)
                                                          cccccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtgg                                                                                        <  6:312890/149‑1 (MQ=255)
                                                          cccccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAAG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtgg                                                                                        >  8:84533/1‑149 (MQ=255)
                                                           ccccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggt                                                                                       <  5:177829/149‑1 (MQ=255)
                                                            cccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCGGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      >  5:81935/1‑149 (MQ=255)
                                                            cccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      <  4:289234/149‑1 (MQ=255)
                                                            cccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      <  4:196011/149‑1 (MQ=255)
                                                            cccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      <  2:325582/149‑1 (MQ=255)
                                                             ccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTgtggtg                                                                                      >  8:306214/1‑148 (MQ=255)
                                                             ccACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGtggtgg                                                                                     <  3:398738/149‑1 (MQ=255)
                                                                   g‑gTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGcc                                                                               <  2:352507/149‑1 (MQ=255)
                                                                   g‑gTCGAGTACGGCGTGTTGGAATTGTGCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGcc                                                                               <  7:154711/149‑1 (MQ=255)
                                                                      tCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCcgc                                                                             <  4:316694/149‑1 (MQ=255)
                                                                      tCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCcgc                                                                             >  7:331858/1‑149 (MQ=255)
                                                                        gattACGGCGTGTTGGAATTGTTCTGTTCGATCTCAGG‑GTTTACGGCCAGACGGCGAATTCATTGATCTGGAACGATTGACCCATATCAATGAATGCTAGGAAGCATAGTCCCGAATGGGCCGGTGGGCCGTGTGGTTTATTGCCTCGt                                                                           <  4:246251/146‑1 (MQ=255)
                                                                           tACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGAc                                                                        <  7:306214/149‑1 (MQ=255)
                                                                             cGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGAcc                                                                       <  4:56726/148‑1 (MQ=255)
                                                                               gCGTGTTGGAATTGATCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCGCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGt                                                                    <  4:112741/149‑1 (MQ=255)
                                                                               gCGTGTTGGAATTGATCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGt                                                                    >  3:112741/1‑149 (MQ=255)
                                                                                cGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCGGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGtt                                                                   <  6:81935/149‑1 (MQ=255)
                                                                                       gAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCa                                                             <  8:100587/148‑1 (MQ=255)
                                                                                       gAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCa                                                             <  4:387176/148‑1 (MQ=255)
                                                                                         aTTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCAGAAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCAccc                                                          >  1:333978/1‑149 (MQ=255)
                                                                                            gttctgttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCctct                                                       >  8:180702/1‑149 (MQ=255)
                                                                                             ttctgttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCctct                                                       <  4:16932/148‑1 (MQ=255)
                                                                                             ttctgttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCtt                                                      >  4:111988/1‑149 (MQ=255)
                                                                                               ctgttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCttt                                                     <  5:246156/148‑1 (MQ=255)
                                                                                               ctgttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCttt                                                     <  5:175268/148‑1 (MQ=255)
                                                                                               ctgttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCttt                                                     >  6:28851/1‑148 (MQ=255)
                                                                                                 gttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGa                                                   <  8:331858/148‑1 (MQ=255)
                                                                                                  ttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATTAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGAt                                                  <  5:164705/148‑1 (MQ=255)
                                                                                                  ttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGAt                                                  >  8:233231/1‑148 (MQ=255)
                                                                                                  ttcGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGAt                                                  <  3:111988/148‑1 (MQ=255)
                                                                                                         gCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTg                                          >  4:57602/1‑149 (MQ=255)
                                                                                                          cAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTg                                          >  4:126533/1‑148 (MQ=255)
                                                                                                             g‑gATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaa                                      <  7:326309/149‑1 (MQ=255)
                                                                                                                  tCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTg                                  >  7:361453/1‑149 (MQ=255)
                                                                                                                    cGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTa                                                >  8:73495/1‑133 (MQ=255)
                                                                                                                    cGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTa                                                <  7:73495/133‑1 (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                                     <  6:126279/130‑1 (MQ=38)
                                                                                                                                 caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     >  5:132063/2‑131 (MQ=38)
                                                                                                                                 caTTCATTGATCTGGGTCAATCGTTCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGaaaa                                     >  5:126279/2‑131 (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)
                                                                                                                                                      |                                                                                                                                                
ATGGTGCTTCCTCTTCAGATAATCCGGCGTCCGGGCTTGTGGCCCCAGTCGCCCACTGCCCCCACTGG‑GTCGAGTACGGCGTGTTGGAATTGTTCTGTTCGATCGCAGG‑GATTCCGGCCAGACGGCCAATTCATTGATCTGGAACGATTGACCCAGATCAATGAATGCTAGGAAGCATAGTCCCGAATGGCCCGGTGGGCCGTGTGGTGGATTGCCGCGTGACCTGTTGCCTCACCCTCTTTGATTAGCAGTGGAAAAGTGCCGACCTGAGCAAACAGTAACCGCTTGTTACA  >  NC_002947/5276695‑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: ‑