Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,481 G→A 5.5% intergenic (+408/‑343) mmsA‑I → / → PP_16SD methylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947697,4810GA 5.5% 262.8 / 3.7 91intergenic (+408/‑343)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (47/38);  new base A (2/3);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-02

ACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGAC  >  NC_002947/697334‑697628
                                                                                                                                                   |                                                                                                                                                    
acCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAgcgc                                                                                                                                                     >  2:12376/1‑149 (MQ=255)
acCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAgcgc                                                                                                                                                     >  6:323326/1‑149 (MQ=255)
  ccTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCtt                                                                                                                                                   >  2:371317/1‑149 (MQ=255)
  ccTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCATTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCtt                                                                                                                                                   >  1:87094/1‑149 (MQ=255)
   ctctTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTc                                                                                                                                                  <  3:139984/149‑1 (MQ=255)
   ctctTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTc                                                                                                                                                  >  1:274565/1‑149 (MQ=255)
        gcaTCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAACCTCCTTGAAATTCACTGAGTTAGCCAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCaa                                                                                                                                               <  2:320652/145‑1 (MQ=255)
        gAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGt                                                                                                                                             >  5:114794/1‑149 (MQ=255)
             cccTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCaa                                                                                                                                               <  8:24521/142‑1 (MQ=255)
             cccTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCaa                                                                                                                                               >  7:24521/1‑142 (MQ=255)
              ccTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGt                                                                                                                                       >  2:52725/1‑148 (MQ=255)
               cTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTAt                                                                                                                                     >  4:354044/1‑149 (MQ=255)
               cTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTAt                                                                                                                                     >  7:311347/1‑149 (MQ=255)
                  ataATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                    <  1:52725/147‑4 (MQ=25)
                  ataATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                    <  6:196405/147‑4 (MQ=25)
                  aaaaTGCGCCCCAATTCCAGCGACAACGGAACTACAAACTCCTTTAAATTAAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                    <  7:281860/145‑4 (MQ=25)
                      tccgcCCCACATCCAGCGACACCGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGc                                                                                                                              <  4:110362/147‑1 (MQ=255)
                      tGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGc                                                                                                                                          <  3:342376/137‑1 (MQ=255)
                      tGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTCTAAGGCCTAAGATCAAATGATCGCCAGCCGATGAGATGATGGTTGACATCGCTTACAAGTGc                                                                                                                                          >  4:342376/1‑137 (MQ=25)
                         gcCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTc                                                                                                                           <  2:284502/149‑1 (MQ=255)
                           cccACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTccc                                                                                                                         >  4:190692/1‑149 (MQ=255)
                                ttCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTAc                                                                                                                    >  6:232346/1‑149 (MQ=255)
                                   cAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAGATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                  <  5:300671/148‑1 (MQ=255)
                                   cAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                  <  8:311347/148‑1 (MQ=255)
                                    aGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgaga                                                                                                                <  2:192218/149‑1 (MQ=255)
                                         cAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCg                                                                                                           >  4:332795/1‑149 (MQ=255)
                                          aaCGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAAATCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGAAGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCg                                                                                                           >  1:332817/1‑148 (MQ=255)
                                                 tgtctaaCTCATTGAATTTCAAGGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGGTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGgagcgagc                                                                                                   >  4:203707/6‑149 (MQ=18)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGc                                                                                                                              >  5:86154/1‑112 (MQ=255)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGc                                                                                                                              <  6:86154/112‑1 (MQ=255)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTAc                                                                                                                    >  4:142564/1‑122 (MQ=255)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                          >  5:245308/1‑148 (MQ=255)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                          >  8:301579/1‑148 (MQ=255)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                          >  1:150846/1‑148 (MQ=255)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                          >  3:69985/1‑148 (MQ=255)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTg                                                                                         <  3:175101/149‑1 (MQ=255)
                                                           cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGATTCCAAGTGCTGTATGATTCGCCTCCAGCTACGAGAGATCGCAGCTAGCCAAGTGttt                                                                                          >  8:256478/1‑148 (MQ=255)
                                                           cTTGAAATTAAAAGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                          >  4:177550/1‑148 (MQ=255)
                                                            ttGAAATTCAATGAGTTAGACAATTCTGAACGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTCCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTg                                                                                         >  8:80533/1‑148 (MQ=255)
                                                              gAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGa                                                                                        <  2:274565/147‑1 (MQ=255)
                                                                    cAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGTTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTaaa                                                                                 <  4:69985/148‑1 (MQ=255)
                                                                     aaTGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAAc                                                                                <  5:100709/148‑1 (MQ=255)
                                                                     aaTGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAAc                                                                                >  6:10281/1‑148 (MQ=255)
                                                                           ttAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTtc                                                                         <  1:150405/149‑1 (MQ=255)
                                                                              gACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTc                                                                                                                                <  8:10932/91‑1 (MQ=255)
                                                                              gACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTc                                                                                                                                >  7:10932/1‑91 (MQ=255)
                                                                              gACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCg                                                                      <  6:114794/149‑1 (MQ=255)
                                                                               aCAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCg                                                                      >  4:282188/1‑148 (MQ=255)
                                                                                 aaTTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGaa                                                                                       <  2:304816/129‑1 (MQ=255)
                                                                                 aaTTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCa                                                                    >  5:164992/1‑148 (MQ=255)
                                                                                 aaTTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCCCTACGAGAGATCGCAGCGAGCCAAGTGTTTGaa                                                                                       >  1:304816/1‑129 (MQ=255)
                                                                                   ttCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCaaa                                                                  >  7:76961/1‑148 (MQ=255)
                                                                                   ttCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCaaa                                                                  >  8:116582/1‑148 (MQ=255)
                                                                                          ggTGCTGTAAGGGCTACGATCAAATGATAGCCAGCGGTCGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCa                                                          <  7:114957/149‑1 (MQ=255)
                                                                                            tGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGTAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCaaa                                                        <  1:278915/149‑1 (MQ=255)
                                                                                              cTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCAGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCaaaa                                                       >  5:279782/1‑148 (MQ=255)
                                                                                               tGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCACCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAAt                                                      >  3:240085/1‑148 (MQ=255)
                                                                                                   aGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCTCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAAc                                                  <  1:12376/148‑1 (MQ=255)
                                                                                                   aGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACg                                                 >  2:248025/1‑149 (MQ=39)
                                                                                                    gggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGGGTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAAc                                                  >  1:293598/1‑148 (MQ=255)
                                                                                                    gggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGc                                                <  7:116582/149‑1 (MQ=39)
                                                                                                    gggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGc                                                >  4:287928/1‑149 (MQ=39)
                                                                                                      gCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCtt                                              >  6:253670/1‑149 (MQ=39)
                                                                                                          cGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGa                                            <  1:248025/147‑1 (MQ=39)
                                                                                                            aTCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAgg                                        >  3:283422/1‑149 (MQ=39)
                                                                                                                aaTGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGctct                                    >  3:126766/1‑149 (MQ=37)
                                                                                                                 aTGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGctct                                    >  8:275622/1‑148 (MQ=37)
                                                                                                                 aTGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTg                                   <  6:245308/149‑1 (MQ=37)
                                                                                                                  tGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCATCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGa                                  <  7:256478/149‑1 (MQ=37)
                                                                                                                  tGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGa                                  <  3:282188/149‑1 (MQ=37)
                                                                                                                  tGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGa                                  <  3:190692/149‑1 (MQ=37)
                                                                                                                    aTCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAg                                 <  7:80533/148‑1 (MQ=37)
                                                                                                                     tCGCCAGCGGTTGAGATGATGCTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCTCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAgg                                <  2:246067/148‑1 (MQ=25)
                                                                                                                        ccAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAg                            >  3:120173/1‑149 (MQ=33)
                                                                                                                        ccAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAg                            >  3:120219/1‑149 (MQ=33)
                                                                                                                        ccAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAg                            >  1:78964/1‑149 (MQ=33)
                                                                                                                         cAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATaa                                                    <  2:301915/124‑1 (MQ=33)
                                                                                                                         cAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATaa                                                    >  1:301915/1‑124 (MQ=33)
                                                                                                                         cAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCAGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGc                           <  6:279782/149‑1 (MQ=17)
                                                                                                                          aGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCg                          <  1:371317/149‑1 (MQ=32)
                                                                                                                            cGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTa                        <  3:177550/149‑1 (MQ=32)
                                                                                                                            cGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTa                        <  7:301579/149‑1 (MQ=32)
                                                                                                                                   gatgatGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATgcg                 >  6:208920/1‑149 (MQ=32)
                                                                                                                                    atgatgGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATgcgc                <  3:187582/149‑1 (MQ=32)
                                                                                                                                    atgatgGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATgcgc                <  3:354044/149‑1 (MQ=32)
                                                                                                                                            ataCAGCACTTGGAAGCGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGt        >  5:335190/3‑149 (MQ=255)
                                                                                                                                             gACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAACCTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt       >  2:235839/1‑149 (MQ=17)
                                                                                                                                              aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt       <  2:78964/148‑1 (MQ=32)
                                                                                                                                              aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt       <  5:232346/148‑1 (MQ=32)
                                                                                                                                               cAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt       >  6:362906/1‑147 (MQ=32)
                                                                                                                                                aGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCt                                               >  2:208951/1‑106 (MQ=32)
                                                                                                                                                aGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt       >  3:118982/1‑146 (MQ=32)
                                                                                                                                                aGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTATCTCGCAAAAAACTTCAAAATAAACGCt                                               <  1:208951/106‑1 (MQ=17)
                                                                                                                                                   gcTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGACAGCGTAGAATGCGCGCCTCGGTTGAGAc  >  2:359506/1‑148 (MQ=17)
                                                                                                                                                   |                                                                                                                                                    
ACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGAC  >  NC_002947/697334‑697628

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: