Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 391,595 Δ1 bp 6.1% coding (1869/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,5950C. 6.1% 264.4 / 10.1 82coding (1869/3015 nt)soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (35/41);  new base . (3/2);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

TTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGACCCTGG  >  NC_002947/391454‑391743
                                                                                                                                              |                                                                                                                                                    
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:129081/1‑149 (MQ=255)
             aCCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGtt                                                                                                                                  >  3:124082/1‑149 (MQ=255)
               cATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACg                                                                                                                                            >  4:289859/1‑137 (MQ=255)
               cATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACg                                                                                                                                            <  3:289859/137‑1 (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                                                                                                   >  7:285697/1‑149 (MQ=255)
                                            gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                                   <  4:124082/149‑1 (MQ=255)
                                            gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                                   <  4:129081/149‑1 (MQ=255)
                                            gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                                   <  1:144243/149‑1 (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                                                                                  <  5:328776/149‑1 (MQ=255)
                                                             gTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCgg                                                                                  <  7:11477/149‑1 (MQ=255)
                                                             gTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCgg                                                                                  >  5:180515/1‑149 (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                                                      <  4:293793/149‑1 (MQ=255)
                                                                                         aCTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTcgc                                                      <  8:169283/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                                                 >  6:53254/1‑148 (MQ=255)
                                                                                               gTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATg                                                 <  2:308193/148‑1 (MQ=255)
                                                                                                  aCCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAg                                             >  2:218157/1‑149 (MQ=255)
                                                                                                   ccTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAgg                                            <  6:272927/149‑1 (MQ=255)
                                                                                                     tGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGcc                                          >  3:29586/1‑149 (MQ=255)
                                                                                                         cGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                                                  >  5:112885/1‑137 (MQ=255)
                                                                                                         cGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                                                  <  6:112885/137‑1 (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                                                        <  2:14312/114‑1 (MQ=255)
                                                                                                                          cTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                                        >  1:14312/1‑114 (MQ=255)
                                                                                                                              cccTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCctgctg                  <  2:183510/149‑1 (MQ=255)
                                                                                                                              cccTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCctgct                   <  1:218157/148‑1 (MQ=255)
                                                                                                                                 tGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGAc               >  6:50804/1‑149 (MQ=255)
                                                                                                                                   cATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGAc               <  4:29586/147‑1 (MQ=255)
                                                                                                                                   cATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGc             <  4:72823/149‑1 (MQ=255)
                                                                                                                                    aTGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGcc            <  3:13012/149‑1 (MQ=255)
                                                                                                                                      gCCTGGCACATGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTc          >  4:326690/1‑149 (MQ=255)
                                                                                                                                         tGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGAc       >  3:10844/1‑149 (MQ=255)
                                                                                                                                             aCGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGa        <  5:14687/144‑1 (MQ=255)
                                                                                                                                             aCGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGa        >  6:14687/1‑144 (MQ=255)
                                                                                                                                             aCGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGACCCTg   <  8:285697/149‑1 (MQ=255)
                                                                                                                                              cGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGACCCTgg  <  8:2771/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGACCCTGG  >  NC_002947/391454‑391743

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: