Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,903,155 C→T G183S (GGC→AGC)  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,1550CT98.5% 255.3 / NA 65G183S (GGC→AGC) bglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/35);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACCAGCCCACCGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAG  >  NC_000913/3903008‑3903294
                                                                                                                                                   |                                                                                                                                           
aaaCCAGCCCACCGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGt                                                                                                                                          <  1:765803/151‑1 (MQ=255)
   ccAGCCCACCGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGa                                                                                                                                        >  2:876279/1‑150 (MQ=255)
    cAGCCCACCGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAAt                                                                                                                                      <  1:550644/151‑1 (MQ=255)
        ccACCGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGgcg                                                                                                                                  >  2:241882/1‑151 (MQ=255)
         cACCGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGgcgc                                                                                                                                 <  1:792654/151‑1 (MQ=255)
          aCCGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGcgcg                                                                                                                                >  1:248201/1‑151 (MQ=255)
           ccGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCgt                                                                                                                               <  2:275743/151‑1 (MQ=255)
           ccGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCgt                                                                                                                               >  1:508243/1‑151 (MQ=255)
              agaAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCgtgt                                                                                                                             >  1:316584/1‑150 (MQ=255)
               gaAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCgtgt                                                                                                                             <  2:193164/149‑1 (MQ=255)
                 aGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACg                                                                                                                         <  2:508189/151‑1 (MQ=255)
                      atTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACAt                                                                                                                     <  2:104517/150‑1 (MQ=255)
                         gCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTg                                                                                                                 <  1:224573/151‑1 (MQ=255)
                         gCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTg                                                                                                                 <  1:373805/151‑1 (MQ=255)
                           cGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGAt                                                                                                               <  1:398112/151‑1 (MQ=255)
                            gATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGAtt                                                                                                              >  2:19661/1‑151 (MQ=255)
                             aTTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGAttt                                                                                                             >  1:478625/1‑151 (MQ=255)
                              ttttCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTc                                                                                                            >  1:536796/1‑151 (MQ=255)
                                  cGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCAtt                                                                                                         <  1:588633/150‑1 (MQ=255)
                                  cGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCAtt                                                                                                         >  1:218048/1‑150 (MQ=255)
                                       cGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGt                                                                                                   <  2:781058/151‑1 (MQ=255)
                                         gggAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCa                                                                                                 <  1:483323/151‑1 (MQ=255)
                                         gggAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGACAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCa                                                                                                 <  1:153755/151‑1 (MQ=255)
                                          ggAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTc                                                                                                  >  2:42448/1‑149 (MQ=255)
                                           gAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGc                                                                                               >  2:808928/1‑151 (MQ=255)
                                             gCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCa                                                                                             <  2:514654/151‑1 (MQ=255)
                                              cAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCAt                                                                                            >  1:617121/1‑151 (MQ=255)
                                                ggCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATaa                                                                                          <  1:148730/151‑1 (MQ=255)
                                                ggCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATaa                                                                                          <  1:992062/151‑1 (MQ=255)
                                                 gCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAAc                                                                                         <  1:944678/151‑1 (MQ=255)
                                                    aTGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGcc                                                                                      >  2:771433/1‑151 (MQ=255)
                                                        cAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACtt                                                                                  <  2:519879/151‑1 (MQ=255)
                                                        cAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACtt                                                                                  >  1:668199/1‑151 (MQ=255)
                                                        cAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACtt                                                                                  >  1:513666/1‑151 (MQ=255)
                                                         aGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACttt                                                                                 <  1:107740/151‑1 (MQ=255)
                                                          ggCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACttt                                                                                 >  2:453747/1‑150 (MQ=255)
                                                           gCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTAt                                                                               >  1:598207/1‑151 (MQ=255)
                                                              ttAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGtt                                                                            >  2:171872/1‑151 (MQ=255)
                                                              ttAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGtt                                                                            >  2:698077/1‑151 (MQ=255)
                                                              ttAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGt                                                                             <  1:656056/150‑1 (MQ=255)
                                                                aaCTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTgg                                                                          <  1:241799/151‑1 (MQ=255)
                                                                aaCTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTg                                                                           <  2:638649/150‑1 (MQ=255)
                                                                 aCTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGt                                                                         >  1:791165/1‑151 (MQ=255)
                                                                    gCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGc                                                                      <  2:409156/151‑1 (MQ=255)
                                                                     cccGCGCCCTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAg                                                                       <  1:32668/149‑1 (MQ=255)
                                                                      ccGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCg                                                                     >  1:97254/1‑150 (MQ=255)
                                                                        gcgcACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGt                                                                  >  2:832978/1‑151 (MQ=255)
                                                                        gcgcACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGt                                                                  <  1:579641/151‑1 (MQ=255)
                                                                           cACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTaa                                                                >  1:695941/1‑150 (MQ=255)
                                                                               ggCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGAct                                                           <  2:319848/151‑2 (MQ=255)
                                                                                 cAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGt                                                         >  1:394814/1‑151 (MQ=255)
                                                                                         tGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCg                                                  <  1:348002/150‑1 (MQ=255)
                                                                                                  aTAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGa                                        <  2:44344/151‑1 (MQ=255)
                                                                                                    aGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAg                                      <  2:99628/151‑1 (MQ=255)
                                                                                                     gCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGt                                     <  2:481080/151‑1 (MQ=255)
                                                                                                         gATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATc                                 >  1:453496/1‑151 (MQ=255)
                                                                                                            aaaCTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGAt                              <  2:208296/151‑1 (MQ=255)
                                                                                                            aaaCTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGAt                              >  1:56015/1‑151 (MQ=255)
                                                                                                             aaCTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATg                             <  1:170110/151‑1 (MQ=255)
                                                                                                              aCTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATg                             >  2:46250/1‑150 (MQ=255)
                                                                                                                ttCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGAcc                          >  1:82421/1‑151 (MQ=255)
                                                                                                                                 cttctGCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCa            <  2:369881/148‑1 (MQ=255)
                                                                                                                                      gCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAAccgcc      <  2:676701/149‑1 (MQ=255)
                                                                                                                                      gCCAGCCCCACGCTCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGt    >  2:8614/1‑151 (MQ=255)
                                                                                                                                        cacccccGCGCCCGTGAATGGCGCTTGAAACAACTTGTTGTTTTCATTAAACGTCAGCCATAACCCCCTTTTTTTTTGGTGGCGATTAAAGCCCGTGCGGGCGTAATGCTCGAAGTTATCGATGCCCGCTCAATTACCCCAACCGCCGTAg  <  2:248174/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
AAACCAGCCCACCGAGAAGCATATTGCCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAG  >  NC_000913/3903008‑3903294

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: ‑