Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,899,957 (C)5→6 coding (1373/1413 nt) bglB ← phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,899,9521.C99.6% 765.7 / NA 247L460V (CTG→GTG) bglBphospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (114/131);  minor base A (1/0);  total (116/131)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

ATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATC  >  NC_000913/3899860‑3900039
                                                                                              |                                                                                        
aTGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAg                                                                                           >  1:3734793/1‑93 (MQ=255)
 tGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGc                                                                                         >  2:4468275/1‑93 (MQ=255)
   gTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTTAATGACAGcc                                                                                       <  1:1044472/93‑1 (MQ=255)
   gTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccc                                                                                      <  1:4093037/93‑1 (MQ=255)
    tGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGcccc                                                                                     >  2:3851203/1‑93 (MQ=255)
     gCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccc                                                                                    >  1:868127/1‑93 (MQ=255)
        aagaaagaATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccccgc                                                                                 >  1:3895456/1‑90 (MQ=255)
           aaagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTc                                                                               <  1:717988/93‑1 (MQ=255)
           aaagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTc                                                                               <  1:2745191/93‑1 (MQ=255)
            aagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCt                                                                              <  1:3355156/93‑1 (MQ=255)
             agaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCtt                                                                             >  2:1087079/1‑93 (MQ=255)
                aTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGAt                                                                          >  1:955582/1‑93 (MQ=255)
                aTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGAt                                                                          <  1:2168339/93‑1 (MQ=255)
                 ttAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATc                                                                         <  1:2345369/93‑1 (MQ=255)
                 ttAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATc                                                                         <  2:4577442/93‑1 (MQ=255)
                 ttAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGACCCCCGCGTCTTGATc                                                                         >  2:372632/1‑93 (MQ=255)
                    aCGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCAcc                                                                      >  2:299594/1‑93 (MQ=255)
                    aCGCATCGCATCCAGACTGTTCTGAATGCGACGATAATCAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATGAcc                                                                      <  2:4132984/93‑1 (MQ=255)
                         tcgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACctctac                                                                 >  2:1984612/1‑91 (MQ=255)
                         tcgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACctctac                                                                 <  2:4437492/93‑3 (MQ=255)
                         tcgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCGTGATCACctctac                                                                 <  2:4606044/93‑3 (MQ=255)
                          cgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCCCCTCTACa                                                                <  2:3220065/93‑1 (MQ=255)
                          cgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACa                                                                <  1:739464/93‑1 (MQ=255)
                           gcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACAt                                                               <  2:1587424/93‑1 (MQ=255)
                              tcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATAcc                                                            <  1:2116684/93‑1 (MQ=255)
                              tcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCACCCCGCGTCTTGATCACCTCTACATAcc                                                            >  2:2706434/1‑93 (MQ=255)
                                cAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCAt                                                          >  2:2571951/1‑93 (MQ=255)
                                cAGACTGTTCTGAATGCGACCATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCAt                                                          <  2:4216438/93‑1 (MQ=255)
                                 aGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATc                                                         >  2:1542191/1‑93 (MQ=255)
                                  gACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATcc                                                        >  1:1189871/1‑93 (MQ=255)
                                    cTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGa                                                      >  1:2151114/1‑93 (MQ=255)
                                     tGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGaa                                                     <  1:1665871/93‑1 (MQ=255)
                                     tGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGaa                                                     >  1:2297562/1‑93 (MQ=255)
                                      gTTCTGAATGCGCCGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAg                                                    <  2:2127244/93‑1 (MQ=255)
                                       ttCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGc                                                   <  1:4650627/93‑1 (MQ=255)
                                        tCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGc                                                   >  2:434347/1‑93 (MQ=255)
                                          tGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCttt                                                <  2:1403260/93‑1 (MQ=255)
                                          tGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCCTTGTCACCCCTCCATACCATCCCAAACttt                                                >  2:866895/1‑93 (MQ=255)
                                            aaTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTc                                              <  1:4155416/93‑1 (MQ=255)
                                            aaTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTc                                              >  1:2389711/1‑93 (MQ=255)
                                             aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                             <  2:1897395/93‑1 (MQ=255)
                                              tGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCtt                                            <  1:3970843/93‑1 (MQ=255)
                                                 gacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCg                                         >  1:3088452/1‑93 (MQ=255)
                                                 gacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCCTGCg                                         >  2:4147519/1‑93 (MQ=255)
                                                  acgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCgt                                        <  1:2572053/93‑1 (MQ=255)
                                                  acgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCgt                                        <  1:2720593/93‑1 (MQ=255)
                                                    gaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCgtgt                                      <  2:3689286/93‑1 (MQ=255)
                                                    gaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCgtgt                                      >  2:3650896/1‑93 (MQ=255)
                                                    gaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATAACATCCGAAGCTTTTCTTGCgtgt                                      >  1:2473415/1‑93 (MQ=255)
                                                     aTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGtt                                     >  2:3115525/1‑93 (MQ=255)
                                                     aTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGtt                                     >  2:3326156/1‑93 (MQ=255)
                                                      tAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTc                                    >  1:4417114/1‑93 (MQ=255)
                                                      tAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTc                                    <  2:1370884/93‑1 (MQ=255)
                                                      tAATTAAGGGGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGATTTTCTTGCGTGTTc                                    >  1:2070209/1‑93 (MQ=255)
                                                       aaTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCt                                   >  1:1313158/1‑93 (MQ=255)
                                                       aaTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCt                                   <  1:356517/93‑1 (MQ=255)
                                                       aaTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCt                                   <  2:2650330/93‑1 (MQ=255)
                                                       aaTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATAACATCCGAAGCTTTTCTTGCGTGTTCt                                   >  1:1619398/1‑93 (MQ=255)
                                                        aTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCtt                                  <  2:1032926/93‑1 (MQ=255)
                                                         ttAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtg                                 <  2:2820495/93‑1 (MQ=255)
                                                         ttAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtg                                 <  1:41703/93‑1 (MQ=255)
                                                         ttAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtg                                 <  1:55942/93‑1 (MQ=255)
                                                           aaGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtgtg                               <  2:1461615/93‑1 (MQ=255)
                                                            aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtgtgg                              >  2:1654547/1‑92 (MQ=255)
                                                            aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                              >  1:3523057/1‑93 (MQ=255)
                                                             ggTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAg                             >  2:2206073/1‑93 (MQ=255)
                                                               tGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGc                           >  2:1926979/1‑93 (MQ=255)
                                                                gCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGct                          >  2:496860/1‑93 (MQ=255)
                                                                gCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGCCTTGATTACCTCTACATACCATCCGAAGCTTTTTTTGCGTGTTCTTGTGAGGct                          >  2:3283563/1‑93 (MQ=255)
                                                                  tttAATTGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCCCCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttc                        <  1:1003761/93‑1 (MQ=255)
                                                                  tttAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttc                        >  1:788170/1‑93 (MQ=255)
                                                                   ttAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcc                       <  2:2053619/93‑1 (MQ=255)
                                                                     taTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                     <  1:461197/92‑1 (MQ=255)
                                                                       tGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGCTCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCg                   >  2:769766/1‑93 (MQ=255)
                                                                        ggTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGTGGCTTCCTTCGc                  >  1:4297172/1‑93 (MQ=255)
                                                                        ggTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGc                  >  1:4540985/1‑93 (MQ=255)
                                                                        ggTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGc                  >  2:3565208/1‑93 (MQ=255)
                                                                         gTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGcc                 <  1:1407534/93‑1 (MQ=255)
                                                                          ttATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCa                >  2:3390571/1‑93 (MQ=255)
                                                                          ttATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCa                >  1:1357019/1‑93 (MQ=255)
                                                                          ttATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCa                <  2:4020417/93‑1 (MQ=255)
                                                                          ttATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCa                <  2:3858187/93‑1 (MQ=255)
                                                                            aTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAtt              <  2:3099523/93‑1 (MQ=255)
                                                                            aTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAtt              >  2:1216481/1‑93 (MQ=255)
                                                                            aTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAtt              <  2:3494967/93‑1 (MQ=255)
                                                                            aTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAtt              >  1:1628705/1‑93 (MQ=255)
                                                                              ttttttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTa             <  2:1712964/93‑1 (MQ=255)
                                                                               tttttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTat            <  1:3875543/93‑1 (MQ=255)
                                                                               tttttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTat            <  1:1767476/93‑1 (MQ=255)
                                                                                ttttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatc           <  1:3039535/93‑1 (MQ=255)
                                                                                  ttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatcat         <  1:3156492/93‑1 (MQ=255)
                                                                                    ttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCTCCATTAtcatca       >  2:1281559/1‑93 (MQ=255)
                                                                                    ttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTAtcatca       <  2:2416575/93‑1 (MQ=255)
                                                                                    ttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTAtcatca       >  2:582912/1‑93 (MQ=255)
                                                                                    ttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTAtcatca       >  1:3654593/1‑93 (MQ=255)
                                                                                     tAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAc      >  1:485158/1‑93 (MQ=255)
                                                                                      aaTGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACg     <  1:2835904/93‑1 (MQ=255)
                                                                                      aaTGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACg     >  1:1956762/1‑93 (MQ=255)
                                                                                      aaTGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACg     >  1:1503407/1‑93 (MQ=255)
                                                                                       aTGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGa    <  1:2367261/93‑1 (MQ=255)
                                                                                              |                                                                                        
ATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATC  >  NC_000913/3899860‑3900039

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