Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 87,956 G→C 5.4% R454R (CGC→CGG betC ← choline‑sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294787,9560GC 5.4% 358.9 / 8.7 111R454R (CGC→CGGbetCcholine‑sulfatase
Reads supporting (aligned to +/- strand):  ref base G (57/47);  new base C (3/3);  total (61/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 120 total aligned reads displayed.

AGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAG  >  NC_002947/87807‑88105
                                                                                                                                                     |                                                                                                                                                     
aGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcgg                                                                                                                                                       >  2:1905609/1‑150 (MQ=255)
  gggggCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcggcg                                                                                                                                                     <  1:1094206/150‑1 (MQ=255)
  gggggCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcggcg                                                                                                                                                     >  2:516779/1‑150 (MQ=255)
    gggcgggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCCTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCTATGGCCAGCGCTTCGGCCACCAGCCggcggcggc                                                                                                                                                   <  1:817833/150‑1 (MQ=255)
         ggcTCTGGATAACGTGCCTTGCGCTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                              <  1:936536/150‑1 (MQ=255)
          gcTGTGGATAACGTGCCTTGCGTGCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggc                                                                                                                                             >  2:1188445/1‑150 (MQ=255)
            tgtgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggctg                                                                                                                                           >  1:1713243/1‑150 (MQ=255)
                  tAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGt                                                                                                                                     >  1:2256838/1‑150 (MQ=255)
                   aaCGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGTGGCGGAGCTGGCTGGCCAGTa                                                                                                                                    >  2:2086247/1‑150 (MQ=255)
                      gTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcggcgg                                                                                                                                                    <  2:711866/131‑1 (MQ=255)
                      gTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcggcgg                                                                                                                                                    >  1:711866/1‑131 (MQ=255)
                      gTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACct                                                                                                                                 >  2:1573829/1‑150 (MQ=255)
                      gTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACct                                                                                                                                 <  2:1835474/150‑1 (MQ=255)
                         ccTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACctgct                                                                                                                              >  1:485308/1‑150 (MQ=255)
                          cTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACctgttg                                                                                                                             <  2:1862568/150‑1 (MQ=255)
                            tGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGgc                                                                                                                           <  2:107294/150‑1 (MQ=255)
                            tGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGACGGAGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGgc                                                                                                                           >  1:626572/1‑150 (MQ=255)
                                       tCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTg                                                                                                                >  2:1573172/1‑150 (MQ=255)
                                         accGAGATCGCTGTGGGTGCGCATGTACTGTTGGCTGGCGACCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGAt                                                                                                              >  2:2439916/3‑150 (MQ=255)
                                         aTCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGAt                                                                                                              >  1:1555159/1‑150 (MQ=255)
                                              gaTCGATGTGGTCGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTccc                                                                                                         >  1:1795918/1‑150 (MQ=255)
                                               aTCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCa                                                                                                        <  1:1905609/150‑1 (MQ=255)
                                               aTCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCa                                                                                                        >  1:1267711/1‑150 (MQ=255)
                                                        gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTg                                                                                               >  1:502439/1‑150 (MQ=255)
                                                        gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTg                                                                                               >  2:1475261/1‑150 (MQ=255)
                                                         ttGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTgc                                                                                              <  1:1067282/150‑1 (MQ=255)
                                                         ttGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTgc                                                                                              <  2:1713243/150‑1 (MQ=255)
                                                             gcATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGTTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCtt                                                                                          <  1:356998/150‑1 (MQ=255)
                                                               aTGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCa                                                                                        >  1:1110571/1‑150 (MQ=255)
                                                                 gTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGTCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTg                                                                                                                   <  1:1432397/121‑1 (MQ=255)
                                                                 gTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGTCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTg                                                                                                                   >  2:1432397/1‑121 (MQ=255)
                                                                    cTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGAc                                                                                   <  1:1573172/150‑1 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTTGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                  >  2:2418479/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCGCAGCGCTGTTGCGCTTCATCGACa                                                                                  >  2:244730/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                      >  1:2214053/1‑146 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                      <  2:2214053/146‑1 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                  >  2:1136312/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                  >  2:1627730/1‑150 (MQ=255)
                                                                       ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACaaa                                                                                <  1:1136312/150‑1 (MQ=255)
                                                                       ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACaaa                                                                                <  1:1136867/150‑1 (MQ=255)
                                                                        tggctggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTGCAGCTTGCCGATGGCCAGCGCTGCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGCTGTCCCAGCGCTGTTGCGCTTCATCGACAAAt                                                                               <  2:146363/150‑1 (MQ=255)
                                                                           ctggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGcc                                                                            <  1:516779/150‑1 (MQ=255)
                                                                            tggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCt                                                                           >  2:2240199/1‑150 (MQ=255)
                                                                              gcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgc                                                                         >  2:1566556/1‑150 (MQ=255)
                                                                                gTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGTCCGCGCGTGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcag                                                                       >  1:159702/1‑150 (MQ=255)
                                                                                gTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcag                                                                       >  2:673094/1‑150 (MQ=255)
                                                                                gTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcag                                                                       >  1:1536240/1‑150 (MQ=255)
                                                                                   caccaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagcag                                                                    >  2:765150/1‑150 (MQ=255)
                                                                                     ccaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggctggc                                                                                                                                         <  1:281557/79‑1 (MQ=255)
                                                                                     ccaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggctggc                                                                                                                                         >  2:281557/1‑79 (MQ=255)
                                                                                      caGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCAc                                                                 >  1:2265860/1‑150 (MQ=255)
                                                                                       aGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGGAGCAcc                                                                >  1:1185674/1‑150 (MQ=255)
                                                                                       aGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCAAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCAcc                                                                >  1:155379/1‑150 (MQ=255)
                                                                                         tggttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtg                                                              <  1:572416/150‑1 (MQ=255)
                                                                                          ggttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                             <  1:1566556/150‑1 (MQ=255)
                                                                                          ggttggtCGTCACAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCGGTCGCCACCGGATCTTCCAGCCGATGTCCCGCTGGCGCAGACTGGCGCTGTAACAGCGATGAAGCGCTTCATTGTCAAATGCCTGCAGCAGCACCtgg                                                             <  1:1125695/150‑1 (MQ=255)
                                                                                           gttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggt                                                            <  1:2439916/150‑1 (MQ=255)
                                                                                           gttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggt                                                            >  1:1139166/1‑150 (MQ=255)
                                                                                                  gtCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCgg                                                     <  2:438684/150‑1 (MQ=255)
                                                                                                    cccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAACGGTATTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTACGGGc                                                   >  1:522379/1‑150 (MQ=255)
                                                                                                    cccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGACGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGc                                                   >  1:224221/1‑150 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgc                                                                         <  2:1933419/127‑1 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgc                                                                         >  1:1933419/1‑127 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCt                                                  >  2:4804/1‑150 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCt                                                  >  1:1180224/1‑150 (MQ=255)
                                                                                                      cAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTg                                                 >  2:2183455/1‑150 (MQ=255)
                                                                                                       aGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGc                                                >  1:1083918/1‑150 (MQ=255)
                                                                                                       aGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGc                                                <  2:1180224/150‑1 (MQ=255)
                                                                                                             tCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                          <  2:1439330/150‑1 (MQ=255)
                                                                                                              cAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                         <  1:1670467/150‑1 (MQ=255)
                                                                                                               aGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCTGGCTGCCGGTGAgg                                        <  1:1627730/150‑1 (MQ=255)
                                                                                                               aGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAgg                                        >  1:1648997/1‑150 (MQ=255)
                                                                                                                gCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCACAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                      >  2:576182/1‑103 (MQ=255)
                                                                                                                gCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCACAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                      <  1:576182/103‑1 (MQ=255)
                                                                                                                    gCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTtctc                                   <  2:1536240/150‑1 (MQ=255)
                                                                                                                       gATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTgg                                                                                                                            >  2:1132326/1‑58 (MQ=255)
                                                                                                                       gATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGt                                            <  2:2374480/138‑1 (MQ=255)
                                                                                                                         tGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCt                              >  1:296804/1‑150 (MQ=255)
                                                                                                                           gCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCg                            >  1:154441/1‑150 (MQ=255)
                                                                                                                               gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTgcg                                                                                             >  1:239925/1‑81 (MQ=255)
                                                                                                                               gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTgcg                                                                                             <  2:239925/81‑1 (MQ=255)
                                                                                                                                cgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgg                       <  2:1648997/150‑1 (MQ=255)
                                                                                                                                cgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgg                       <  2:1887851/150‑1 (MQ=255)
                                                                                                                                cgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgg                       <  1:1987659/150‑1 (MQ=255)
                                                                                                                                 gcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGACGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCggg                      <  2:224221/150‑1 (MQ=255)
                                                                                                                                  cTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCACGAGCggga                     <  1:1568994/150‑2 (MQ=255)
                                                                                                                                       gCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgc                  <  2:2497/148‑1 (MQ=255)
                                                                                                                                       gCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgc                  >  1:2497/1‑148 (MQ=255)
                                                                                                                                       gCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCgg                <  1:2183455/150‑1 (MQ=255)
                                                                                                                                         ggccagccagCGCCGCCGCCGGCTGGTCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCt              >  1:1954229/3‑150 (MQ=255)
                                                                                                                                         caccaGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgg                       <  2:1544712/141‑1 (MQ=255)
                                                                                                                                          gccagccagCGCCGCCGCCGGCTGGTCAGTACCTGCTGGCGCAGGCTGGTGATAAACCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTc             >  1:1356618/2‑150 (MQ=255)
                                                                                                                                          gccagccagCGCCGCCGCCGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                         <  2:672546/137‑1 (MQ=255)
                                                                                                                                          gccagccagCGCCGCCGCCGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                         >  1:672546/2‑138 (MQ=255)
                                                                                                                                          accaGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTc             >  2:484383/1‑150 (MQ=255)
                                                                                                                                            caGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAg           >  2:752475/1‑150 (MQ=255)
                                                                                                                                               ccGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTc        <  2:1187391/150‑1 (MQ=255)
                                                                                                                                                     gcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGGGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATagag  >  2:376132/1‑150 (MQ=255)
                                                                                                                                                     gcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGGGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATagag  >  2:375993/1‑150 (MQ=255)
                                                                                                                                                     gcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATagag  <  2:343506/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAG  >  NC_002947/87807‑88105

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: