Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,481 G→A 6.7% 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 6.7% 214.0 / 4.1 75intergenic (+408/‑343)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base A (2/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGAC  >  NC_002947/697336‑697628
                                                                                                                                                 |                                                                                                                                                   
ccTCTTGAATCCCCTTATAATGAGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATCATCGCCAGCGGGTGAGATGATGGTTGACAGCGCtt                                                                                                                                                  >  5:254141/1‑149 (MQ=255)
      gAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTTTAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCaa                                                                                                                                              <  7:299459/147‑1 (MQ=255)
            ccTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTa                                                                                                                                      <  1:227529/149‑1 (MQ=255)
             cTTATAATGCGCCCCACTTCCTGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACTATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGTTGGTTGACGGCGCTTCCATGTGCTGTAt                                                                                                                                     >  4:30653/1‑149 (MQ=21)
                ataATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                    <  3:64579/147‑4 (MQ=25)
                ataATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAAAGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGAGCTTCGAAGTGCTGTATGAt                                                                                                                                  >  8:253582/1‑149 (MQ=255)
                   aTGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATCCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTc                                                                                                                                <  4:10039/148‑1 (MQ=255)
                   aTGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATCCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTc                                                                                                                                >  3:10039/1‑148 (MQ=255)
                    tGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGTGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGc                                                                                                                              >  6:145800/1‑149 (MQ=255)
                         cccACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTa                                                                                                                                      >  6:89416/1‑136 (MQ=255)
                         cccACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTa                                                                                                                                      <  5:89416/136‑1 (MQ=255)
                            aCTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCt                                                                                                                      <  5:219003/149‑1 (MQ=255)
                             cTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                    <  3:62351/134‑4 (MQ=25)
                             cTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                    >  4:62351/1‑131 (MQ=25)
                              ttCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTAc                                                                                                                    <  1:199729/149‑1 (MQ=255)
                              ttCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                    >  4:199038/1‑130 (MQ=25)
                              ttCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                    <  3:199038/133‑4 (MQ=25)
                                 cAGCGACAACGGAACGACAAACTCCTTGAAATTCATTGAGTTAGACAATTCTGAAGGTGCTGTAAGGGATACGATCAAATGATATCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCaa                                                                                                                                              >  8:206724/1‑120 (MQ=37)
                                 cAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgag                                                                                                                 <  1:58674/149‑1 (MQ=255)
                                 cAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgag                                                                                                                 <  8:142189/149‑1 (MQ=255)
                                 cAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                  <  5:145800/148‑1 (MQ=255)
                                     gACAACGGAACGCCAAAATCCTTGAAATTGAATTAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCTCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGAt                                                                                                             <  7:253582/149‑1 (MQ=21)
                                        aaCGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTccc                                                                                                                         >  4:143786/1‑134 (MQ=255)
                                        aaCGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTccc                                                                                                                         <  3:143786/134‑1 (MQ=255)
                                           ggAACGACAAACTCCTTGAAATTCAATGAGTCAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATagcagc                                                                                                       >  1:194208/1‑149 (MQ=255)
                                           ggAACGACAAACTACTTGAAATTCAATGAGTTAGACAAATCTGAAGCTGCTGTAAGGGCTACGATCAAATGATAGCCAGCGGTTGAGATGATGGTTGACAGCTCTTACAAGTGCTGTATGATTCGCCTCCCTCTACGAGAGATCGCat                                                                                                        >  8:213420/1‑147 (MQ=21)
                                                gACAAACTCCTTTAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGcc                                                                                                  <  3:30653/149‑1 (MQ=255)
                                                gACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGcc                                                                                                  <  1:88768/149‑1 (MQ=255)
                                                   aaaCTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAg                                                                                               >  4:148263/1‑149 (MQ=255)
                                                    aaCTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAgt                                                                                              <  7:85583/149‑1 (MQ=255)
                                                      cTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAgtgt                                                                                            >  6:164502/1‑149 (MQ=255)
                                                       tCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                  >  1:282200/1‑126 (MQ=255)
                                                       tCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                  <  2:282200/126‑1 (MQ=255)
                                                          ttGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTg                                                                                         >  7:193837/1‑148 (MQ=255)
                                                          ttGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATAGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTg                                                                                         >  6:74668/1‑148 (MQ=255)
                                                            gAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGTGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGAAAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGaa                                                                                       <  7:213420/148‑1 (MQ=255)
                                                            gAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTg                                                                                         >  3:190837/1‑146 (MQ=255)
                                                             aaaTTCAATGAGTTAGACAATTATGAATGTGCTGTAATGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAg                                                                                      <  6:259913/148‑1 (MQ=255)
                                                                    aTGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAAc                                                                                >  7:260109/1‑147 (MQ=255)
                                                                       aGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGt                                                                            <  3:148263/148‑1 (MQ=255)
                                                                           aGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                          >  1:157914/1‑130 (MQ=255)
                                                                           aGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGttt                                                                                          <  2:157914/130‑1 (MQ=255)
                                                                            gACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTtctc                                                                       <  6:67418/148‑1 (MQ=255)
                                                                             aCAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCg                                                                      >  7:172280/1‑148 (MQ=255)
                                                                                     gAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACt                                                             >  5:216706/1‑149 (MQ=255)
                                                                                        ggTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCa                                                          <  4:190837/149‑1 (MQ=255)
                                                                                         gTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCGAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCa                                                          <  6:254141/148‑1 (MQ=255)
                                                                                              gTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGa                                                                                                     >  2:258357/1‑100 (MQ=255)
                                                                                              gTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGa                                                                                                     <  1:258357/100‑1 (MQ=255)
                                                                                               tAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAAc                                                              >  3:92190/1‑138 (MQ=40)
                                                                                               tAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAAc                                                              <  4:92190/138‑1 (MQ=40)
                                                                                                 aGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAAc                                                  <  8:172280/148‑1 (MQ=39)
                                                                                                  gggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                  <  1:208817/83‑1 (MQ=255)
                                                                                                  gggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                  >  2:208817/1‑83 (MQ=255)
                                                                                                   ggCTACGATCAAATGATCGCCAGCTGTTTAGATTATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCt                                               <  6:273792/149‑1 (MQ=255)
                                                                                                   ggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTa                                                                                   <  1:311754/113‑1 (MQ=39)
                                                                                                   ggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTa                                                                                   >  2:311754/1‑113 (MQ=39)
                                                                                                    gCTACGATCAAATGATCTCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCtt                                              <  5:216633/149‑1 (MQ=255)
                                                                                                     cTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTg                                             >  6:58322/1‑149 (MQ=39)
                                                                                                        cGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGa                                            <  5:74668/147‑1 (MQ=39)
                                                                                                        cGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACa                                          <  8:193837/149‑1 (MQ=39)
                                                                                                         gATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgaga                                                                                                                <  5:139041/78‑1 (MQ=255)
                                                                                                         gATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgaga                                                                                                                >  6:139041/1‑78 (MQ=255)
                                                                                                            cAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGAGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGct                                      >  8:46861/1‑149 (MQ=38)
                                                                                                             aaaTGATCGCCAGCGGTTGAGATGATGGTTCACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGAGAGCCAAGTGTTTGAAGCTAAACGAGTTTCACGCAAAAAACTTAAAAATAAACTCTAGACAGGct                                      >  3:62145/1‑148 (MQ=11)
                                                                                                              aaTGTTCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCATGTGTTTGAAGCTAAACGAGTTTCTCGCAATAAACTTCAAAAt                                                      <  1:143916/131‑1 (MQ=255)
                                                                                                              aaTGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGAGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGATTCTCGCAAAAAACTACAAAAt                                                      >  2:143916/1‑131 (MQ=255)
                                                                                                               aTGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCACAATAAACGCTTGACAGGCTCTg                                   <  7:46861/149‑1 (MQ=37)
                                                                                                               aTGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGctct                                    >  2:200142/1‑148 (MQ=37)
                                                                                                                  aTCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAg                                 <  5:58322/148‑1 (MQ=37)
                                                                                                                     gCCAGCTGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGaaa                             <  3:23313/149‑1 (MQ=21)
                                                                                                                          cGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTa                        <  2:194208/149‑1 (MQ=32)
                                                                                                                                        gTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCg          >  2:190947/1‑149 (MQ=32)
                                                                                                                                            aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt       >  5:68614/1‑148 (MQ=32)
                                                                                                                                              aGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt       >  8:140243/1‑146 (MQ=32)
                                                                                                                                                cgcTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCATCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTa                        <  7:232944/127‑1 (MQ=17)
                                                                                                                                                cgcTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTa                        >  8:232944/1‑127 (MQ=32)
                                                                                                                                                 gcTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGAc  >  6:40661/1‑148 (MQ=32)
                                                                                                                                                 |                                                                                                                                                   
CCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGAC  >  NC_002947/697336‑697628

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