Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,486 C→G 8.6% intergenic (+413/‑338) 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,4860CG 8.6% 219.3 / 8.1 78intergenic (+413/‑338)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base G (3/4);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCCCCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAG  >  NC_002947/697342‑697634
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                                            <  5:89416/136‑1 (MQ=255)
                   cccACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTa                                                                                                                                            >  6:89416/1‑136 (MQ=255)
                      aCTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCt                                                                                                                            <  5:219003/149‑1 (MQ=255)
                       cTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                          >  4:62351/1‑131 (MQ=25)
                       cTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCACTTGGAAGCGCTGtcaa                                                                                                                                          <  3:62351/134‑4 (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                                                                                                                       <  8:142189/149‑1 (MQ=255)
                           cAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgag                                                                                                                       <  1:58674/149‑1 (MQ=255)
                           cAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                        <  5:145800/148‑1 (MQ=255)
                               gACAACGGAACGCCAAAATCCTTGAAATTGAATTAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCTCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGAt                                                                                                                   <  7:253582/149‑1 (MQ=21)
                                  aaCGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTccc                                                                                                                               <  3:143786/134‑1 (MQ=255)
                                  aaCGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTccc                                                                                                                               >  4:143786/1‑134 (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                                                                                                                        <  2:282200/126‑1 (MQ=255)
                                                 tCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                        >  1:282200/1‑126 (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                                                                                                           <  1:258357/100‑1 (MQ=255)
                                                                                        gTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGa                                                                                                           >  2:258357/1‑100 (MQ=255)
                                                                                         tAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAAc                                                                    <  4:92190/138‑1 (MQ=40)
                                                                                         tAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAAc                                                                    >  3:92190/1‑138 (MQ=40)
                                                                                           aGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAAc                                                        <  8:172280/148‑1 (MQ=39)
                                                                                            gggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                        >  2:208817/1‑83 (MQ=255)
                                                                                            gggCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                                        <  1:208817/83‑1 (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                                                                                                                      >  6:139041/1‑78 (MQ=255)
                                                                                                   gATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgaga                                                                                                                      <  5:139041/78‑1 (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)
                                                                                                                                             ttCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGa      >  3:135835/1‑148 (MQ=32)
                                                                                                                                                cAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAg  >  2:197359/1‑149 (MQ=32)
                                                                                                                                                cAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAg  >  1:172222/1‑149 (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 < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: