Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,228 C→G 8.9% intergenic (+155/‑596) 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,2280CG 8.9% 246.3 / 13.1 79intergenic (+155/‑596)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base G (4/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

GATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGC  >  NC_002947/697082‑697372
                                                                                                                                                  |                                                                                                                                                   
gATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgc                                                                                                                                                   >  5:333045/1‑149 (MQ=255)
  tGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                                 <  4:162335/149‑1 (MQ=255)
    atgCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcggc                                                                                                                                                <  2:318025/146‑1 (MQ=255)
         gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                          >  3:137241/1‑148 (MQ=255)
         gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                          <  4:137241/149‑2 (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                                                                                                            >  6:342800/1‑109 (MQ=255)
                                                                            tGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                            <  5:342800/109‑1 (MQ=255)
                                                                             gagaGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                          >  5:21014/1‑80 (MQ=17)
                                                                             gagaGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                          <  6:21014/81‑2 (MQ=17)
                                                                                      ttGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                            >  4:374348/1‑99 (MQ=255)
                                                                                      ttGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGc                                                                                                            <  3:374348/99‑1 (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                                                     <  5:369244/139‑1 (MQ=255)
                                                                                                     gtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                                     >  6:369244/1‑139 (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                          >  7:126107/1‑149 (MQ=255)
                                                                                                                      cAGTGCCTGTGAGATCGAGCGCCGCCCGCGTGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATccc                          <  8:126107/149‑1 (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)
                                                                                                                                      gAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcc                           <  2:345874/132‑1 (MQ=255)
                                                                                                                                      gAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCa          >  1:384287/1‑149 (MQ=255)
                                                                                                                                            cgcccgcGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTc      >  8:229585/1‑147 (MQ=255)
                                                                                                                                              cccgcGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGc  <  3:312104/149‑1 (MQ=255)
                                                                                                                                              cccgcGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGc  <  2:384287/149‑1 (MQ=255)
                                                                                                                                               ccgcGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGc  >  3:67185/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCAC‑‑AGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGC  >  NC_002947/697082‑697372

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: