Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,981,166 G→A 9.3% A179V (GCC→GTC)  cobW ← adenosylcobalamin biosynthesis protein CobW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,981,1660GA 9.3% 178.7 / 5.5 64A179V (GCC→GTC) cobWadenosylcobalamin biosynthesis protein CobW
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (3/3);  total (31/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 = 1.00e+00

CACCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGA  >  NC_002947/3981024‑3981313
                                                                                                                                              |                                                                                                                                                   
caCCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCt                                                                                                                                               <  3:212928/149‑1 (MQ=255)
    tCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGt                                                                                                                                            >  5:273609/1‑148 (MQ=255)
         tGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCAGCAGGTCGGCGCTGACCAGCTGGCCAGCgc                                                                                                                                      <  6:143277/149‑2 (MQ=255)
           gCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGaa                                                                                                                                     >  6:171201/1‑148 (MQ=255)
               cTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCt                                                                                                                                <  3:225586/149‑1 (MQ=255)
                  gCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGt                                                                                                                             >  5:236733/1‑149 (MQ=255)
                         cttgCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAgggg                                                                                                                      <  8:275115/146‑1 (MQ=255)
                         ctgGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAgggg                                                                                                                      <  2:89164/149‑1 (MQ=255)
                         ctgGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAggg                                                                                                                       <  6:189820/148‑1 (MQ=255)
                              cTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGAt                                                                                                                  >  3:36322/1‑148 (MQ=255)
                                 gATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCa                                                                                                               <  5:60997/148‑1 (MQ=255)
                                  aTCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATg                                                                                                             <  5:171201/149‑1 (MQ=255)
                                     aCTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTc                                                                                                          <  3:26163/149‑1 (MQ=255)
                                      cTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTc                                                                                                          >  8:172122/1‑148 (MQ=255)
                                       tttCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGCTTCATGGTCCa                                                                                                        >  8:80681/1‑149 (MQ=255)
                                       tttCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATgg                                                                                                            <  4:154935/145‑1 (MQ=255)
                                       tttCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATgg                                                                                                            >  3:154935/1‑145 (MQ=255)
                                             cgccgGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAg                                                                                                  >  6:66682/1‑149 (MQ=255)
                                                   cggcAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGa                                                                                            <  5:143151/149‑1 (MQ=255)
                                                     gcAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGc                                                                                          >  6:315198/1‑149 (MQ=255)
                                                        gTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTg                                                                                       <  8:316497/149‑1 (MQ=255)
                                                         ttcttcGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTgc                                                                                      >  5:131329/1‑149 (MQ=255)
                                                               ggCCACCTCGGCACGCACCTTGTCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                 <  4:36322/148‑1 (MQ=255)
                                                               ggCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                 <  5:66682/148‑1 (MQ=255)
                                                                   aCCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATc                                                                            <  7:172122/149‑1 (MQ=255)
                                                                     cTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGa                                                                          >  2:105551/1‑149 (MQ=255)
                                                                        ggCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGAc                                                                         <  6:131329/147‑1 (MQ=255)
                                                                           acgcacCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCACCACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGt                                                                    >  5:182721/1‑149 (MQ=255)
                                                                            cgcacCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTc                                                                   <  1:105551/149‑1 (MQ=255)
                                                                             gcacCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTcc                                                                  >  7:157671/1‑149 (MQ=255)
                                                                                ccTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCggg                                                               <  6:273609/149‑1 (MQ=255)
                                                                                  ttGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                 <  6:352028/129‑1 (MQ=255)
                                                                                  ttGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                 >  5:352028/1‑129 (MQ=255)
                                                                                          gCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCa                                                      >  3:187340/1‑148 (MQ=255)
                                                                                             ttCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGTTCGGCGCTGGCCAGCAGGTCAGCGAACAGCTCGTGCAGGGGCCATTCATGGTCCAGCTTACGCTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAgg                                                  >  8:292017/1‑149 (MQ=255)
                                                                                               cGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGt                                                 <  8:157671/148‑1 (MQ=255)
                                                                                                 gCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAg                                                   <  5:120408/144‑1 (MQ=255)
                                                                                                 gCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAg                                                   >  6:120408/1‑144 (MQ=255)
                                                                                                      gATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCgg                                         >  5:137380/1‑149 (MQ=255)
                                                                                                       aTCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGc                                        >  5:58155/1‑149 (MQ=255)
                                                                                                        tCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGcc                                       >  7:45951/1‑149 (MQ=255)
                                                                                                           ggTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACg                                    >  2:11420/1‑149 (MQ=255)
                                                                                                               ggCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCg                                <  5:315198/149‑1 (MQ=255)
                                                                                                                gCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCgg                               <  6:58155/149‑1 (MQ=255)
                                                                                                                 ccTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCggg                              <  3:204473/149‑1 (MQ=255)
                                                                                                                     gTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCt                            <  7:80681/147‑1 (MQ=255)
                                                                                                                      ttCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCATGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                         <  6:236733/149‑1 (MQ=255)
                                                                                                                         aaCACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCa                                             >  5:316372/1‑126 (MQ=255)
                                                                                                                         aaCACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCa                                             <  6:316372/126‑1 (MQ=255)
                                                                                                                          acacCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGAc                      >  3:331261/1‑148 (MQ=255)
                                                                                                                            accaccAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGCCCTGGTCCGGGTAGGCGGCATAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCAc                   >  6:6254/1‑149 (MQ=255)
                                                                                                                             ccaccaGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACg                  >  2:230916/1‑149 (MQ=255)
                                                                                                                              caccaGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCt                            >  8:60167/1‑138 (MQ=255)
                                                                                                                              caccaGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCt                            <  7:60167/138‑1 (MQ=255)
                                                                                                                              caccaGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACgg                 >  2:169986/1‑149 (MQ=255)
                                                                                                                               accaGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGt                <  8:45951/149‑1 (MQ=255)
                                                                                                                                    gTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGa           <  4:331261/149‑1 (MQ=255)
                                                                                                                                    ctgttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGAGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGa           >  5:249853/6‑149 (MQ=255)
                                                                                                                                      gttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                         >  1:203889/4‑133 (MQ=255)
                                                                                                                                      gttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                         >  3:240735/4‑133 (MQ=255)
                                                                                                                                      gttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                         <  2:203889/130‑1 (MQ=255)
                                                                                                                                      gttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                         <  4:240735/130‑1 (MQ=255)
                                                                                                                                          gcTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCAt     >  3:349487/1‑149 (MQ=255)
                                                                                                                                          gcTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCAt     <  1:11420/149‑1 (MQ=255)
                                                                                                                                              gCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGa  <  6:137380/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CACCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGA  >  NC_002947/3981024‑3981313

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑