Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,981,175 T→C 9.5% D176G (GAC→GGC)  cobW ← adenosylcobalamin biosynthesis protein CobW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,981,1750TC 9.5% 168.8 / 7.4 63D176G (GAC→GGC) cobWadenosylcobalamin biosynthesis protein CobW
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (3/3);  total (30/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 = 9.79e-01

TCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGACGGTGCAGGC  >  NC_002947/3981028‑3981323
                                                                                                                                                   |                                                                                                                                                    
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                                                                                           >  5:352028/1‑129 (MQ=255)
                                                                              ttGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                           <  6:352028/129‑1 (MQ=255)
                                                                                      gCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCa                                                                >  3:187340/1‑148 (MQ=255)
                                                                                         ttCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGTTCGGCGCTGGCCAGCAGGTCAGCGAACAGCTCGTGCAGGGGCCATTCATGGTCCAGCTTACGCTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAgg                                                            >  8:292017/1‑149 (MQ=255)
                                                                                           cGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGt                                                           <  8:157671/148‑1 (MQ=255)
                                                                                             gCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAg                                                             >  6:120408/1‑144 (MQ=255)
                                                                                             gCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAg                                                             <  5:120408/144‑1 (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                                   <  4:240735/130‑1 (MQ=255)
                                                                                                                                  gttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                                   >  3:240735/4‑133 (MQ=255)
                                                                                                                                  gttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                                   <  2:203889/130‑1 (MQ=255)
                                                                                                                                  gttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                                   >  1:203889/4‑133 (MQ=255)
                                                                                                                                      gcTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCAt               <  1:11420/149‑1 (MQ=255)
                                                                                                                                      gcTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCAt               >  3:349487/1‑149 (MQ=255)
                                                                                                                                          gCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGa            <  6:137380/148‑1 (MQ=255)
                                                                                                                                                   tCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGACGGTGCAGGc  >  5:4709/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGACGGTGCAGGC  >  NC_002947/3981028‑3981323

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: ‑