Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,984,228 C→A 7.7% intergenic (+35/+124) aslB → / ← aslA putative anaerobic sulfatase maturation enzyme AslB/arylsulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,984,2280CA 7.7% 506.9 / 20.9 155intergenic (+35/+124)aslB/aslAputative anaerobic sulfatase maturation enzyme AslB/arylsulfatase
Reads supporting (aligned to +/- strand):  ref base C (69/74);  new base A (4/8);  total (73/82)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

Only 100 of 166 total aligned reads displayed.

CAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTG  >  NC_000913/3984080‑3984377
                                                                                                                                                    |                                                                                                                                                     
 aGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcgg                                                                                                                                                     >  2:1210866/1‑150 (MQ=255)
 aGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcgg                                                                                                                                                     >  1:766599/1‑150 (MQ=255)
 aGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcgg                                                                                                                                                     >  1:552408/1‑150 (MQ=255)
    gTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATAAGGCGt                                                                                                                                                  >  1:1085912/1‑150 (MQ=255)
      tATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATAAGGCGttt                                                                                                                                                <  2:1223513/150‑3 (MQ=255)
           ccGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGcc                                                                                                                                           <  1:969576/150‑1 (MQ=255)
           ccGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGcc                                                                                                                                           <  1:107234/150‑1 (MQ=255)
                  cTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                    >  1:834420/1‑150 (MQ=255)
                    aCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                  <  1:1020072/150‑1 (MQ=255)
                    aCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                  <  1:439976/150‑1 (MQ=255)
                    aCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGGCGCCCGGCGAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCTCCTGCGGTTGATTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                  <  1:181091/150‑1 (MQ=255)
                     ccGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                 >  1:1227873/1‑150 (MQ=255)
                      cGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcgg                                                                                                                                                     <  2:331971/129‑1 (MQ=255)
                       gCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                               >  2:40238/1‑150 (MQ=255)
                       gCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                               >  1:223003/1‑150 (MQ=255)
                       gCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATAAGGCGTTTACGCCGCATCCGGcaat                                                                                                                               <  2:137533/150‑4 (MQ=255)
                         cATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                                                                             <  1:643027/150‑1 (MQ=255)
                         cATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                                                                             >  1:859390/1‑150 (MQ=255)
                          cagcaaatgcgcatgcataatgtcgctggccgggcgaccgtgcgccagcGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAt                                                                                                                            >  1:965262‑M2/50‑150 (MQ=255)
                           tatCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                              <  2:1025437/147‑1 (MQ=255)
                             tCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGAt                                                                                                                         >  1:259924/1‑150 (MQ=255)
                             tCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGAt                                                                                                                         >  1:655180/1‑150 (MQ=255)
                               ttAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCg                                                                                                                       >  1:351984/1‑150 (MQ=255)
                                 aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTg                                                                                                                     <  2:286587/150‑1 (MQ=255)
                                  aaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGc                                                                                                                    >  2:1253342/1‑150 (MQ=255)
                                  aaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGc                                                                                                                    <  2:216560/150‑1 (MQ=255)
                                    gCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaa                                                                                                                  >  1:1127058/1‑150 (MQ=255)
                                      aaTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAAt                                                                                                                <  1:741174/150‑1 (MQ=255)
                                      aaTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAAt                                                                                                                <  1:1277616/150‑1 (MQ=255)
                                      aaTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAAt                                                                                                                >  2:761000/1‑150 (MQ=255)
                                      aaTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAAt                                                                                                                <  1:435162/150‑1 (MQ=255)
                                       aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAAtt                                                                                                               <  2:1150892/150‑1 (MQ=255)
                                        tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTc                                                                                                              <  1:277858/150‑1 (MQ=255)
                                           cTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAAt                                                                                                           <  1:935576/150‑1 (MQ=255)
                                           cTGATTTGCTGGAGCACTGTCTCCCTGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCATGCGGTTGATTGCCGTATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAAt                                                                                                           <  1:1072494/150‑1 (MQ=38)
                                              aTTTGCTGGCGCACGGTCGCCCGGCCAGCTACATTATGCAGGCGCATTTGCTGGTGTGGAGTAAGTAGAAAGCGGCGGCCTGCTGATGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATaaa                                                                                                        <  2:68378/150‑1 (MQ=37)
                                                  GCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGc                                                                                                    <  2:1086426‑M1/150‑1 (MQ=255)
                                                     GGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                                                                                 <  2:929985‑M1/150‑1 (MQ=255)
                                                         cACGGGCACCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGGTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTc                                                                                             <  2:344130/150‑1 (MQ=255)
                                                           cGGTCGCCCGTCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTg                                                                                           <  2:243788/150‑1 (MQ=255)
                                                           cGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTg                                                                                           <  2:851600/150‑1 (MQ=255)
                                                           cGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTg                                                                                           >  1:1340591/1‑150 (MQ=255)
                                                           cGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTg                                                                                           >  1:243788/1‑150 (MQ=255)
                                                            ggTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGt                                                                                          <  2:462910/150‑1 (MQ=255)
                                                               cGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAgg                                                                                       <  1:722265/150‑1 (MQ=255)
                                                               cGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAgg                                                                                       >  1:268599/1‑150 (MQ=255)
                                                                  ccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTg                                                                                    <  1:234257/150‑1 (MQ=255)
                                                                     gccagcGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGAt                                                                                 <  1:536204‑M2/144‑1 (MQ=255)
                                                                      ccagcGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATa                                                                                <  2:1127058‑M2/145‑1 (MQ=255)
                                                                      ccagcGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATa                                                                                <  2:141788‑M2/145‑1 (MQ=255)
                                                                       cagcGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATaa                                                                               >  1:610325‑M2/5‑150 (MQ=255)
                                                                         gcGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGa                                                                             >  2:896294‑M2/3‑150 (MQ=255)
                                                                          cGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGAt                                                                            >  1:344607‑M2/2‑150 (MQ=255)
                                                                             cATTATGCATGCGCATTTGCGGGGGGTGAGTAAGTAGACATCGGCGGCCGCCTGCGGTGGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGCTCGTGCAAATTCAATAAATTTCAGCGTTCTGTAGGCGGGATAAGATGCg                                                                         <  2:712218/150‑1 (MQ=37)
                                                                               ttATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTc                                                                       <  2:321167/150‑1 (MQ=255)
                                                                               ttATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTc                                                                       <  2:859390/150‑1 (MQ=255)
                                                                                 attcatgcGCATTTGCTTGTTGTTAGTAAGTAGAATTCTGCGTCCGCCTGCGTTTGATTGCCTGATGCGTCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATATATTGCAGCGTTCTTTAGGCTGTATAAGATGCGTCAg                                                                     <  2:866711/147‑1 (MQ=25)
                                                                                    catgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcat                                                                  >  1:837144/1‑150 (MQ=255)
                                                                                     atgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatc                                                                 >  2:798779/1‑150 (MQ=255)
                                                                                      tgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcg                                                                >  2:526312/1‑150 (MQ=255)
                                                                                        cgcATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcgca                                                              >  1:38764/1‑150 (MQ=255)
                                                                                          cATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcgcatc                                                            >  1:896761/1‑150 (MQ=255)
                                                                                           aTTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCTGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATcc                                                           <  2:1079998/150‑1 (MQ=255)
                                                                                               gCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCa                                                       >  1:1257945/1‑150 (MQ=255)
                                                                                               gCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCa                                                       <  2:259924/150‑1 (MQ=255)
                                                                                                  ggtggtGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAg                                                    <  2:655180/150‑1 (MQ=255)
                                                                                                   gtggtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAgg                                                   <  2:261810/150‑1 (MQ=255)
                                                                                                   gtggtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAgg                                                   <  1:1253342/150‑1 (MQ=255)
                                                                                                      gtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAg                                                >  1:1384216/1‑150 (MQ=255)
                                                                                                       tgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGa                                               >  1:1041490/1‑150 (MQ=255)
                                                                                                       tgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGa                                               <  2:165147/150‑1 (MQ=255)
                                                                                                        gAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAt                                              <  1:264022/150‑1 (MQ=255)
                                                                                                        gAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTCTAGGCTGGATAAGATCCGTCAGCATCGCATCCCGCaaacccacaa                                              >  2:473390/1‑143 (MQ=38)
                                                                                                         agtaagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAtc                                             <  1:453199/150‑1 (MQ=255)
                                                                                                          gtaagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAtct                                            >  1:145468/1‑150 (MQ=255)
                                                                                                            aagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCa                                          >  1:716989/1‑150 (MQ=255)
                                                                                                             agtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAg                                         >  2:267630/1‑150 (MQ=255)
                                                                                                             agtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAg                                         >  2:576276/1‑150 (MQ=255)
                                                                                                              gtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGc                                        <  1:810969/150‑1 (MQ=255)
                                                                                                              gtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGc                                        <  1:324512/150‑1 (MQ=255)
                                                                                                               taGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCg                                       <  2:834420/150‑1 (MQ=255)
                                                                                                               taGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCg                                       >  2:1396532/1‑150 (MQ=255)
                                                                                                                 gAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGAt                                     >  1:98844/1‑150 (MQ=255)
                                                                                                                  aaaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATa                                    >  1:104676/1‑150 (MQ=255)
                                                                                                                   aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                                   >  2:562751/1‑150 (MQ=255)
                                                                                                                   aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                                   >  1:715545/1‑150 (MQ=255)
                                                                                                                         cggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCgg                             >  2:1196772/1‑150 (MQ=255)
                                                                                                                            ccgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTAAGAGATAGCGACGTCtt                          <  2:223003/150‑1 (MQ=255)
                                                                                                                               ccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGt                       >  2:302799/1‑150 (MQ=255)
                                                                                                                                    ggTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAt                  >  2:801053/1‑150 (MQ=255)
                                                                                                                                    ggTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAt                  <  1:480884/150‑1 (MQ=255)
                                                                                                                                          ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATc            <  1:562751/150‑1 (MQ=255)
                                                                                                                                          ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATc            <  1:827911/150‑1 (MQ=255)
                                                                                                                                           tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCt           >  1:1029958/1‑150 (MQ=255)
                                                                                                                                            gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTg          >  1:510040/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgc         >  2:170355/1‑150 (MQ=39)
                                                                                                                                              cGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgcg        >  1:1009074/1‑150 (MQ=255)
                                                                                                                                                gATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgcgcg      <  2:1369649/150‑1 (MQ=255)
                                                                                                                                                  tGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGGTCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGcgcgcg    <  2:961213/150‑1 (MQ=255)
                                                                                                                                                    cggcgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtg  >  1:1173394/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTG  >  NC_000913/3984080‑3984377

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