Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 391,581 C→A 6.3% L619M (CTG→ATG) ‡ soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947391,5810CA 6.3% 256.5 / 5.9 79L619M (CTG→ATG) ‡soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base A (3/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

CTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTG  >  NC_002947/391437‑391727
                                                                                                                                                   |                                                                                                                                                  
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                                                                                               >  6:182768/1‑148 (MQ=255)
                                                  caccaTGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                               >  3:335352/1‑148 (MQ=255)
                                                  caccaTGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                               >  3:1228/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                                 <  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                                  >  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)
                                                                                                                                                   |                                                                                                                                                  
CTGGCCATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTACTTTCGGCGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACTGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTG  >  NC_002947/391437‑391727

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: