Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,903,631 T→C Q24R (CAG→CGG)  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,6310TC98.3% 427.8 / ‑1.4 115Q24R (CAG→CGG) bglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (57/56);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 115 total aligned reads displayed.

CAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGT  >  NC_000913/3903485‑3903777
                                                                                                                                                   |                                                                                                                                                  
cAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGcc                                                                                                                                                 >  2:678474/1‑150 (MQ=255)
cAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCa                                                                                                                                                <  1:810760/151‑1 (MQ=255)
 aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGcc                                                                                                                                                 >  1:1175834/1‑149 (MQ=255)
 aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCa                                                                                                                                                <  1:693789/150‑1 (MQ=255)
 aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAg                                                                                                                                               <  2:1266691/151‑1 (MQ=255)
  aaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAgg                                                                                                                                              >  1:435997/1‑151 (MQ=255)
   aTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGc                                                                                                                                             >  1:1259181/1‑151 (MQ=255)
   aTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGc                                                                                                                                             <  1:1609957/151‑1 (MQ=255)
        gcgATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCtt                                                                                                                                        >  2:1409187/1‑151 (MQ=255)
          gATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTcc                                                                                                                                      >  2:1076839/1‑151 (MQ=255)
           atatCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTcc                                                                                                                                      >  1:1024186/1‑150 (MQ=255)
            tatCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCAc                                                                                                                                    >  2:668769/1‑151 (MQ=255)
               cTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTg                                                                                                                                 >  1:1473873/1‑151 (MQ=255)
                 tCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTg                                                                                                                                 <  2:340563/149‑1 (MQ=255)
                  ccGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGAtt                                                                                                                              <  2:69897/151‑1 (MQ=255)
                  ccGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGAt                                                                                                                               <  1:1282704/150‑1 (MQ=255)
                    gggTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTgg                                                                                                                            <  2:625700/151‑1 (MQ=255)
                      gTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGc                                                                                                                           >  2:581887/1‑150 (MQ=255)
                       tAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGc                                                                                                                           >  1:1685368/1‑149 (MQ=255)
                       tAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTGCCATCCTTCCCGCAGGCCCCTTCCACCTTAtttg                                                                                                                            >  2:938047/1‑145 (MQ=255)
                          cGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTg                                                                                                                      <  1:777937/151‑1 (MQ=255)
                          cGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTg                                                                                                                      <  1:1163290/151‑1 (MQ=255)
                            gTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGtt                                                                                                                    <  2:131830/151‑1 (MQ=255)
                            gTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGtt                                                                                                                    <  1:1158930/151‑1 (MQ=255)
                             tGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTg                                                                                                                   >  1:656971/1‑151 (MQ=255)
                              gATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTg                                                                                                                   <  2:905989/150‑1 (MQ=255)
                               aTAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTgcg                                                                                                                 >  2:1665468/1‑151 (MQ=255)
                                  aaaaTCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGcc                                                                                                               >  2:696447/1‑150 (MQ=255)
                                   aaaTCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCAc                                                                                                             >  2:1552306/1‑151 (MQ=255)
                                   aaaTCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCAc                                                                                                             >  1:703369/1‑151 (MQ=255)
                                       cGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCa                                                                                                         <  2:1670369/151‑1 (MQ=255)
                                           ggCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCaa                                                                                                        >  2:47178/1‑148 (MQ=255)
                                            gCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaa                                                                                                    <  1:563547/151‑1 (MQ=255)
                                             cGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaa                                                                                                    >  1:877567/1‑150 (MQ=255)
                                               aCATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATg                                                                                                 <  1:1409095/151‑1 (MQ=255)
                                                cATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGt                                                                                                >  2:1255421/1‑151 (MQ=255)
                                                cATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGt                                                                                                <  1:1608779/151‑1 (MQ=255)
                                                 aTCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGtt                                                                                               >  1:26394/1‑151 (MQ=255)
                                                   cTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTc                                                                                             >  1:180579/1‑151 (MQ=255)
                                                     ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGttt                                                                                              <  2:1460399/148‑1 (MQ=255)
                                                     ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGttt                                                                                              <  2:913243/148‑1 (MQ=255)
                                                       gATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCGCCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTgg                                                                                          <  1:366088/150‑1 (MQ=255)
                                                       gATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGa                                                                                         <  2:1425258/151‑1 (MQ=255)
                                                        atatTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTTCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCGCCCCAAAGAAATGTTTCTGGaa                                                                                        <  1:476778/151‑1 (MQ=255)
                                                        atatTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGaa                                                                                        <  2:656692/151‑1 (MQ=255)
                                                         tatTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGaa                                                                                        >  2:424414/1‑150 (MQ=255)
                                                            tctcTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGc                                                                                    >  1:183647/1‑151 (MQ=255)
                                                               ctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCt                                                                                   <  1:1185997/149‑1 (MQ=255)
                                                                  tCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGGGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTc                                                                                >  2:1264682/1‑149 (MQ=255)
                                                                   ccccAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCAt                                                                              >  1:993584/1‑150 (MQ=255)
                                                                   ccccAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATa                                                                             <  2:352781/151‑1 (MQ=255)
                                                                    cccAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCCGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACACCAAAGAAATGTTTCTGGAAATGCTTTCATa                                                                             >  2:1384299/1‑150 (MQ=255)
                                                                    cccAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATa                                                                             >  1:1081271/1‑150 (MQ=255)
                                                                     ccAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAAt                                                                           <  1:445528/151‑1 (MQ=255)
                                                                          aTGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACt                                                                      <  2:1319493/151‑1 (MQ=255)
                                                                          aTGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACt                                                                      >  1:1496408/1‑151 (MQ=255)
                                                                               cgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCttt                                                                 >  2:584756/1‑151 (MQ=255)
                                                                               cgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCttt                                                                 <  1:678416/151‑1 (MQ=255)
                                                                                    ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCg                                                            >  2:203665/1‑151 (MQ=255)
                                                                                    ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCg                                                            <  2:1018386/151‑1 (MQ=255)
                                                                                    ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCg                                                            >  1:1112262/1‑151 (MQ=255)
                                                                                     cATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGt                                                           >  2:832697/1‑151 (MQ=255)
                                                                                     cATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGt                                                           >  1:1502126/1‑151 (MQ=255)
                                                                                      aTTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGtt                                                          >  1:514269/1‑151 (MQ=255)
                                                                                       tttttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGtt                                                          >  2:1396963/1‑150 (MQ=255)
                                                                                          ttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGc                                                       <  2:1318105/150‑1 (MQ=255)
                                                                                            cccATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGaa                                                    <  1:64859/151‑1 (MQ=255)
                                                                                            cccATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGa                                                     >  2:640276/1‑150 (MQ=255)
                                                                                              cATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatg                                                  >  2:861929/1‑151 (MQ=255)
                                                                                               aTTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatga                                                 <  2:389583/151‑1 (MQ=255)
                                                                                                 tACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatgatg                                               >  2:1274095/1‑151 (MQ=255)
                                                                                                      cATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATaa                                          <  2:1218571/151‑1 (MQ=255)
                                                                                                      cATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATa                                           >  1:1616735/1‑150 (MQ=255)
                                                                                                      cATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATa                                           >  1:1345474/1‑150 (MQ=255)
                                                                                                        tGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAAc                                         >  2:1677874/1‑150 (MQ=255)
                                                                                                        tGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAAc                                         >  2:244967/1‑150 (MQ=255)
                                                                                                        tGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAAc                                         >  1:913277/1‑150 (MQ=255)
                                                                                                          aGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGc                                      <  2:1166706/151‑1 (MQ=255)
                                                                                                           ggCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCg                                     <  1:1076753/151‑1 (MQ=255)
                                                                                                             cTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGt                                   >  1:463299/1‑151 (MQ=255)
                                                                                                              tGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGtt                                  >  1:798868/1‑151 (MQ=255)
                                                                                                               gTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGt                                   <  1:1589220/149‑1 (MQ=255)
                                                                                                               gTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTc                                 <  1:1552210/151‑1 (MQ=255)
                                                                                                                tAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTc                                 >  2:508735/1‑150 (MQ=255)
                                                                                                                 aaaTCTGAGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCa                                <  1:1380920/150‑1 (MQ=255)
                                                                                                                        aGGTCGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCt                        <  1:584699/151‑1 (MQ=255)
                                                                                                                            cGAGATCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTAt                     <  2:1556116/150‑1 (MQ=255)
                                                                                                                                aTCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTgcg                <  1:1563464/151‑1 (MQ=255)
                                                                                                                                 tCCCTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTgcgc               >  1:676672/1‑151 (MQ=255)
                                                                                                                                   ccTTTGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccg             <  2:1260577/151‑1 (MQ=255)
                                                                                                                                     tttGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtg           >  1:1148589/1‑151 (MQ=255)
                                                                                                                                     tttGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtg           >  2:1152544/1‑151 (MQ=255)
                                                                                                                                      ttGCCATCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtg           >  2:479405/1‑150 (MQ=255)
                                                                                                                                           aTCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGccgtgccgtg      <  2:801341/150‑1 (MQ=255)
                                                                                                                                           aTCTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTgg     >  2:1253237/1‑151 (MQ=255)
                                                                                                                                              cTTCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTgggg   >  2:460575/1‑151 (MQ=255)
                                                                                                                                               ttCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGt  >  2:1599444/1‑151 (MQ=255)
                                                                                                                                                tCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCCCCTGCGCTTATCTGCGCCGTGCCGTggg    <  2:1008778/148‑1 (MQ=255)
                                                                                                                                                tCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTgggg   <  1:1235376/149‑1 (MQ=255)
                                                                                                                                                tCCCGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGt  <  2:877635/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGT  >  NC_000913/3903485‑3903777

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