Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,980,284 T→C 7.8% D116G (GAT→GGT)  cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,980,2840TC 7.8% 215.6 / 7.3 78D116G (GAT→GGT) cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base T (41/30);  new base C (3/3);  total (45/33)
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.67e-01

CGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCG  >  NC_002947/3980141‑3980431
                                                                                                                                                |                                                                                                                                                    
cgcgGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGccc                                                                                                                                                 >  8:242092/1‑149 (MQ=255)
cgcgGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCGCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGccc                                                                                                                                                 >  7:301136/1‑149 (MQ=255)
cgcgGCCCAGCCATTGGCTGGCCAGGCAGTGGAACCGGTTGAGGCCGCTGGCCGTGCCGCCCTGGCGCAGGGAATGCCTAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGGCAGCTCCGCGTCCGGGCGGTCATCGccc                                                                                                                                                 <  6:69801/149‑1 (MQ=255)
 gcgGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGccc                                                                                                                                                 <  5:151990/148‑1 (MQ=255)
  ctggccAGCCAATGGCTGGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCg                                                                                                                                                >  8:95020/5‑149 (MQ=255)
   ggCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCTCCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGc                                                                                                                                              <  6:330991/149‑1 (MQ=255)
    gCCCAGCCATTGGCCGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCa                                                                                                                                             <  2:366073/149‑1 (MQ=255)
      ccAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCAc                                                                                                                                            >  5:233692/1‑148 (MQ=255)
       cAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCAc                                                                                                                                            >  5:243574/1‑147 (MQ=255)
           cATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAAt                                                                                                                                     >  3:357418/1‑149 (MQ=255)
               ggctggcCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCCCCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGATGACCGCCCgga                                                                                                                                              >  8:368728/1‑136 (MQ=255)
               ggctggcCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGATGACCGCCCgga                                                                                                                                              <  7:368728/137‑2 (MQ=255)
                      ccAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACg                                                                                                                           >  7:218454/1‑149 (MQ=255)
                       cAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGc                                                                                                                          >  3:16118/1‑149 (MQ=255)
                        aGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTATTGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCAcc                                                                                                                                           >  2:332989/1‑132 (MQ=255)
                         ggCAGTTAAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCa                                                                                                                        >  4:90609/1‑149 (MQ=255)
                          gCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCAc                                                                                                                       <  6:233692/149‑1 (MQ=255)
                            aGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACgg                                                                                                                     <  7:217305/149‑1 (MQ=255)
                                      gTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCACAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCt                                                                                                           <  6:390728/149‑1 (MQ=255)
                                      gTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCCGCt                                                                                                           >  3:46825/1‑149 (MQ=255)
                                       ttGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTc                                                                                                          <  2:269444/149‑1 (MQ=255)
                                       ttGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTc                                                                                                          <  3:202406/149‑1 (MQ=255)
                                       ttGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTc                                                                                                          >  3:43006/1‑149 (MQ=255)
                                         gATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCg                                                                                                         <  3:60581/148‑1 (MQ=255)
                                          aTGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGAc                                                                                                       <  1:106083/149‑1 (MQ=255)
                                             gCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCa                                                                                                     <  8:218454/148‑1 (MQ=255)
                                               gTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAgtt                                                                                                  <  3:90609/149‑1 (MQ=255)
                                                  ggCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAgttgt                                                                                                <  5:227641/148‑1 (MQ=255)
                                                      ttGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGAc                                                                                           >  4:112676/1‑149 (MQ=255)
                                                      ttGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGAc                                                                                           >  2:102050/1‑149 (MQ=255)
                                                        gccgccCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGAc                                                                                           <  8:301136/147‑1 (MQ=255)
                                                        gccgccCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACg                                                                                          <  7:95020/148‑1 (MQ=255)
                                                             ccTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAg                                                                                    <  8:287646/149‑1 (MQ=255)
                                                                ggCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAgcgc                                                                                 <  6:243574/149‑1 (MQ=255)
                                                                     aGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTa                                                                            <  7:242092/149‑1 (MQ=255)
                                                                      ggTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAg                                                                           <  4:43006/149‑1 (MQ=255)
                                                                       gTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGc                                                                          >  5:87131/1‑149 (MQ=255)
                                                                          aTGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCaa                                                                       >  3:127399/1‑149 (MQ=255)
                                                                             ccAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGc                                                                    >  7:219605/1‑149 (MQ=255)
                                                                               aaaGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGc                                                                          <  1:374832/141‑1 (MQ=255)
                                                                               aaaGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGc                                                                          >  2:374832/1‑141 (MQ=255)
                                                                                                    cgcgGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCGGTAGCCAACGCCACCGTGCAGCGACACCAGGATg                                             >  2:332824/1‑149 (MQ=255)
                                                                                                    cgcgGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                              >  7:30913/1‑148 (MQ=255)
                                                                                                     gcgGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGa                                            >  3:5816/1‑149 (MQ=255)
                                                                                                         aCGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCtt                                        <  2:177827/149‑1 (MQ=255)
                                                                                                               cccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgt                                   <  4:357418/148‑1 (MQ=255)
                                                                                                                ccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCTGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                                 >  2:296617/1‑149 (MQ=255)
                                                                                                                ccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                                 >  5:59302/1‑149 (MQ=255)
                                                                                                                ccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                                 >  6:353381/1‑149 (MQ=255)
                                                                                                                ccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                                 >  7:254683/1‑149 (MQ=255)
                                                                                                                  aGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtggcg                               >  1:270290/1‑149 (MQ=255)
                                                                                                                  aGGCTGGTCAGCTCCGGGTCCGAGAGGTCAACGCCCGGAACCAGAATCAGCTGCACGCAACGGACGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACATtggcgtggcg                               >  7:59031/1‑149 (MQ=255)
                                                                                                                     cTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCgcag                            <  6:237097/149‑1 (MQ=255)
                                                                                                                         tCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCAcc                        <  8:30913/149‑1 (MQ=255)
                                                                                                                         tCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCAcc                        >  4:274965/1‑149 (MQ=255)
                                                                                                                         tCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCAcc                        >  3:58884/1‑149 (MQ=255)
                                                                                                                           aGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTg                      >  6:82388/1‑149 (MQ=255)
                                                                                                                           aGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTg                      >  7:236188/1‑149 (MQ=255)
                                                                                                                            gCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTgg                     <  4:16118/149‑1 (MQ=255)
                                                                                                                             ctctGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGt                    >  3:111377/1‑149 (MQ=255)
                                                                                                                             cTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGt                    >  7:282758/1‑149 (MQ=255)
                                                                                                                              tccgggtccgggCGGTCATCGCCCGGCACCGAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGCCCAGTTGTTCGACGCCGTAGCGCCAGTAGCCACCGCCCCCGTGCAGCGCCACCAGGATGACCTTGGCGTGGCGCAGCCCCTGGt                    >  1:235288/1‑149 (MQ=255)
                                                                                                                              tccgggtccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTAGGCGTGGCGCAGCACCTGCTc                   >  8:225527/1‑149 (MQ=255)
                                                                                                                              tccgggtccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGAGCAGCGACACCAGGATGCCCTTGGCGTGGCGCAGCACCTGGTc                   >  8:365462/1‑149 (MQ=255)
                                                                                                                                cgggtccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTc                   <  4:127399/147‑1 (MQ=255)
                                                                                                                                cgggtccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGa                 >  2:354747/1‑149 (MQ=255)
                                                                                                                                  ggtccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACa               >  2:5599/1‑149 (MQ=255)
                                                                                                                                     ccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATAc            <  4:111377/149‑1 (MQ=255)
                                                                                                                                      cgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATa             >  7:352462/1‑147 (MQ=255)
                                                                                                                                      cgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACa           <  8:236188/149‑1 (MQ=255)
                                                                                                                                      cgggCGATGACCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                              >  1:277962/1‑114 (MQ=255)
                                                                                                                                      cgggCGATGACCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                              <  2:272494/114‑1 (MQ=255)
                                                                                                                                      cgggCGATGACCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                              >  1:272494/1‑114 (MQ=255)
                                                                                                                                      cgggCGATGACCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                              <  2:277962/114‑1 (MQ=255)
                                                                                                                                          cGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTc       <  8:219605/149‑1 (MQ=255)
                                                                                                                                             tCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGAc    >  7:158236/1‑149 (MQ=255)
                                                                                                                                               aTCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGAcc   >  1:326128/1‑148 (MQ=255)
                                                                                                                                               aTCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCg  >  7:225007/1‑149 (MQ=255)
                                                                                                                                                tCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGAcc   >  7:13480/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCG  >  NC_002947/3980141‑3980431

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: