Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,903,621 G→A 9.2% G27G (GGC→GGT bglB ← cryptic phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,903,6210GA 9.2% 515.7 / 30.0 163G27G (GGC→GGTbglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base G (78/70);  new base A (8/7);  total (86/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

AAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGT  >  NC_000913/3903472‑3903767
                                                                                                                                                     |                                                                                                                                                  
aaGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTg                                                                                                                                                    <  2:2890952/150‑1 (MQ=255)
aaGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTg                                                                                                                                                    <  2:2890967/150‑1 (MQ=255)
aaGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTg                                                                                                                                                    <  1:1448895/150‑1 (MQ=255)
  gCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGcc                                                                                                                                                  <  1:2350265/150‑1 (MQ=255)
  gCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGcc                                                                                                                                                  >  2:1888130/1‑150 (MQ=255)
  gCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGcc                                                                                                                                                  <  2:3037562/150‑1 (MQ=255)
    ccATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCAt                                                                                                                                                >  1:2035164/1‑150 (MQ=255)
    ccATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCAt                                                                                                                                                >  2:2986613/1‑150 (MQ=255)
        gtcGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTc                                                                                                                                            <  2:2152478/149‑1 (MQ=255)
        ctcGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTc                                                                                                                                            >  1:274230/1‑150 (MQ=255)
          cGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCt                                                                                                                                          <  2:2205747/150‑1 (MQ=255)
           ggCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTg                                                                                                                                         <  2:3065283/150‑1 (MQ=255)
            gCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGc                                                                                                                                        <  2:2814555/150‑1 (MQ=255)
                aTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAgg                                                                                                                                    >  1:3033696/1‑150 (MQ=255)
                  aaCGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCa                                                                                                                                  <  1:1332611/150‑1 (MQ=255)
                      cgATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCtt                                                                                                                              >  1:766659/1‑150 (MQ=255)
                      cgATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCtt                                                                                                                              >  1:766648/1‑150 (MQ=255)
                       gATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTc                                                                                                                             <  2:725488/150‑1 (MQ=255)
                       gATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTc                                                                                                                             >  2:1524888/1‑150 (MQ=255)
                       gATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCGTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTc                                                                                                                             <  2:2599164/150‑1 (MQ=255)
                         tatCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCa                                                                                                                           <  1:1788551/150‑1 (MQ=255)
                          atCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCAc                                                                                                                          <  1:188508/150‑1 (MQ=255)
                            cTTCCGTGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATTCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCt                                                                                                                        <  2:2000650/150‑1 (MQ=255)
                            cTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCt                                                                                                                        <  2:1187572/150‑1 (MQ=255)
                              tCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGa                                                                                                                      <  1:1846014/150‑1 (MQ=255)
                               ccGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGAt                                                                                                                     <  1:1674590/150‑1 (MQ=255)
                                cGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGAtt                                                                                                                    <  2:435665/150‑1 (MQ=255)
                                   gTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGc                                                                                                                 >  2:2419177/1‑150 (MQ=255)
                                   gTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGc                                                                                                                 >  2:2317641/1‑150 (MQ=255)
                                   gTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGc                                                                                                                 >  2:1042606/1‑150 (MQ=255)
                                   gTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGc                                                                                                                 >  2:1042623/1‑150 (MQ=255)
                                    tAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCa                                                                                                                <  1:679741/150‑1 (MQ=255)
                                      aCGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGc                                                                                                              <  2:2731109/150‑1 (MQ=255)
                                        ggTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGcc                                                                                                                                       <  1:464321/123‑1 (MQ=255)
                                        ggTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTg                                                                                                            >  1:803534/1‑150 (MQ=255)
                                          tGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGtt                                                                                                          <  2:2266573/150‑1 (MQ=255)
                                               aaaaTCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGcc                                                                                                     <  1:2340831/150‑1 (MQ=255)
                                                   tCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCAccc                                                                                                 <  1:973551/150‑1 (MQ=255)
                                                      aTGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTGTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCaa                                                                                              >  2:1170993/1‑150 (MQ=255)
                                                         gCGACATCTTTGATATTCTCTTGCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGa                                                                                           <  1:772325/150‑1 (MQ=255)
                                                           gACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaaa                                                                                         <  1:2178764/150‑1 (MQ=255)
                                                           gACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaaa                                                                                         >  2:1998374/1‑150 (MQ=255)
                                                           gACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTCCGCCACCCCAAAGaaa                                                                                         >  2:1653358/1‑150 (MQ=255)
                                                            aCATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGCCAGCTGTTGCGCCACCCCAAAGAAAt                                                                                        >  2:88097/1‑150 (MQ=255)
                                                               tCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGtt                                                                                     >  1:42131/1‑150 (MQ=255)
                                                                 tttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTc                                                                                   >  2:1743688/1‑150 (MQ=255)
                                                                 tttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTc                                                                                   >  2:1602390/1‑150 (MQ=255)
                                                                  ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATGTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                                                                  >  2:879400/1‑150 (MQ=255)
                                                                  ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                                                                  <  2:594836/150‑1 (MQ=255)
                                                                  ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                                                                  <  2:529679/150‑1 (MQ=255)
                                                                  ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                                                                  <  2:462279/150‑1 (MQ=255)
                                                                   tGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTg                                                                                 >  1:132338/1‑150 (MQ=255)
                                                                     atatTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGa                                                                               >  2:1049572/1‑150 (MQ=255)
                                                                        ttCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAt                                                                            <  1:1995510/150‑1 (MQ=255)
                                                                           tctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCt                                                                         >  1:1017589/1‑150 (MQ=255)
                                                                           tctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCt                                                                         <  1:1599667/150‑1 (MQ=255)
                                                                               tCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCa                                                                     <  2:2567840/150‑1 (MQ=255)
                                                                               tCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCa                                                                     <  2:2559693/150‑1 (MQ=255)
                                                                                   cAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAAt                                                                 >  2:79168/1‑150 (MQ=255)
                                                                                    aGGATGCGTGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAAtt                                                                <  1:2271515/150‑1 (MQ=255)
                                                                                     ggATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTa                                                               <  1:2693608/150‑1 (MQ=255)
                                                                                      gATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTaa                                                              <  1:2852691/150‑1 (MQ=255)
                                                                                       aTGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                                                             >  1:1028475/1‑150 (MQ=255)
                                                                                       aTGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                                                             >  1:657120/1‑150 (MQ=255)
                                                                                       aTGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                                                             <  1:363138/150‑1 (MQ=255)
                                                                                           gcgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCt                                                         >  2:2380391/1‑150 (MQ=255)
                                                                                            cgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtt                                                        <  1:154462/150‑1 (MQ=255)
                                                                                            cgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtt                                                        <  1:746775/150‑1 (MQ=255)
                                                                                             ggTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCttt                                                       >  1:1501822/1‑150 (MQ=255)
                                                                                             ggTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCttt                                                       <  1:835676/150‑1 (MQ=255)
                                                                                             ggTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCttt                                                       >  1:1501823/1‑150 (MQ=255)
                                                                                             ggTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCAGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAATCCttt                                                       >  1:2036437/1‑150 (MQ=255)
                                                                                               ttCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTa                                                     >  1:929063/1‑150 (MQ=255)
                                                                                               ttCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTa                                                     >  2:3054643/1‑150 (MQ=255)
                                                                                                    tttttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGtt                                                >  2:2990229/1‑150 (MQ=255)
                                                                                                     ttttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTa                                               <  2:1252566/150‑1 (MQ=255)
                                                                                                     ttttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTa                                               <  1:2356669/150‑1 (MQ=255)
                                                                                                      tttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAg                                              >  1:580725/1‑150 (MQ=255)
                                                                                                      tttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAg                                              >  1:271014/1‑150 (MQ=255)
                                                                                                       ttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGc                                             >  2:1417223/1‑150 (MQ=255)
                                                                                                       ttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGc                                             >  1:2028891/1‑150 (MQ=255)
                                                                                                              tACGCCATGAGGCTGTAACTCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatgat                                      <  1:1868616/150‑1 (MQ=255)
                                                                                                                cGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATgg                                    >  2:184697/1‑150 (MQ=255)
                                                                                                                cGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATgg                                    >  1:2113775/1‑150 (MQ=255)
                                                                                                                cGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATgg                                    >  2:184693/1‑150 (MQ=255)
                                                                                                                      gAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACa                              <  2:32988/150‑1 (MQ=255)
                                                                                                                      gAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACa                              <  2:1185367/150‑1 (MQ=255)
                                                                                                                       aGGCTGTAAATCTGAGGTCGAGATCCCTTTGTCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAg                             <  1:3057556/150‑1 (MQ=255)
                                                                                                                       aGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAg                             <  1:2806928/150‑1 (MQ=255)
                                                                                                                        ggCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGc                            >  2:2318701/1‑150 (MQ=255)
                                                                                                                          cTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCgg                          >  2:2753274/1‑150 (MQ=255)
                                                                                                                          cTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCgg                          <  2:467988/150‑1 (MQ=255)
                                                                                                                                   tGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgc                 >  2:2956983/1‑150 (MQ=255)
                                                                                                                                     aGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcgc               >  2:780219/1‑150 (MQ=255)
                                                                                                                                       gTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCtt             >  1:2188136/1‑150 (MQ=255)
                                                                                                                                          gagaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATc          >  1:2706499/1‑150 (MQ=255)
                                                                                                                                             aTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTgc       <  1:110910/150‑1 (MQ=255)
                                                                                                                                               cccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTgcgc     <  1:184693/150‑1 (MQ=255)
                                                                                                                                                ccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgcc    >  2:205819/1‑150 (MQ=255)
                                                                                                                                                 cTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccg   >  1:2624199/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGT  >  NC_000913/3903472‑3903767

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