Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,980,280 G→A 7.9% D117D (GAC→GAT) ‡ cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,980,2800GA 7.9% 218.1 / 4.8 75D117D (GAC→GAT) ‡cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base G (39/30);  new base A (3/3);  total (42/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.36e-01

CGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTC  >  NC_002947/3980141‑3980426
                                                                                                                                            |                                                                                                                                                   
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                                                                     >  2:374832/1‑141 (MQ=255)
                                                                               aaaGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGc                                                                     <  1:374832/141‑1 (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                   >  3:58884/1‑149 (MQ=255)
                                                                                                                         tCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCAcc                   <  8:30913/149‑1 (MQ=255)
                                                                                                                         tCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCAcc                   >  4:274965/1‑149 (MQ=255)
                                                                                                                           aGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTg                 >  7:236188/1‑149 (MQ=255)
                                                                                                                           aGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTg                 >  6:82388/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                                         <  2:272494/114‑1 (MQ=255)
                                                                                                                                      cgggCGATGACCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                         >  1:277962/1‑114 (MQ=255)
                                                                                                                                      cgggCGATGACCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                         >  1:272494/1‑114 (MQ=255)
                                                                                                                                      cgggCGATGACCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                         <  2:277962/114‑1 (MQ=255)
                                                                                                                                          cGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTc  <  8:219605/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTC  >  NC_002947/3980141‑3980426

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: