Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 87,956 G→C 7.1% 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 7.1% 273.0 / 10.2 85R454R (CGC→CGGbetCcholine‑sulfatase
Reads supporting (aligned to +/- strand):  ref base G (45/34);  new base C (4/2);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

AGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGA  >  NC_002947/87807‑88104
                                                                                                                                                     |                                                                                                                                                       
aGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcgg                                                                                                                                                         >  2:682178/1‑150 (MQ=255)
     ggcgggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGgcggcggcg                                                                                                                                                    >  2:996156/1‑150 (MQ=255)
        gggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCAGGAGGCGGCgctg                                                                                                                                                 >  2:101823/1‑150 (MQ=255)
         ggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                                <  1:867548/150‑1 (MQ=255)
            tgtgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                                <  2:430808/147‑1 (MQ=255)
            tgtgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                                >  1:430808/1‑147 (MQ=255)
                gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCCGACTGCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCa                                                                                                                                      >  1:1071392/1‑150 (MQ=255)
                  tAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGt                                                                                                                                    >  2:388039/1‑150 (MQ=255)
                          cTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACctgctg                                                                                                                            >  2:411309/1‑150 (MQ=255)
                           ttGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTgg                                                                                                                           >  1:310358/1‑150 (MQ=255)
                            tGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTATTGGCTGGCATCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGgc                                                                                                                          >  1:1038711/1‑150 (MQ=255)
                                     ggTCATCGAGATCGAGGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                                 <  2:249702/150‑1 (MQ=255)
                                         aTCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACctgc                                                                                                                              <  2:679814/133‑1 (MQ=255)
                                         aTCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACctgc                                                                                                                              >  1:679814/1‑133 (MQ=255)
                                          tCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATg                                                                                                            >  1:814478/1‑150 (MQ=255)
                                          tCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATg                                                                                                            >  2:411766/1‑150 (MQ=255)
                                           cGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCACTACCTGCTGGCGCAGGCTGGTGATGt                                                                                                           <  2:15883/150‑1 (MQ=255)
                                                  gATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAgcg                                                                                                    >  2:1281674/1‑150 (MQ=255)
                                                       ggTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                               <  1:1320878/150‑1 (MQ=255)
                                                       ggTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                               <  1:1209866/150‑1 (MQ=255)
                                                       ggTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                               <  1:411766/150‑1 (MQ=255)
                                                        gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                                 <  1:884128/131‑1 (MQ=255)
                                                        gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                                 >  2:884128/1‑131 (MQ=255)
                                                             gcATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCa                                                                                                       >  2:1154447/1‑136 (MQ=255)
                                                             gcATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCa                                                                                                       <  1:1154447/136‑1 (MQ=255)
                                                              cATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTc                                                                                        <  1:682178/150‑1 (MQ=255)
                                                              cATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTc                                                                                        <  2:846704/150‑1 (MQ=255)
                                                                  tACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggct‑‑‑gg                                                                                                                                         <  1:551449/97‑1 (MQ=255)
                                                                  tACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggct‑‑‑gg                                                                                                                                         >  2:551449/1‑97 (MQ=255)
                                                                   aCTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCGTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGa                                                                                   >  1:608516/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTTAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                 >  1:625750/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                 >  1:898696/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                 >  2:216860/1‑150 (MQ=255)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggct‑‑‑gg                                                                                                                                         >  2:806926/1‑94 (MQ=38)
                                                                     tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggct‑‑‑gg                                                                                                                                         <  1:806926/94‑1 (MQ=38)
                                                                       ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggct‑‑‑ggc                                                                                                                                        >  2:761152/1‑93 (MQ=255)
                                                                       ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggct‑‑‑ggc                                                                                                                                        <  1:761152/93‑1 (MQ=255)
                                                                            tggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCt                                                                          >  1:138474/1‑150 (MQ=255)
                                                                            tggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCt                                                                          <  1:996156/150‑1 (MQ=255)
                                                                            aggcGTGCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCt                                                                          >  2:1065482/2‑150 (MQ=255)
                                                                                    accaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagcagc                                                                  >  2:405431/1‑150 (MQ=255)
                                                                                     ccaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagcagca                                                                 >  1:1220419/1‑150 (MQ=255)
                                                                                      caGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCAc                                                                >  2:1092827/1‑150 (MQ=255)
                                                                                        gTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGa                                                                                                              >  1:17521/1‑102 (MQ=255)
                                                                                        gTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGa                                                                                                              <  2:17521/102‑1 (MQ=255)
                                                                                         tggttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtg                                                             <  2:608516/150‑1 (MQ=255)
                                                                                           gttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggt                                                           >  2:99451/1‑150 (MQ=255)
                                                                                            ttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggtg                                                          <  2:138474/150‑1 (MQ=255)
                                                                                                  gtCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCgg                                                    <  1:216860/150‑1 (MQ=255)
                                                                                                   tCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTAATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   >  1:946141/1‑150 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTa                                                                                                                                   >  1:1388833/1‑68 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTa                                                                                                                                   <  2:1388833/68‑1 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                      <  2:120161/113‑1 (MQ=255)
                                                                                                     ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                      >  1:120161/1‑113 (MQ=255)
                                                                                                      cAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   >  1:882845/1‑147 (MQ=255)
                                                                                                      cAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   <  2:882845/147‑1 (MQ=255)
                                                                                                       aGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCGGTTGCGCTGCATCGACAAAAGCCTGCAGCAGCACCTGGTGGTCCGGGCTgt                                               >  2:654163/1‑149 (MQ=255)
                                                                                                          ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCgg                                            >  2:355166/1‑150 (MQ=255)
                                                                                                            ttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTg                                          >  1:45853/1‑150 (MQ=255)
                                                                                                             tCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                         <  2:1038711/150‑1 (MQ=255)
                                                                                                             tCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                         <  1:411309/150‑1 (MQ=255)
                                                                                                              cAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCAAGCGCTATTGCGCTTCATCGACAAAGACCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                        <  1:388039/150‑1 (MQ=255)
                                                                                                               aGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAgg                                       >  2:1348175/1‑150 (MQ=255)
                                                                                                                     ccGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTcc                                 <  2:1241011/150‑1 (MQ=255)
                                                                                                                      cGATGGCCAGCGCTTCGAAAACAAGCCGGCGGCGGCGCTGTCT‑‑‑GTCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCCGCGATGTTACGCTTCATCGACAAATGTCGCGATCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTccc                                >  2:1020983/1‑150 (MQ=255)
                                                                                                                       gATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCg                               >  1:707732/1‑150 (MQ=255)
                                                                                                                        aTGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGc                              <  2:814478/150‑1 (MQ=255)
                                                                                                                            ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                            <  2:894137/116‑1 (MQ=255)
                                                                                                                            ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                            >  1:894137/1‑116 (MQ=255)
                                                                                                                            ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGt                          >  2:25721/1‑150 (MQ=255)
                                                                                                                             cAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTg                         >  2:1130791/1‑150 (MQ=255)
                                                                                                                              aGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                        <  2:898696/150‑1 (MQ=255)
                                                                                                                              aGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                        <  1:1224803/150‑1 (MQ=255)
                                                                                                                               gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTtc                                    >  1:755243/1‑137 (MQ=255)
                                                                                                                               gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTtc                                    <  2:755243/137‑1 (MQ=255)
                                                                                                                                       gCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCgg               <  2:1071392/150‑1 (MQ=255)
                                                                                                                                         ggccagccagCGCCGCCGCCGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCt             >  2:276604/3‑150 (MQ=255)
                                                                                                                                         caccaGCCGGCGGCGGCGCTGGCTGGCGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgcg                >  2:1258255/1‑150 (MQ=255)
                                                                                                                                         caccaGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCt             >  1:53860/1‑150 (MQ=255)
                                                                                                                                         caccaGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCACGCTCGTGCGGGGCGCGGCt             >  2:1327261/1‑150 (MQ=255)
                                                                                                                                          accaGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTc            >  2:915190/1‑150 (MQ=255)
                                                                                                                                            cagccagcGCCGCCGCCGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAg          >  1:130311/1‑150 (MQ=255)
                                                                                                                                             aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   <  2:1263987/108‑1 (MQ=255)
                                                                                                                                             aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   >  1:1263987/1‑108 (MQ=255)
                                                                                                                                             aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAgg         <  2:583917/150‑1 (MQ=255)
                                                                                                                                             aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAgg         <  1:1281674/150‑1 (MQ=255)
                                                                                                                                             aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAgg         >  1:1261863/1‑150 (MQ=255)
                                                                                                                                                cggcggcggCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCa      >  1:1150388/1‑150 (MQ=255)
                                                                                                                                                  gcggcggcgCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATa    <  1:654163/150‑1 (MQ=255)
                                                                                                                                                    ggcggcgCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATaga  <  1:355166/150‑1 (MQ=255)
                                                                                                                                                    ggcggcgCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATaga  <  2:310358/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AGGGGGGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGA  >  NC_002947/87807‑88104

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: