Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,490 T→C 5.8% intergenic (+417/‑334) 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,4900TC 5.8% 251.6 / 4.3 86intergenic (+417/‑334)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (43/38);  new base C (2/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

GAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAG  >  NC_002947/697342‑697634
                                                                                                                                                    |                                                                                                                                                 
gAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGt                                                                                                                                                   >  5:114794/1‑149 (MQ=255)
      ccTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGt                                                                                                                                             >  2:52725/1‑148 (MQ=255)
       cTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTAt                                                                                                                                           >  7:311347/1‑149 (MQ=255)
       cTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTAt                                                                                                                                           >  4:354044/1‑149 (MQ=255)
          ataATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                          <  6:196405/147‑4 (MQ=25)
          ataATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                          <  1:52725/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                                                                                                >  3:69985/1‑148 (MQ=255)
                                                   cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                                >  5:245308/1‑148 (MQ=255)
                                                   cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                                >  8:301579/1‑148 (MQ=255)
                                                   cTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                                >  1:150846/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                                                                                      >  6:10281/1‑148 (MQ=255)
                                                             aaTGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAAc                                                                                      <  5:100709/148‑1 (MQ=255)
                                                                   ttAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTtc                                                                               <  1:150405/149‑1 (MQ=255)
                                                                      gACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTc                                                                                                                                      >  7:10932/1‑91 (MQ=255)
                                                                      gACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTc                                                                                                                                      <  8:10932/91‑1 (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                                                                        >  8:116582/1‑148 (MQ=255)
                                                                           ttCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCaaa                                                                        >  7:76961/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                                  >  1:78964/1‑149 (MQ=33)
                                                                                                                ccAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAg                                  >  3:120173/1‑149 (MQ=33)
                                                                                                                ccAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAg                                  >  3:120219/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                              <  7:301579/149‑1 (MQ=32)
                                                                                                                    cGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTa                              <  3:177550/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             <  5:232346/148‑1 (MQ=32)
                                                                                                                                      aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt             <  2:78964/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)
                                                                                                                                               ccAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAACCTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGaaa    <  3:203707/148‑1 (MQ=17)
                                                                                                                                                cAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAg  <  4:120219/149‑1 (MQ=32)
                                                                                                                                                    |                                                                                                                                                 
GAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAG  >  NC_002947/697342‑697634

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: