Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,515,926 G→C 5.7% P34A (CCC→GCC)  cbpA ← curved DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,515,9260GC 5.7% 225.4 / 3.5 70P34A (CCC→GCC) cbpAcurved DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (28/38);  new base C (2/2);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TGGTGCCTGGAACGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGATAGGGGT  >  NC_002947/5515780‑5516073
                                                                                                                                                     |                                                                                                                                                     
tggtgCCTGGAACGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGt                                                                                                                                                     <  3:41509/149‑1 (MQ=255)
tggtgCCTGGAACGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGt                                                                                                                                                     <  7:360475/149‑1 (MQ=255)
tggtgCCTGGAACGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGt                                                                                                                                                     <  4:171329/149‑1 (MQ=255)
  gtgCCTGGAAAGGCCGGGCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTGTTCGGGGTCCTCCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCggg                                                                                                                                                       <  1:330509/148‑1 (MQ=255)
     ccTGGAACGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATAc                                                                                                                                               >  4:358485/1‑149 (MQ=255)
         gAACGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACtt                                                                                                                                             <  4:14113/147‑1 (MQ=255)
           aCGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTgcgc                                                                                                                                         <  1:355667/149‑1 (MQ=255)
           aCGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTgcgc                                                                                                                                         <  5:23078/149‑1 (MQ=255)
               ccggccATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTgcgcgc                                                                                                                                       >  4:331039/1‑147 (MQ=255)
               ccggccATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTgcgcgc                                                                                                                                       >  5:207073/1‑147 (MQ=255)
                     aTGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGc                                                                                                                               <  6:207073/149‑1 (MQ=255)
                      tGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCg                                                                                                                              <  3:331039/149‑1 (MQ=255)
                      tGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCg                                                                                                                              >  3:68455/1‑149 (MQ=255)
                         tGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTa                                                                                                                           >  3:42615/1‑149 (MQ=255)
                           gccgccGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTAcg                                                                                                                         >  6:199196/1‑149 (MQ=255)
                              gccgTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCgg                                                                                                                      <  7:101370/149‑1 (MQ=255)
                               ccgTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGc                                                                                                                     >  3:354477/1‑149 (MQ=255)
                               ccgTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGc                                                                                                                     <  7:223774/149‑1 (MQ=255)
                                 gTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCt                                                                                                                   <  3:358485/149‑1 (MQ=255)
                                  tACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTCCGCGCCAGTTTGCGGTACGCGGCCtt                                                                                                                  <  1:12297/149‑1 (MQ=255)
                                     ttGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGAt                                                                                                               >  1:53195/1‑149 (MQ=255)
                                       gCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATc                                                                                                              <  8:222892/148‑1 (MQ=255)
                                         ggATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGcc                                                                                                           >  2:220469/1‑149 (MQ=255)
                                          gATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCt                                                                                                          >  1:300341/1‑149 (MQ=255)
                                           aTTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCtt                                                                                                         <  3:113029/149‑1 (MQ=255)
                                           aTTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCtt                                                                                                         <  3:149515/149‑1 (MQ=255)
                                              tCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCtttt                                                                                                       <  4:330508/148‑1 (MQ=255)
                                               cATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTc                                                                                                      <  2:229291/148‑1 (MQ=255)
                                               cATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCg                                                                                                     >  3:356503/1‑149 (MQ=255)
                                                aTCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGt                                                                                                    >  4:365340/1‑149 (MQ=255)
                                                 tCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTc                                                                                                   <  4:354477/149‑1 (MQ=255)
                                                   gAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGCATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTcc                                                                                                 <  1:95154/149‑1 (MQ=255)
                                                   gAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTcc                                                                                                 <  3:334022/148‑1 (MQ=255)
                                                   gAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCg                                                                                                <  1:220469/149‑1 (MQ=255)
                                                         ggCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTgg                                                                                          >  6:368651/1‑149 (MQ=255)
                                                          gCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTc                                                                                                      >  8:353810/1‑137 (MQ=255)
                                                          gCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTc                                                                                                      <  7:353810/137‑1 (MQ=255)
                                                               tttCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTa                                                                                    <  4:42615/149‑1 (MQ=255)
                                                                     cgcgTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGc                                                                                 <  8:32480/146‑1 (MQ=255)
                                                                     cgcgTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGa                                                                              <  5:199196/149‑1 (MQ=255)
                                                                         gTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTg                                                                                                                                            >  8:96370/1‑85 (MQ=255)
                                                                         gTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTg                                                                                                                                            <  7:96370/85‑1 (MQ=255)
                                                                           c‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATGGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTTGGCTCTACGCCGAGTATc                                                                         <  1:250049/149‑1 (MQ=255)
                                                                              gCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATc                                                                         <  2:300341/148‑1 (MQ=255)
                                                                              gCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCt                                                                        >  2:80208/1‑149 (MQ=255)
                                                                                ccAGCACCACAAAGGCCACGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGAGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTaacatg                                                                      >  8:309854/1‑144 (MQ=255)
                                                                                 cAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTg                                                                      >  2:304177/1‑148 (MQ=255)
                                                                                       cTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTg                                                                                                                                >  5:163734/1‑85 (MQ=255)
                                                                                       cTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTg                                                                                                                                <  6:163734/85‑1 (MQ=255)
                                                                                       cTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCCTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGt                                                               <  4:68455/149‑1 (MQ=255)
                                                                                             ggCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTg                                                          <  6:218183/148‑1 (MQ=255)
                                                                                              gCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTg                                                          >  2:346587/1‑147 (MQ=255)
                                                                                               ccTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAg                                                       >  7:223528/1‑149 (MQ=255)
                                                                                                       gCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTc                                                                                                      <  3:174364/95‑1 (MQ=255)
                                                                                                       gCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTc                                                                                                      >  4:174364/1‑95 (MQ=255)
                                                                                                         ctcctTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCTCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTAt                                              <  1:304177/149‑1 (MQ=255)
                                                                                                            ctTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATgg                                            <  8:286308/147‑1 (MQ=255)
                                                                                                            ctTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATgg                                            <  1:62810/147‑1 (MQ=255)
                                                                                                            ctTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGCGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGAt                                          >  6:186667/1‑149 (MQ=255)
                                                                                                            ctTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGCGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGAt                                          >  6:180655/1‑149 (MQ=255)
                                                                                                                         tcGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCaatgtgaat                             >  2:351625/1‑149 (MQ=255)
                                                                                                                           ggCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCaatgtgaatgt                           <  1:80208/149‑1 (MQ=255)
                                                                                                                           ggCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCaatgtgaat                             <  3:365340/147‑1 (MQ=255)
                                                                                                                           ggCGTCGCGCTCCTTACTGACATCGGCGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCaatgtgaat                             <  5:180655/147‑1 (MQ=255)
                                                                                                                           ggCGTCGCGCTCCTTACTGACATCGGCGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCaatgtgaat                             <  5:186667/147‑1 (MQ=255)
                                                                                                                                  cgcTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGtt                    >  7:267786/1‑149 (MQ=255)
                                                                                                                                    cTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGc                  <  4:356503/149‑1 (MQ=255)
                                                                                                                                     tCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGc                  >  1:63995/1‑148 (MQ=255)
                                                                                                                                           cTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGaa           >  7:218716/1‑149 (MQ=255)
                                                                                                                                              aCATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGa         >  6:279363/1‑148 (MQ=255)
                                                                                                                                                 tCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGATAgg     >  3:165161/1‑149 (MQ=255)
                                                                                                                                                   ggggTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGATAgggg   >  4:33351/1‑149 (MQ=255)
                                                                                                                                                   ggggTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGATAggg    <  1:351625/148‑1 (MQ=255)
                                                                                                                                                   ggggTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGATAggg    <  8:223528/148‑1 (MQ=255)
                                                                                                                                                   ggggTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGATAggg    <  2:19106/148‑1 (MQ=255)
                                                                                                                                                    gggTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGATAGGGGt  >  2:403836/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTGCCTGGAACGGCCGGCCATGCTGGCCGCCGTACTTGCGGATTTCATCGAACTCGGCACGTTTCTGCGCGTC‑‑GCCCAGCACCTCATAGGCCTCGTTGGCCTCCTTGAATTTCTCCTCGGCGTCGCGCTCCTTACTGACATCGGGGTGATACTTGCGCGCCAGTTTGCGGTACGCGGCCTTGATCGCCTTTTCGTCCGCAGTGGGCTCTACGCCGAGTATCTTGTAATAGTCTTTGAAGTCCATTTATGGATCACCAATGTGAATGTGCATGTTGCAACTGAAGATAGGGGT  >  NC_002947/5515780‑5516073

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: