Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,105 A→C 7.7% intergenic (+29/‑722) 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,1050AC 7.7% 419.4 / 10.7 150intergenic (+29/‑722)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (67/65);  new base C (4/10);  total (71/79)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

GACCTGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGG  >  NC_002947/696957‑697251
                                                                                                                                                    |                                                                                                                                                  
gACCTGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAg                                                                                                                                                   <  1:3208744/150‑1 (MQ=255)
gACCTGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAg                                                                                                                                                   <  2:3254365/150‑1 (MQ=255)
   cTGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGcc                                                                                                                                                >  1:1495929/1‑150 (MQ=255)
   cTGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGcc                                                                                                                                                >  2:2293760/1‑150 (MQ=255)
    tGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCt                                                                                                                                               >  1:3453985/1‑150 (MQ=255)
    tGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCt                                                                                                                                               >  1:3454142/1‑150 (MQ=255)
    tGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCt                                                                                                                                               >  1:1716469/1‑150 (MQ=255)
       aTGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                            <  1:4306457/150‑1 (MQ=255)
         gCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCTGACGGCCTCCCCg                                                                                                                                          <  2:2354233/150‑1 (MQ=255)
          ccTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGACGCGGGGAGGCCGTCAGGCCTCCCCGc                                                                                                                                         <  2:2199462/150‑1 (MQ=255)
          ccTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCTGACGGCCTCCCCGc                                                                                                                                         <  1:4370690/150‑1 (MQ=255)
           cTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCt                                                                                                                                        >  2:126919/1‑150 (MQ=255)
             aTGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCttt                                                                                                                                      <  2:3421473/150‑1 (MQ=255)
                   cGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGttttt                                                                                                                                >  2:3334280/1‑150 (MQ=255)
                   cGGATGGTGTGCGCTTCCATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGttttt                                                                                                                                >  1:1799677/1‑150 (MQ=255)
                    ggttGGTGTGCGCTTCTATCCCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCTTGACGCGTCGCTGTTTGCTTTCCCTAGCCTGTACGGGAAGTATATGACGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTg                                                                                                                               <  1:4341626/147‑1 (MQ=255)
                       tGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTgggg                                                                                                                            >  2:1013309/1‑150 (MQ=255)
                               gcTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGgttggttg                                                                                                                    >  1:4138060/1‑150 (MQ=255)
                                cTTCTATACCCGTCGCACGTCCATTACCCCGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCTGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATTAAGCGGGGAGGCCGTAAGGCATCCCCGCTTTGTTTTTGGGGGTTGgttgt                                                                                                                   <  2:3124900/150‑1 (MQ=255)
                                 ttCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGgttgtt                                                                                                                  <  2:2390454/150‑1 (MQ=255)
                                 ttCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGgttgtt                                                                                                                  >  2:3706660/1‑150 (MQ=255)
                                      taCCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGt                                                                                                             <  1:4061040/150‑1 (MQ=255)
                                        ccggTCGCATGGCATTTACCCTGGGTGGGCCGCTGGGGGCAAGCATTTATGCGTCGCTTTTTGAATTCCCTCGCCTTTAAGGTATTTAGATGACGGGGGGTGGCCGTCCGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTg                                                                                                           <  1:1664502/147‑1 (MQ=16)
                                        cccGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTGGGTTGTTGCTGTTg                                                                                                           >  2:3633473/1‑150 (MQ=255)
                                         ccGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTgg                                                                                                          >  2:4498454/1‑150 (MQ=255)
                                         ccGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTgg                                                                                                          >  2:661786/1‑150 (MQ=255)
                                         ccGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTgg                                                                                                          >  1:2122500/1‑150 (MQ=255)
                                                   gCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCGTGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTg                                                                                                >  2:4307304/1‑150 (MQ=255)
                                                     aaTTCCCCAGGGTTGGCCGCTGCGGGCCAGCCATGATGGGTCGCTGTTTGCATTCCCTCGCCTGTATGGGAATTTGTTGACGGGGGGCGGCGGTCCGGCCCCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGtt                                                                                              <  1:2333644/149‑1 (MQ=255)
                                                         aCCCAGCGTTGGCCGCAGCGGGCCAGCCATTACGGGTCGCTGTTTGCATTCCCTAGCCTGTAAGGGAAGTTGATGAAGGGGGTAGGCCGTCTGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTga                                                                                          <  2:189160/150‑1 (MQ=38)
                                                           ccAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTgaga                                                                                        >  2:1561509/1‑150 (MQ=255)
                                                           ccAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTGTTGGGGGTTGGTGGTTGCTGTTGGGGAGGAACTGTTTGgaga                                                                                        >  2:4369299/1‑150 (MQ=255)
                                                             tgggTTGGCCGCCGCGGGCCATCAATTATGCGTCGCTTTTTTATTTCCCTAGCTTGTATGGTATTTTTATGATGGGGGGTGGCGGTCTGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAgg                                                                                      <  1:1235177/147‑1 (MQ=255)
                                                             tgggTTGGCCGCAGCGGGCCTGCCTTTAAGCGTCGCTTTTTGAATTCCCTAGCCTGTATGGTATTTTGTTGATGGGGGGAGGCCTTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAgg                                                                                      <  1:3198077/147‑1 (MQ=25)
                                                                 ttGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTGGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGCCGTCAGGCCTACCCGCGTTGTTTTTGGGGGTTGGTTTTTTCTGTTGGGGGGGAACTGTTTTGTGAGGgggt                                                                                  >  1:3081567/1‑150 (MQ=255)
                                                                 ttGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTGGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGGTTTGAGAGGgtgt                                                                                  >  1:1800068/1‑150 (MQ=255)
                                                                  tGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTAAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGtt                                                                                 <  1:1823094/150‑1 (MQ=255)
                                                                    gCCGCAGCGGGCCAGCCATTAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGTAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTgg                                                                               <  2:2521668/150‑1 (MQ=255)
                                                                    gCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTgg                                                                               >  1:4283122/1‑150 (MQ=255)
                                                                     ccGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGGGTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGc                                                                              >  2:3866835/1‑150 (MQ=255)
                                                                      cGCAGCGGGCCAGCCATTATGCGTCGCTGTTTGCATTCCCTAGCCTGTACGGGATGTTGTTGACGGGGGGAGGCCGTCCGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCt                                                                             <  1:2374249/150‑1 (MQ=37)
                                                                      cGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCt                                                                             >  2:3556680/1‑150 (MQ=255)
                                                                          gCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTga                                                                         <  1:2594824/150‑1 (MQ=255)
                                                                          gCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTga                                                                         <  2:2069629/150‑1 (MQ=255)
                                                                          gCGGGCCAGCATTGACGGGTCGCTTTTGGCTTCCCCTAGCCTTTACGGTATTTTTATGACGGGGGGTGGCCGTCTGGCCCCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTga                                                                         <  1:410283/150‑1 (MQ=255)
                                                                             ggCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTACGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGAt                                                                      <  2:3430833/150‑1 (MQ=255)
                                                                             ggCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGAt                                                                      >  1:2740324/1‑150 (MQ=255)
                                                                               ccagcatTGATGCGTCGCTGTTTTATTTCCCTAGCCTTTATGGTATTTTTATGATGGGGGGTGGCCTTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGAtgt                                                                    <  1:1172135/150‑1 (MQ=21)
                                                                                 agccaTGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATgtgt                                                                  >  2:3920589/1‑150 (MQ=255)
                                                                                     aTGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTACCCGCTGTGGTTTTGGGGGGTGGTTGTTGCGTTTGGGGAGGAACTGTTTTTAGAGGGTGTTGGCTGTGAGATGTGTCCGc                                                              >  1:1033011/1‑150 (MQ=255)
                                                                                      tGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGcc                                                             <  2:3620631/150‑1 (MQ=255)
                                                                                      tGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGcc                                                             >  2:1519225/1‑150 (MQ=255)
                                                                                      tGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGcc                                                             >  2:1330265/1‑150 (MQ=255)
                                                                                            gTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAAcc                                                       <  1:2624654/150‑1 (MQ=255)
                                                                                             tCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAAccc                                                      >  1:2188113/1‑150 (MQ=255)
                                                                                              cGCCGTTTGCAGTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCt                                                     <  2:1885305/150‑1 (MQ=255)
                                                                                              cGCAGTTTGCATTCCCTAGCCTGTAGGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGTTGTGTCCGCCACAACCCt                                                     >  1:1552264/1‑150 (MQ=255)
                                                                                               gCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCtt                                                    <  1:2171758/150‑1 (MQ=255)
                                                                                                ctttttGCATTCCCTCGCCTTTACGGTAAGTAGATGAAGCGGGGCGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTg                                                   <  1:3777458/147‑1 (MQ=25)
                                                                                                  gTTTGCATTCCCTAGCCTGTAAGGGAAGTTGTTGATGCGGGGTGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTGGGCTGTGAGATGTGTCCGCCACAACCCTTGCa                                                 <  2:473607/150‑1 (MQ=25)
                                                                                                      gCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGc                                             <  2:3109503/150‑1 (MQ=255)
                                                                                                       cATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGcc                                            <  1:3552154/150‑1 (MQ=255)
                                                                                                       cATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGcc                                            >  2:2130036/1‑150 (MQ=255)
                                                                                                            ccTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTga                                       <  2:1716469/150‑1 (MQ=255)
                                                                                                            ccTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTga                                       >  2:4424538/1‑150 (MQ=255)
                                                                                                            ccTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTga                                       <  2:684767/150‑1 (MQ=255)
                                                                                                            ccTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCACCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACAGTTTTTAGAGGGTGTTGGCAGTGAGATTTGTCCGCCACAAACCATTCAGTGACTGTga                                       >  1:562136/1‑150 (MQ=255)
                                                                                                              tAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCATGCCTACCCGCTTTGTTTTTTGGGGTTTGTTGTTTCTGTTTGGGAGGTACTGTTTTTAGATGGTGTTGTGTGTGTGATTTTTACGCCACAAACCTTGCAGTTGCTGTgagt                                     >  1:2347011/1‑149 (MQ=255)
                                                                                                              tAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTgaga                                     >  2:2226562/1‑150 (MQ=255)
                                                                                                              tAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTgaga                                     >  2:4197627/1‑150 (MQ=255)
                                                                                                              tAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCCTCAGGCCTCCCCCCTTTGGTTTTGGGGGGTGGTTGTTGCTGTTGGGGAGGAACAGTTTTGAGAGGGTGTTGGCGGGGAGATGTGTCCGCCACAACCCATGCAGTGCCAGTgaga                                     >  1:4180975/1‑150 (MQ=25)
                                                                                                                gCCTGTATGGCAAGTAGATGACGGGGGGAGGCGGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATc                                   <  1:2296609/150‑1 (MQ=21)
                                                                                                                gCCTGTACGGTAAGTAGATGACGCGGGGCGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATc                                   <  1:3043069/150‑1 (MQ=255)
                                                                                                                  cTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGa                                 <  2:1025324/150‑1 (MQ=255)
                                                                                                                  cTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTATTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGa                                 <  2:3542554/150‑1 (MQ=255)
                                                                                                                   tGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAg                                <  2:2188113/150‑1 (MQ=255)
                                                                                                                    gTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGaa                                                                                                   <  1:3732553/82‑1 (MQ=255)
                                                                                                                    gTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGaa                                                                                                   >  2:3732553/1‑82 (MQ=255)
                                                                                                                    gTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAgc                               <  2:3454142/150‑1 (MQ=255)
                                                                                                                     tAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAgcg                              >  1:3762927/1‑150 (MQ=255)
                                                                                                                          gAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCgccgcc                         >  2:4219232/1‑150 (MQ=255)
                                                                                                                              tAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGGTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCAGTGAGATCGAGCGCCGCCcgcg                     >  2:1714913/1‑150 (MQ=255)
                                                                                                                               ataTGACGGGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgc                    <  1:1972853/148‑1 (MQ=255)
                                                                                                                                  tGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcggc                 <  2:3803066/150‑1 (MQ=255)
                                                                                                                                     cgCGGGGAGGCCGTCCGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt              <  1:2420879/149‑1 (MQ=255)
                                                                                                                                      gCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc             >  1:1838831/1‑150 (MQ=255)
                                                                                                                                        ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg           >  2:3010007/1‑150 (MQ=255)
                                                                                                                                        ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGGTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg           >  2:3009654/1‑150 (MQ=18)
                                                                                                                                        ggggAGGCCTGACGGCCACCCCGCTTTGTTTTTGGGGGGTGGTTGTTGCTGTTTGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg           >  1:860062/1‑150 (MQ=18)
                                                                                                                                        ggggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg           <  2:2298779/150‑1 (MQ=255)
                                                                                                                                        ggggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGGTGGGGAGGACCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgcggcgcgcgg           >  2:1999079/1‑150 (MQ=255)
                                                                                                                                        ggggAGGCCGTCAGGCCTCCCCGCTTTGGTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg           >  2:4405971/1‑150 (MQ=255)
                                                                                                                                          ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACGACAGTTTCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtc         <  2:2122500/150‑1 (MQ=255)
                                                                                                                                          ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGAtc         <  1:2830128/150‑1 (MQ=18)
                                                                                                                                            aGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtctc       <  1:2419802/150‑1 (MQ=255)
                                                                                                                                             ggCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTca      <  2:2498545/150‑1 (MQ=255)
                                                                                                                                              gCCGTCTGGCCTCCCCGCTTTGTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTcaca    <  2:737298/150‑1 (MQ=255)
                                                                                                                                              gCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTcac     >  2:859797/1‑150 (MQ=255)
                                                                                                                                                 gTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAgg  >  2:2221814/1‑150 (MQ=255)
                                                                                                                                                 gTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAgg  >  1:1002889/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GACCTGCATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGG  >  NC_002947/696957‑697251

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: