Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,481 G→A 7.4% 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 7.4% 220.8 / 6.1 81intergenic (+408/‑343)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (33/42);  new base A (2/4);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

CCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGA  >  NC_002947/697335‑697627
                                                                                                                                                    |                                                                                                                                                  
cccTCTTGAATCCCCTTATACTGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCt                                                                                                                                                  <  8:51197/149‑1 (MQ=255)
   tctTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCACttg                                                                                                                                                >  5:247949/1‑147 (MQ=255)
    ctggAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACATACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTcc                                                                                                                                               <  5:235790/145‑1 (MQ=255)
    ctTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAAAGAGTTAGACAATTCTGAAGGTCCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑CTGGAGATGAAGGTTGACAGCGCTAcc                                                                                                                                               >  6:235790/1‑148 (MQ=25)
     ttGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGGGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTcc                                                                                                                                               <  6:242904/147‑1 (MQ=255)
             ccTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtc                                                                                                                                     <  6:247949/149‑2 (MQ=255)
                 ataATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaac                                                                                                                                  <  5:409053/148‑5 (MQ=21)
                 ataATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                   <  5:55793/147‑4 (MQ=25)
                          cccACTTCCAGCGACAACGTAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGAtt                                                                                                                                <  5:63585/141‑1 (MQ=255)
                          cccACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGAtt                                                                                                                                >  6:63585/1‑141 (MQ=255)
                            cACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGc                                                                                                                      <  7:154780/149‑1 (MQ=255)
                             aCTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGGTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTa                                                                                                                                     >  6:157110/1‑133 (MQ=255)
                             aCTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTa                                                                                                                                     <  5:157110/133‑1 (MQ=255)
                               ttCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGATTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTAc                                                                                                                   <  4:21786/149‑1 (MQ=255)
                                   aGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgaga                                                                                                               <  8:206459/149‑1 (MQ=255)
                                    gCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgagag                                                                                                              <  5:277789/149‑1 (MQ=255)
                                    gCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                 >  3:302083/1‑146 (MQ=255)
                                    gCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                 >  3:296740/1‑146 (MQ=255)
                                      gACAACGGAACGACAAATTCCTTGAAATTCAATGAGTTAGACAATACTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGAt                                                                                                            <  2:186256/149‑1 (MQ=255)
                                       aCAAAGGAACGAAAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATc                                                                                                           <  4:225062/149‑1 (MQ=255)
                                         aaCGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCg                                                                                                          >  4:336715/1‑148 (MQ=255)
                                            ggAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGc                                                                                                      <  3:208115/149‑1 (MQ=255)
                                                 gACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGc                                                                                                      >  3:194656/1‑144 (MQ=255)
                                                 gACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGc                                                                                                      <  4:194656/144‑1 (MQ=255)
                                                 gACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGc                                                                                                  >  6:416220/1‑148 (MQ=255)
                                                   cAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGACTCCCGCTACGAGAGATCGCAGCGAGcc                                                                                                 >  4:116261/1‑147 (MQ=255)
                                                   cAAACTCCTTGAAATTCAATGAGTTAGACAAATCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGcc                                                                                                 >  2:77121/1‑147 (MQ=255)
                                                    aaaCTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAg                                                                                              >  7:254225/1‑149 (MQ=255)
                                                        tCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGtt                                                                                          >  2:31432/1‑149 (MQ=255)
                                                           ttGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTg                                                                                        >  5:261400/1‑148 (MQ=255)
                                                            tGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGAGTTGAGATGATGGTTGACAGCGCTTCCAAGTGATGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                         >  1:120351/1‑148 (MQ=255)
                                                             gAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGaa                                                                                      <  7:182183/148‑1 (MQ=255)
                                                              aaaTTCAATGAGTTATAGAATTCTGAACGTGCTGTAAGGTCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGTTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGc                                                                                    <  1:75855/149‑1 (MQ=17)
                                                                aTTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCt                                                                                   >  5:112113/1‑148 (MQ=255)
                                                                    aaTGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAAc                                                                               <  4:302083/148‑1 (MQ=255)
                                                                    aaTGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAAc                                                                               <  4:296740/148‑1 (MQ=255)
                                                                     aTGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGa                                                                             <  8:254225/149‑1 (MQ=255)
                                                                      tgtgTTAGACTATTCTGAAGGTGATGTAAGGTCTACTTTCAAATGATCGTCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTTTATGATTCGCCTCCCGCTAAGAGAGATGCCAGCGAGCCAAGTGTTTGAAGCTAAACGAg                                                                            <  1:34084/146‑1 (MQ=17)
                                                                      tGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAg                                                                            >  5:336514/1‑149 (MQ=255)
                                                                          ttAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTtc                                                                        <  2:334345/149‑1 (MQ=255)
                                                                             gACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTtctc                                                                      <  5:416220/148‑1 (MQ=255)
                                                                              aCAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCg                                                                     >  2:321211/1‑148 (MQ=255)
                                                                               cAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGc                                                                    >  8:212214/1‑148 (MQ=255)
                                                                               cAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGATAAACGAGTTTCTCGc                                                                    >  8:206455/1‑148 (MQ=255)
                                                                                aaTTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCa                                                                   >  5:175502/1‑148 (MQ=255)
                                                                                aaTTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATCCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCa                                                                   >  8:108247/1‑148 (MQ=255)
                                                                                       aaGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCaaaaaa                                                              <  7:14499/146‑1 (MQ=255)
                                                                                       aaGGGGCTGTAAGGCCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTCCTGTATGATTCGCCTACCGCTACGAGATATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACtt                                                           <  8:74262/149‑1 (MQ=14)
                                                                                         ggTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATCCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCa                                                         <  7:108247/149‑1 (MQ=255)
                                                                                            gCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACg                                                                                                                  >  1:361746/1‑89 (MQ=255)
                                                                                            gCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACg                                                                                                                  <  2:361746/89‑1 (MQ=255)
                                                                                                  aGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACtt                                                           >  5:254737/1‑138 (MQ=40)
                                                                                                  aGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACtt                                                           <  6:254737/138‑1 (MQ=40)
                                                                                                   gggCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGc                                               >  6:345462/1‑149 (MQ=39)
                                                                                                         cGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGa                                           >  5:252055/1‑147 (MQ=39)
                                                                                                         cGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACCAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACa                                         <  4:390589/149‑1 (MQ=255)
                                                                                                          gATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAg                                        >  6:303736/1‑149 (MQ=39)
                                                                                                           aTCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAATCGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAgg                                       >  7:291238/1‑149 (MQ=255)
                                                                                                           aTCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAgg                                       <  6:261400/149‑1 (MQ=39)
                                                                                                               aaTGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGctct                                   >  8:345830/1‑149 (MQ=37)
                                                                                                                   aTCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAg                                <  5:303736/148‑1 (MQ=37)
                                                                                                                      gCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGaaa                            <  5:345462/149‑1 (MQ=34)
                                                                                                                        cAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGc                          <  6:336514/149‑1 (MQ=32)
                                                                                                                          gCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAAAGCTTGACAGGCTCTGAGGAAAGCGt                        <  2:116698/149‑1 (MQ=17)
                                                                                                                           cG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAATCGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTa                       <  8:291238/149‑1 (MQ=17)
                                                                                                                           cG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTa                       <  1:185622/149‑1 (MQ=32)
                                                                                                                               gTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGa                     <  2:325495/149‑1 (MQ=32)
                                                                                                                                  gagaTGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCGCTGAGGAAAGCGTAGAAt                   <  5:37372/148‑1 (MQ=17)
                                                                                                                                    gatgatGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATgcg                >  5:252905/1‑149 (MQ=32)
                                                                                                                                       gatgGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATgcgcgc             >  2:264631/1‑149 (MQ=32)
                                                                                                                                       gatgGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATgcgcgc             >  8:49053/1‑149 (MQ=32)
                                                                                                                                           gTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCg         <  1:77121/149‑1 (MQ=32)
                                                                                                                                            ttGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCTAAATAAACGCTTGACAGGCTCTGAGGaaa                            <  1:237528/129‑1 (MQ=17)
                                                                                                                                            ttGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGaaa                            >  2:237528/1‑129 (MQ=32)
                                                                                                                                             tGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGGTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAAc                                                 >  5:248865/1‑107 (MQ=17)
                                                                                                                                             tGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGCTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAAc                                                 <  6:248865/107‑1 (MQ=17)
                                                                                                                                             ataCAGCACTTGGAAGCGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGt       >  4:143362/3‑149 (MQ=255)
                                                                                                                                              taCAGCACTTGGAAGCGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt      <  3:143362/148‑1 (MQ=255)
                                                                                                                                               aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTg     <  6:112113/149‑1 (MQ=32)
                                                                                                                                                cAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTga    <  6:175502/149‑1 (MQ=32)
                                                                                                                                                cAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTga    <  6:252905/149‑1 (MQ=32)
                                                                                                                                                  gcgcTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTgaga  >  8:32338/1‑149 (MQ=32)
                                                                                                                                                  gcgcTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTgaga  <  6:252055/149‑1 (MQ=32)
                                                                                                                                                    |                                                                                                                                                  
CCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCG‑‑GTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGA  >  NC_002947/697335‑697627

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: