Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,214 T→C 6.2% intergenic (+141/‑610) 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,2140TC 6.2% 237.2 / 4.7 81intergenic (+141/‑610)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (39/37);  new base C (3/2);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

TGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCC  >  NC_002947/697069‑697362
                                                                                                                                                 |                                                                                                                                                       
tGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATc                                                                                                                                                        >  5:376030/1‑147 (MQ=255)
 gTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAg                                                                                                                                                     >  7:179338/1‑149 (MQ=255)
             gATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgc                                                                                                                                         >  5:333045/1‑149 (MQ=255)
               tGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                       <  4:162335/149‑1 (MQ=255)
                 atgCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcggc                                                                                                                                      <  2:318025/146‑1 (MQ=255)
                      gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                <  4:137241/149‑2 (MQ=255)
                      gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                >  3:137241/1‑148 (MQ=255)
                           gCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTca                                                                                                                           <  2:276200/149‑1 (MQ=255)
                             cGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTca                                                                                                                           >  5:90927/1‑147 (MQ=255)
                              gTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGCGCGTGTTGGCTGTGCGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑Ag                                                                                                                     <  5:396632/149‑1 (MQ=255)
                              gTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑Ag                                                                                                                     >  1:199077/1‑149 (MQ=255)
                                      tCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCgat                                                                                                                               <  7:7370/134‑3 (MQ=255)
                                      tCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCgat                                                                                                                               >  8:7370/1‑132 (MQ=255)
                                        cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAAc                                                                                                           >  1:211761/1‑149 (MQ=255)
                                         ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcg                                                                                                                                          <  6:271088/120‑1 (MQ=255)
                                         ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcg                                                                                                                                          >  5:271088/1‑120 (MQ=255)
                                         ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACt                                                                                                          <  6:90927/149‑1 (MQ=255)
                                         ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACt                                                                                                          >  2:69005/1‑149 (MQ=255)
                                            cTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTc                                                                                                       <  1:69005/149‑1 (MQ=255)
                                             tttGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCg                                                                                                      >  2:136746/1‑149 (MQ=255)
                                             tttGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCg                                                                                                      >  4:24248/1‑149 (MQ=255)
                                               tGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCg                                                                                                      <  6:376030/147‑1 (MQ=255)
                                               tGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTg                                                                                                    >  2:177890/1‑149 (MQ=255)
                                                      gggggTTGGTTGTTGCTGTTGGGGAGGAAATGTTTTGAGAGGGTTTTGGCTGTGAGATTTGTCCGCAACAACCATTGCAGTGCCTGTCAGATCCAGAGCCGCCCGCGCGGCGATCGTTCTAAC‑‑AGGCGCTCAAAATTACGTGGAGAac                                                                                              >  5:302356/1‑148 (MQ=255)
                                                       ggggTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAacac                                                                                            >  2:29753/1‑149 (MQ=255)
                                                           ttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                        <  1:168599/149‑1 (MQ=255)
                                                           ttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                        <  1:177890/149‑1 (MQ=255)
                                                             ggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                  <  8:47802/137‑1 (MQ=255)
                                                             ggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                  >  7:47802/1‑137 (MQ=255)
                                                              gttgttGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAg                                                                                     >  7:26746/1‑149 (MQ=255)
                                                               ttgttgCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                        >  7:126799/1‑145 (MQ=255)
                                                               ttgttgCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                        <  8:126799/145‑1 (MQ=255)
                                                                    gCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCg                                                                                >  3:1463/1‑148 (MQ=255)
                                                                        ttGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATAc                                                                            <  2:199077/148‑1 (MQ=255)
                                                                        ttGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACg                                                                           <  4:386683/149‑1 (MQ=255)
                                                                          ggggAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACgg                                                                          <  6:302356/148‑1 (MQ=255)
                                                                             gaggaACTGTTTTGAGAGGGTGTGTGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCtt                                                                      <  3:241706/149‑1 (MQ=255)
                                                                               ggaACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCttt                                                                     >  5:35110/1‑148 (MQ=255)
                                                                                gaACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTc                                                                   >  1:405326/1‑149 (MQ=255)
                                                                                  aCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTccc                                                                 <  6:35110/149‑1 (MQ=255)
                                                                                        ttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                          >  7:365160/1‑118 (MQ=255)
                                                                                        ttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                          <  8:365160/118‑1 (MQ=255)
                                                                                         tGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                  <  5:342800/109‑1 (MQ=255)
                                                                                         tGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                  >  6:342800/1‑109 (MQ=255)
                                                                                          gagaGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                >  5:21014/1‑80 (MQ=17)
                                                                                          gagaGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                <  6:21014/81‑2 (MQ=17)
                                                                                                   ttGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                  <  3:374348/99‑1 (MQ=255)
                                                                                                   ttGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                  >  4:374348/1‑99 (MQ=255)
                                                                                                   ttGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCt                                                <  8:179338/149‑1 (MQ=255)
                                                                                                          tgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCAGTTTCTTGAAAt                                          >  6:276365/1‑148 (MQ=255)
                                                                                                          tgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAAt                                          >  5:371812/1‑148 (MQ=255)
                                                                                                            ataTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAA‑‑AGGCGCTGATACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGGCTTTTCCCCTTATATTTCACTTTCTTGAAATTa                                        >  5:293187/3‑148 (MQ=255)
                                                                                                            agaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTc                                                       <  7:115184/133‑1 (MQ=255)
                                                                                                            agaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTc                                                       >  8:115184/1‑133 (MQ=255)
                                                                                                             gaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                      <  2:211761/149‑1 (MQ=255)
                                                                                                              aTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtctc                                                                                                                            >  3:56497/1‑65 (MQ=255)
                                                                                                              aTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtctc                                                                                                                            <  4:56497/65‑1 (MQ=255)
                                                                                                                  gtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                           >  6:369244/1‑139 (MQ=255)
                                                                                                                  gtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                           <  5:369244/139‑1 (MQ=255)
                                                                                                                   tcggCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCTGCACTCGGTCTCAC‑‑AGTCTCTGAAACTGTCGTGGCGAACACTTgt                                                                                        <  1:68938/90‑2 (MQ=25)
                                                                                                                   tCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCCTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGa                                >  2:187333/1‑149 (MQ=255)
                                                                                                                   tCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCACTCGGTCTCAC‑‑AGGCGCTGAAACTGACGTGGCGAACACTTgg                                                                                        >  2:68938/1‑93 (MQ=255)
                                                                                                                      gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                  <  3:250412/144‑1 (MQ=255)
                                                                                                                      gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                  >  4:250412/1‑144 (MQ=255)
                                                                                                                       ccACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATATGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                           <  4:218011/134‑1 (MQ=255)
                                                                                                                       ccACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATAAGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                           >  3:218011/1‑134 (MQ=255)
                                                                                                                           aaCCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATAttt                                                        <  7:378001/117‑1 (MQ=255)
                                                                                                                           aaCCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTACCCCTTATAttt                                                        >  8:378001/1‑117 (MQ=255)
                                                                                                                              ccTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTg                      >  4:337446/1‑148 (MQ=255)
                                                                                                                                tttcAGCGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCAC‑‑AGGCACTGCAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGATATTAAAGGTTGACACCCTCTTGAAt                   >  5:282013/4‑149 (MQ=11)
                                                                                                                                   cAGTGCCTGTGAGATCGAGCGCCGCCCGCGTGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATccc                <  8:126107/149‑1 (MQ=255)
                                                                                                                                   cAGTGCCTGTGAGATCGAGCGCCGCCCGCGTGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATccc                >  7:126107/1‑149 (MQ=255)
                                                                                                                                    aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTccc                                                                 >  7:167471/1‑99 (MQ=255)
                                                                                                                                    aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTccc                                                                 <  8:167471/99‑1 (MQ=255)
                                                                                                                                    aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc               >  4:312104/1‑149 (MQ=255)
                                                                                                                                      tGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCtt             <  3:24248/149‑1 (MQ=255)
                                                                                                                                        ccTGTGATATCGAGCGCCGCCCGCGCGGCGCTCGGTTTCACTAATGGGGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc               <  1:136746/148‑1 (MQ=14)
                                                                                                                                            tgAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAAt        <  2:405326/148‑1 (MQ=255)
                                                                                                                                            tgAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCCTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATg       <  1:187333/149‑1 (MQ=255)
                                                                                                                                               gaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgcgc    <  8:26746/149‑1 (MQ=255)
                                                                                                                                                aTCGAGAGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑ATGCGCTGAAACTGTCGTGGCGAACACTTGGCAGGCCCGATACGGTCTTTTCCCCTTATATTTAACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATACCCTTATAATGCGcc   >  7:39485/1‑149 (MQ=14)
                                                                                                                                                 tCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGccc  <  6:333045/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
TGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCC  >  NC_002947/697069‑697362

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: