Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,490 T→C 6.5% 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 6.5% 219.2 / 4.1 79intergenic (+417/‑334)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (34/38);  new base C (2/3);  total (38/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

CCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAGC  >  NC_002947/697348‑697635
                                                                                                                                              |                                                                                                                                                 
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                                                                                                                                           <  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)
                     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                                                                                                                                >  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                                                                                                                       >  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   >  1:172222/1‑149 (MQ=32)
                                                                                                                                          cAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAg   >  2:197359/1‑149 (MQ=32)
                                                                                                                                           aaGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAGc  <  7:140243/149‑1 (MQ=32)
                                                                                                                                              |                                                                                                                                                 
CCTTATAATGCGCCCCACTTCCAGCGACAACGGAACGACAAACTCCTTGAAATTCAATGAGTTAGACAATTCTGAAGGTGCTGTAAGGGCTACGATCAAATGATCGCCAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAAGC  >  NC_002947/697348‑697635

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: