Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 391,582 T→G 6.2% L619R (CTG→CGG) ‡ soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947391,5820TG 6.2% 260.4 / 7.1 81L619R (CTG→CGG) ‡soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (36/40);  new base G (3/2);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

CTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGAC  >  NC_002947/391437‑391730
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                               >  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:335352/1‑148 (MQ=255)
                                                  caccaTGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                                  >  3:1228/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:124082/149‑1 (MQ=255)
                                                             gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                      <  4:129081/149‑1 (MQ=255)
                                                             gCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                      >  7:285697/1‑149 (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                                                                     <  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                                         <  4:293793/149‑1 (MQ=255)
                                                                                                            aCTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTcgc                                         <  8:169283/149‑1 (MQ=255)
                                                                                                                tgGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGaa                                      <  8:154099/148‑1 (MQ=255)
                                                                                                                tgGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                                     <  6:80924/149‑1 (MQ=255)
                                                                                                                tgGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                                     <  6:180515/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                                     <  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                                      <  4:230158/123‑1 (MQ=255)
                                                                                                                                         ttCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGaa                                      >  3:230158/1‑123 (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)
                                                                                                                                                    |                                                                                                                                                    
CTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGAC  >  NC_002947/391437‑391730

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: