Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 391,576:1 +G 6.3% coding (1850/3015 nt) soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947391,5761.G 6.3% 317.8 / 14.5 79coding (1850/3015 nt)soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base . (37/37);  new base G (3/2);  total (40/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 = 6.45e-01

ATCAACGGCCTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCT  >  NC_002947/391428‑391723
                                                                                                                                                       |                                                                                                                                                   
aTCAACGGCCTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAc                                                                                                                                                      >  8:267957/1‑149 (MQ=255)
  cAACGGCCTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTg                                                                                                                                                   >  1:308193/1‑149 (MQ=255)
  cAACGGCCTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTg                                                                                                                                                   >  4:48352/1‑149 (MQ=255)
         cTGGCCATCGCCGCCCGTTCGATCGGCATCACAATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCa                                                                                                                                            >  5:106960/1‑149 (MQ=255)
          tGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCAt                                                                                                                                           >  5:115138/1‑149 (MQ=255)
           ggCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATg                                                                                                                                          >  7:31368/1‑149 (MQ=255)
            gCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCAt                                                                                                                                           <  4:257967/147‑1 (MQ=255)
            gCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATg                                                                                                                                          <  3:89808/148‑1 (MQ=255)
            gCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGc                                                                                                                                         <  7:258831/149‑1 (MQ=255)
              cATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCt                                                                                                                                       >  5:80924/1‑149 (MQ=255)
                  gccgccCGTTCGATAGGCCTCACCATCCCCGTAATGGGTTACACCATGTTTCGCCCCAAATACACGCCGGTTACTTTTGGGCTGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTgg                                                                                                                                     <  1:102911/149‑1 (MQ=255)
                          ttCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTCCACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACGTGCCAGGCATGCAGGGCAGTgaagc                                                                                                                           <  8:90182/149‑2 (MQ=255)
                               tCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGGTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGgc                                                                                                                      >  8:304342/1‑149 (MQ=255)
                                 ggCATCACCATCCCGGAAAGGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGGCACCTGTTCGAGCCCGTGCGCTTCACGTGCCAGGCATGCAGGGCAGTgaagcgcacgg                                                                                                                     <  1:155988/148‑8 (MQ=37)
                                   cATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGcc                                                                                                                   >  5:248130/1‑148 (MQ=255)
                                    aTCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGa                                                                                                                 <  4:169103/149‑1 (MQ=255)
                                       aCCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGtt                                                                                                              >  3:124082/1‑149 (MQ=255)
                                       aCCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGtt                                                                                                              >  3:129081/1‑149 (MQ=255)
                                         cATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACg                                                                                                                        <  3:289859/137‑1 (MQ=255)
                                         cATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACg                                                                                                                        >  4:289859/1‑137 (MQ=255)
                                           tCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGaa                                                                                                          <  3:48352/149‑1 (MQ=255)
                                             ccGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGa                                                                                                        <  7:267957/149‑1 (MQ=255)
                                              cGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGAc                                                                                                       >  7:169283/1‑149 (MQ=255)
                                                      ggTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCa                                                                                               <  6:133833/149‑1 (MQ=255)
                                                        tACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGt                                                                                             >  3:72823/1‑149 (MQ=255)
                                                           caccaTGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                           >  3:1228/1‑148 (MQ=255)
                                                           caccaTGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                           >  3:335352/1‑148 (MQ=255)
                                                           caccaTGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                           >  6:182768/1‑148 (MQ=255)
                                                            accaTGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGGTGTCACCTGTTCGAGCCCTTGCGCTTCACTGCCCTGCATGCCTGTCACGTGAAGAACGGCGCCTAGTTCGAAGACGTCGGCCAGTGGaa                                                                                         <  6:342164/149‑1 (MQ=255)
                                                              caTGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAgc                                                                                       >  7:154099/1‑149 (MQ=255)
                                                                  ttCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGcc                                                                                    <  6:248130/148‑1 (MQ=255)
                                                                      gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                               <  4:129081/149‑1 (MQ=255)
                                                                      gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                               <  4:124082/149‑1 (MQ=255)
                                                                      gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                               <  1:144243/149‑1 (MQ=255)
                                                                      gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                               >  7:285697/1‑149 (MQ=255)
                                                                        cccAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTa                                                                             <  4:335352/149‑1 (MQ=255)
                                                                           aaCTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACtt                                                                          >  8:11477/1‑149 (MQ=255)
                                                                                  cGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGaaa                                                                   <  1:104474/149‑1 (MQ=255)
                                                                                      ggTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCg                                                               >  6:328776/1‑149 (MQ=255)
                                                                                       gTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCgg                                                              <  7:11477/149‑1 (MQ=255)
                                                                                       gTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCgg                                                              >  5:180515/1‑149 (MQ=255)
                                                                                       gTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCgg                                                              <  5:328776/149‑1 (MQ=255)
                                                                                          aCTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTg                                                            <  4:1228/148‑1 (MQ=255)
                                                                                                 ggCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACAt                                                     <  6:106960/149‑1 (MQ=255)
                                                                                                      cgGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCAt                                                 <  1:121439/149‑1 (MQ=255)
                                                                                                      cgGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTACCGAAAGCCGGTGAAGACATCCAt                                                 >  3:128433/1‑149 (MQ=255)
                                                                                                       ggTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCAtg                                                <  6:115138/149‑1 (MQ=255)
                                                                                                            cGGGGCGTCACTGTTGTCACCTGTTCGATCCCGTGCGCTTCACTTCCCTGCATGCCTGGCACGTTAATAACGGCGCCGATTTCGAATACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGcc                                           <  4:128433/149‑1 (MQ=255)
                                                                                                                   tCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                    >  3:293793/1‑149 (MQ=255)
                                                                                                                   tCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                    >  5:272927/1‑149 (MQ=255)
                                                                                                                     aCTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTcgc                                  <  8:169283/149‑1 (MQ=255)
                                                                                                                     aCTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTcgc                                  <  4:293793/149‑1 (MQ=255)
                                                                                                                         tgGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGaa                               <  8:154099/148‑1 (MQ=255)
                                                                                                                         tgGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                              <  6:180515/149‑1 (MQ=255)
                                                                                                                         tgGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                              <  6:80924/149‑1 (MQ=255)
                                                                                                                           gTCACCTGTTCGAGCCCGTGCGCTTCTCTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATg                             >  1:183510/1‑148 (MQ=255)
                                                                                                                           gTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATg                             <  2:308193/148‑1 (MQ=255)
                                                                                                                           gTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATg                             >  6:53254/1‑148 (MQ=255)
                                                                                                                              aCCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAg                         >  2:218157/1‑149 (MQ=255)
                                                                                                                               ccTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAgg                        <  6:272927/149‑1 (MQ=255)
                                                                                                                                 tGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGcc                      >  3:29586/1‑149 (MQ=255)
                                                                                                                                     cGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                              <  6:112885/137‑1 (MQ=255)
                                                                                                                                     cGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                              >  5:112885/1‑137 (MQ=255)
                                                                                                                                      gAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTgcg                 >  7:62235/1‑149 (MQ=255)
                                                                                                                                        gCCCGTGCGCTTCACTGCCCTGCATGCCTGGCAAGAGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACCTCCATGCGGCCGTGCCTCGCGAAAGCAAGGCGGTGcgcg               >  4:276121/1‑149 (MQ=255)
                                                                                                                                            gTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGATAGCCGGTGATGACATCCATGCTGACGTGACTcgcg                                 >  6:238126/1‑127 (MQ=255)
                                                                                                                                            gTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTcgcg                                 <  5:238126/127‑1 (MQ=255)
                                                                                                                                            gTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAg           >  6:323265/1‑149 (MQ=255)
                                                                                                                                             cgttcTTCACGTGCCAGGCATGCAGGGCAGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGc          >  5:12879/5‑149 (MQ=38)
                                                                                                                                               ttcTTCACGTGCCAGGCATGCAGGGCAGTGAAGAACGGCTCCGAGTTCGAAGAAGTAGGCAAGTGGAATCGCCCTTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGATGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGt        >  2:27966/3‑149 (MQ=255)
                                                                                                                                               ttcTTCACGTGCCAGGCATGCAGGGCAGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGt        >  6:252412/3‑149 (MQ=255)
                                                                                                                                               cgcTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCg         <  5:182768/148‑1 (MQ=255)
                                                                                                                                                gcTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTg       >  4:337326/1‑149 (MQ=255)
                                                                                                                                                gcTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATACAAGGCCGTGCGCGACAGCGTg       >  5:258652/1‑149 (MQ=255)
                                                                                                                                                  ttCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGaa                               >  3:230158/1‑123 (MQ=255)
                                                                                                                                                  ttCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGaa                               <  4:230158/123‑1 (MQ=255)
                                                                                                                                                  ttCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTggg     >  4:143466/1‑149 (MQ=255)
                                                                                                                                                    cACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGc    <  8:31368/148‑1 (MQ=255)
                                                                                                                                                    cACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATACAAGGCCGTGCGCGACAGCGTGGGc    <  6:258652/148‑1 (MQ=255)
                                                                                                                                                     aCTGCCCTGCATGCCTGGCACGTGAAGTACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCTACAGCGTGGGCct  <  3:276121/149‑1 (MQ=255)
                                                                                                                                                      cTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                    >  1:14312/1‑114 (MQ=255)
                                                                                                                                                      cTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                    <  2:14312/114‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
ATCAACGGCCTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCT  >  NC_002947/391428‑391723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: