Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 390,096 A→C 5.0% M123L (ATG→CTG)  soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947390,0960AC 5.0% 310.7 / 6.2 99M123L (ATG→CTG) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (47/47);  new base C (3/2);  total (50/49)
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.51e-01

Only 100 of 104 total aligned reads displayed.

GATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCACTGC  >  NC_002947/389948‑390242
                                                                                                                                                    |                                                                                                                                                   
gATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCatgat                                                                                                                                                    <  1:2061775/150‑1 (MQ=255)
     gCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATgccgc                                                                                                                                               <  1:1036979/150‑1 (MQ=255)
     gCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTCGGCGGCATCATGCTgccgc                                                                                                                                               <  1:740722/150‑1 (MQ=255)
         cAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGTTTGGCCGAACGCACACAACGAAACTCCAGTCACCTTGGTGAAGGTTGGCGGCAGCATGATCCCGCCggg                                                                                                                                           >  1:2340773/1‑150 (MQ=255)
           gCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCt                                                                                                                                         <  1:1144225/150‑1 (MQ=255)
             gAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTc                                                                                                                                       <  1:513793/150‑1 (MQ=255)
                gCTACCCAGACCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctac                                                                                                                                    <  2:1292880/150‑1 (MQ=255)
                  tACCCAGATCCCCAACGTGCGCGCCACCCAATAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctacta                                                                                                                                  <  2:1929092/150‑1 (MQ=255)
                  tACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctacta                                                                                                                                  >  2:1937260/1‑150 (MQ=255)
                     ccAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCGACCTGGACCAACGGTAGTCCTTACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACaa                                                                                                                               >  1:2100382/1‑150 (MQ=255)
                     ccAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACaa                                                                                                                               >  2:221577/1‑150 (MQ=255)
                       aGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACaaaa                                                                                                                             >  2:493909/1‑150 (MQ=255)
                       aGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACaaaa                                                                                                                             >  2:162927/1‑150 (MQ=255)
                            cccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                        >  2:1726038/1‑150 (MQ=255)
                             ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCa                                                                                                                       >  1:1580697/1‑150 (MQ=255)
                              cAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCAt                                                                                                                      <  2:2260847/150‑1 (MQ=255)
                              cAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCAt                                                                                                                      >  1:1094728/1‑150 (MQ=255)
                                 cGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTa                                                                                                                  <  1:1302572/150‑1 (MQ=255)
                                        gcCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCaaa                                                                                                           <  2:1580697/150‑1 (MQ=255)
                                            cccAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGt                                                                                                       >  1:1132220/1‑150 (MQ=255)
                                             ccAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                        >  1:1637731/1‑133 (MQ=255)
                                             ccAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                        <  2:1637731/133‑1 (MQ=255)
                                              cAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTc                                                                                                     >  2:1766211/1‑150 (MQ=255)
                                                   ggCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGAt                                                                                                <  1:493909/150‑1 (MQ=255)
                                                      aCTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCgg                                                                                                                                            <  2:1535688/104‑1 (MQ=255)
                                                      aCTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCgg                                                                                                                                            >  1:1535688/1‑104 (MQ=255)
                                                      aCTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctac                                                                                                                                    >  2:1333771/1‑112 (MQ=255)
                                                      aCTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctac                                                                                                                                    <  1:1333771/112‑1 (MQ=255)
                                                        tGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCt                                                                                                    >  2:1039037/1‑141 (MQ=255)
                                                        tGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCt                                                                                                    <  1:1039037/141‑1 (MQ=255)
                                                         gTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTa                                                                                          <  2:1094728/150‑1 (MQ=255)
                                                           aCGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACg                                                                                        <  2:2100382/150‑1 (MQ=255)
                                                            cgcgGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTACCGCATGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACga                                                                                       <  1:2068284/150‑1 (MQ=255)
                                                              tggccattGTCGCCACCAGCACCAACGGTTGGCCGATCGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACgaga                                                                                     >  2:2210519/7‑150 (MQ=255)
                                                               gggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGACTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGaa                                                                                    >  1:793842/1‑150 (MQ=255)
                                                               gggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGAAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGaa                                                                                    <  1:1747206/150‑1 (MQ=255)
                                                                 ccattGTCGCCACCAGCACCAACGGTTGGCCGATCGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGt                                                                                  <  1:2210519/147‑1 (MQ=255)
                                                                  ccTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTa                                                                                 >  2:2044036/1‑150 (MQ=255)
                                                                    ttGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACa                                                                               >  2:1380416/1‑150 (MQ=255)
                                                                       tCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATcc                                                                            >  1:2192503/1‑150 (MQ=255)
                                                                         gCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTc                                                                                                                                       <  2:1109545/90‑1 (MQ=255)
                                                                         gCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTc                                                                                                                                       >  1:1109545/1‑90 (MQ=255)
                                                                          ccaccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCAGCTCCCAGGCACTGCGCAAGGCAGCTGTGACGCGAAATCGTTCTGGATGACTTACGAGAATTACATCCGTa                                                                         >  1:284609/1‑150 (MQ=255)
                                                                           caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCtt                                                                                                                                        <  2:2250826/87‑1 (MQ=255)
                                                                           caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCtt                                                                                                                                        >  1:2250826/1‑87 (MQ=255)
                                                                           caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTaa                                                                        >  1:1048130/1‑150 (MQ=255)
                                                                            accaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTaaa                                                                       >  2:2274967/1‑150 (MQ=255)
                                                                               aGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgcc                                                                    <  2:1223076/150‑1 (MQ=255)
                                                                                 cACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTcta                                                                                                                                     <  2:1038350/84‑1 (MQ=255)
                                                                                 cACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTcta                                                                                                                                     >  1:1038350/1‑84 (MQ=255)
                                                                                  aCCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAGGTACATCCGTAAAgcagcc                                                                 >  1:1662268/1‑150 (MQ=255)
                                                                                     aaCGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGc                                                              >  2:429877/1‑150 (MQ=255)
                                                                                              ccGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGt                                                     <  2:832148/150‑1 (MQ=255)
                                                                                                gAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgc                                                   <  2:1690433/150‑1 (MQ=255)
                                                                                                   cGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGcgcc                                                >  1:362765/1‑150 (MQ=255)
                                                                                                         caacGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCAt                                                                                                                      <  1:1993495/75‑1 (MQ=255)
                                                                                                            cgacgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGaa                                       <  1:429877/150‑1 (MQ=255)
                                                                                                             gacgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACgag                                                                                      <  1:731419/102‑1 (MQ=255)
                                                                                                             gacgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACgag                                                                                      >  2:731419/1‑102 (MQ=255)
                                                                                                             gacgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGt                                                     >  1:2338187/1‑135 (MQ=255)
                                                                                                             gacgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGt                                                     <  2:2338187/135‑1 (MQ=255)
                                                                                                             gacgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgc                                                   >  2:1585451/1‑137 (MQ=255)
                                                                                                             gacgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgc                                                   <  1:1585451/137‑1 (MQ=255)
                                                                                                             gacgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAAc                                      <  2:362765/150‑1 (MQ=255)
                                                                                                                 tCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATc                                  >  1:1574450/1‑150 (MQ=255)
                                                                                                                 tCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATc                                  <  2:1259556/150‑1 (MQ=255)
                                                                                                                  cATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATcc                                 <  2:284609/150‑1 (MQ=255)
                                                                                                                   aTGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgc                                                                  <  2:1246165/116‑1 (MQ=255)
                                                                                                                   aTGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgc                                                                  >  1:1246165/1‑116 (MQ=255)
                                                                                                                   aTGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCGGCCGGCCGGGGCCGTGCGCCGCTGCAGAACGATCCt                                >  2:760148/1‑150 (MQ=255)
                                                                                                                     gggCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGa                              <  1:344426/150‑1 (MQ=255)
                                                                                                                      ggCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCNGTAAAGNCGCCGGCCTGGGNCGTGCGCCGCTGCAGAACGATCCTGAc                             <  1:1766211/150‑1 (MQ=255)
                                                                                                                       gCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgcc                                                                 <  2:1651201/113‑1 (MQ=255)
                                                                                                                       gCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgcc                                                                 >  1:1651201/1‑113 (MQ=255)
                                                                                                                       gCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACa                            <  2:656212/150‑1 (MQ=255)
                                                                                                                          tCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAATCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCt                         >  2:2098102/1‑150 (MQ=255)
                                                                                                                          tCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGa                                    >  1:390497/1‑139 (MQ=255)
                                                                                                                          tCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGa                                    <  2:390497/139‑1 (MQ=255)
                                                                                                                          tCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCt                         <  1:1726038/150‑1 (MQ=255)
                                                                                                                          tCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCt                         <  1:950695/150‑1 (MQ=255)
                                                                                                                             tCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgc                                                                  >  2:290143/1‑106 (MQ=255)
                                                                                                                             tCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgc                                                                  <  1:290143/106‑1 (MQ=255)
                                                                                                                             tCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACg                      >  1:960870/1‑150 (MQ=255)
                                                                                                                              cGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGa                     <  1:1380416/150‑1 (MQ=255)
                                                                                                                               ggCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGAc                    <  1:162927/150‑1 (MQ=255)
                                                                                                                                gCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACa                                                                               >  1:466475/1‑90 (MQ=255)
                                                                                                                                gCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACa                                                                               <  2:466475/90‑1 (MQ=255)
                                                                                                                                    ggTTGGCGGCAGCATGATGCCGCTGGGCTTCTACTACAAAACCTTCAt                                                                                                                      >  1:1743470/1‑48 (MQ=255)
                                                                                                                                    ggTTGGCGGCAGCATGATGCCGCTGGGCTTCTACTACAAAACCTTCAt                                                                                                                      <  2:1743470/48‑1 (MQ=255)
                                                                                                                                      ccggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGa             >  1:227713/3‑150 (MQ=255)
                                                                                                                                       cggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAg                                                                                   <  1:2343286/78‑1 (MQ=39)
                                                                                                                                       cggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAg                                                                                   >  2:2343286/2‑79 (MQ=39)
                                                                                                                                         gcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAAcc          >  1:953952/1‑150 (MQ=255)
                                                                                                                                         gcggcAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAAcc          >  2:331972/1‑150 (MQ=255)
                                                                                                                                          cggcAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCg         <  1:2044036/150‑1 (MQ=255)
                                                                                                                                           ggcAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCgg        >  2:371361/1‑150 (MQ=255)
                                                                                                                                                catcatGCCGCCGTGCTTCTACTACAAAACCTTCATGTACCCACAATCGTTCTGGATGACTTACGAGAATTACTTCCGTAAAGCCGCCGGCCCGGGCGGTGCGCCCCTGCAGATCGATCCTGACAGCTACGACTACGCGCACCGGCACTg   <  2:799975/150‑1 (MQ=255)
                                                                                                                                                 atgatgCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCACTGc  >  1:1710478/1‑150 (MQ=255)
                                                                                                                                                   gatgCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGaa                                                                                    >  2:499513/1‑66 (MQ=255)
                                                                                                                                                   gatgCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGaa                                                                                    <  1:499513/66‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCACTGC  >  NC_002947/389948‑390242

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