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.1% 704.6 / ‑3.9 225L460V (CTG→GTG) bglBphospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base . (1/1);  major base C (90/131);  minor base A (2/0);  total (93/132)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 266 total aligned reads displayed.

ATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGAT  >  NC_000913/3899860‑3900038
                                                                                              |                                                                                      
 tGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGc                                                                                       >  1:702538/1‑93 (MQ=255)
 tGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGc                                                                                       <  2:2008426/93‑1 (MQ=255)
 tGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGc                                                                                       <  2:39794/93‑1 (MQ=255)
  ggTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGcc                                                                                      >  1:1009383/1‑93 (MQ=255)
    tGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGcccc                                                                                    <  1:2754232/93‑1 (MQ=255)
     gCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccc                                                                                   >  1:4876683/1‑93 (MQ=255)
     gCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccc                                                                                   >  1:5716798/1‑93 (MQ=255)
     gCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccc                                                                                   <  2:3744090/93‑1 (MQ=255)
      cAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGcccccc                                                                                  <  1:5276347/93‑2 (MQ=255)
      cAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGcccccc                                                                                  <  2:4534136/93‑2 (MQ=255)
      cAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGcccccc                                                                                  >  2:643363/1‑92 (MQ=255)
       aaagaaagAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccccg                                                                                 >  2:2039028/1‑91 (MQ=255)
       aaagaaagAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccccg                                                                                 >  1:2935355/1‑91 (MQ=255)
        aagaaagaATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccccgc                                                                                >  2:257152/1‑90 (MQ=255)
        aagaaagaATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGccccccgc                                                                                >  1:2905265/1‑90 (MQ=255)
           aaagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTc                                                                              >  1:5087380/1‑93 (MQ=255)
           aaagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTc                                                                              >  2:2895343/1‑93 (MQ=255)
           aaagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTc                                                                              <  1:2497651/93‑1 (MQ=255)
            aagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCt                                                                             >  1:315988/1‑93 (MQ=255)
            aagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCt                                                                             <  2:753242/93‑1 (MQ=255)
            aagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCt                                                                             <  1:5939046/93‑1 (MQ=255)
            aagaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCt                                                                             <  2:5542836/93‑1 (MQ=255)
             agaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCtt                                                                            <  1:6041041/93‑1 (MQ=255)
             agaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCtt                                                                            >  1:2007072/1‑93 (MQ=255)
              gaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTg                                                                           <  2:2280303/93‑1 (MQ=255)
              gaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTg                                                                           <  2:2666297/93‑1 (MQ=255)
                    aCGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTTATTAcc                                                                     >  2:3271226/1‑93 (MQ=255)
                    aCGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCAcc                                                                     >  1:3871444/1‑93 (MQ=255)
                       catcgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACctct                                                                  <  1:2667879/93‑1 (MQ=255)
                        atcgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGACCCCCGCGTCTTGATCACctcta                                                                 >  2:5109881/1‑92 (MQ=255)
                         tcgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACctctac                                                                >  1:910112/1‑91 (MQ=255)
                         tcgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACctctac                                                                >  2:1783434/1‑91 (MQ=255)
                          cgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACa                                                               >  1:2237133/1‑93 (MQ=255)
                           gcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACAt                                                              <  2:5988275/93‑1 (MQ=255)
                            catcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATa                                                             >  1:144052/1‑93 (MQ=255)
                             atcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATAc                                                            >  2:2743827/1‑93 (MQ=255)
                              tcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATAcc                                                           <  2:2366406/93‑1 (MQ=255)
                                 aGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATc                                                        <  2:2471539/93‑1 (MQ=255)
                                 aGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATc                                                        <  2:2914383/93‑1 (MQ=255)
                                 aGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATc                                                        <  2:322201/93‑1 (MQ=255)
                                  gACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATcc                                                       <  1:4271651/93‑1 (MQ=255)
                                    cTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGa                                                     <  1:4507710/93‑1 (MQ=255)
                                    cTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGa                                                     >  1:1853161/1‑93 (MQ=255)
                                     tGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGaa                                                    <  1:5470047/93‑1 (MQ=255)
                                      gTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCATCTCTACATACCATCCGAAg                                                   <  2:6083469/93‑1 (MQ=255)
                                        tCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCt                                                 <  1:1032921/93‑1 (MQ=255)
                                         cTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCCCCTCTACATGCCATCCGAAGCtt                                                <  1:3458301/93‑1 (MQ=255)
                                         cTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCtt                                                <  2:1262077/93‑1 (MQ=255)
                                         cTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCtt                                                >  1:1885149/1‑93 (MQ=255)
                                          tGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCttt                                               <  2:4527366/93‑1 (MQ=255)
                                          tGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCttt                                               <  1:5168283/93‑1 (MQ=255)
                                            aaTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTc                                             >  2:1491303/1‑93 (MQ=255)
                                             aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACGTCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                            <  2:500107/93‑1 (MQ=255)
                                             aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                            <  1:6259550/93‑1 (MQ=255)
                                             aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                            <  1:1420408/93‑1 (MQ=255)
                                              tGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCtt                                           <  2:3934260/93‑1 (MQ=255)
                                               gCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTg                                          <  2:4487088/93‑1 (MQ=255)
                                               gCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTg                                          <  2:1005047/93‑1 (MQ=255)
                                                cgacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGc                                         <  2:3672332/93‑1 (MQ=255)
                                                    gaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGGCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTCCgtgt                                     >  1:4096154/1‑93 (MQ=255)
                                                     aTAATTAAGGGGGTTTTATGGGTATTTTTTTTAATGACAGCCCCCCGCGTCCTGATCACCTCTTCATTACATTCGAAGCTTTTCTTGTGtttt                                    >  2:4834182/1‑93 (MQ=255)
                                                      tAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTc                                   <  2:5864785/93‑1 (MQ=255)
                                                       aaTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCt                                  >  2:5652613/1‑93 (MQ=255)
                                                       aaTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGACCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCt                                  >  2:582176/1‑93 (MQ=255)
                                                        aTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCtt                                 <  1:4748764/93‑1 (MQ=255)
                                                        aTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCtt                                 <  2:1333444/93‑1 (MQ=255)
                                                         ttAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtg                                <  2:4182892/93‑1 (MQ=255)
                                                          tAAGGTTGTTTAATGGTTTTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtgt                               <  2:5602877/93‑1 (MQ=255)
                                                            aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                             <  2:308588/93‑1 (MQ=255)
                                                            aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                             <  1:1166974/93‑1 (MQ=255)
                                                            aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                             <  1:5405247/93‑1 (MQ=255)
                                                            aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                             <  2:4505037/93‑1 (MQ=255)
                                                             agggCTTTTAAGGGTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAg                            <  2:1907876/90‑1 (MQ=255)
                                                                gCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGct                         <  2:4731745/93‑1 (MQ=255)
                                                                gCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGct                         >  2:436605/1‑93 (MQ=255)
                                                                gCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGct                         <  1:1326693/93‑1 (MQ=255)
                                                                 cTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGctt                        >  2:1975461/1‑93 (MQ=255)
                                                                  tttAATGTTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCCCCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttc                       <  2:1649686/93‑1 (MQ=255)
                                                                  tttAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttc                       >  1:6176899/1‑93 (MQ=255)
                                                                  tttAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttc                       >  1:3294696/1‑93 (MQ=255)
                                                                   ttAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcc                      <  2:1824844/93‑1 (MQ=255)
                                                                     aaTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                    >  1:3040339/1‑93 (MQ=255)
                                                                     aaTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                    >  1:5625482/1‑93 (MQ=255)
                                                                     aaTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                    <  2:2201849/93‑1 (MQ=255)
                                                                     aaTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                    <  1:2567652/93‑1 (MQ=255)
                                                                     aaTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                    <  1:4645994/93‑1 (MQ=255)
                                                                      aTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttccttc                   >  1:113484/1‑93 (MQ=255)
                                                                       tGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttccttcc                  >  1:5183004/1‑92 (MQ=255)
                                                                           tATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAt              <  2:4896629/93‑1 (MQ=255)
                                                                            aTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAtt             <  1:4181213/93‑1 (MQ=255)
                                                                             ttttttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCCTTa            >  2:5194043/1‑93 (MQ=255)
                                                                             ttttttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTa            <  1:2156055/93‑1 (MQ=255)
                                                                                tttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatca         >  2:5357175/1‑93 (MQ=255)
                                                                                 ttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatcat        <  1:5209320/93‑1 (MQ=255)
                                                                                  tttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatcatc       >  1:1554960/1‑93 (MQ=255)
                                                                                    tAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAc     >  1:3825990/1‑93 (MQ=255)
                                                                                    tAATGACAGCCCCCCGCGGCTTGAGCACCTCTACATACCATCCGAAGGTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCCATATCAGCAc     >  2:3526380/1‑93 (MQ=255)
                                                                                     aaTGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACg    >  1:3957303/1‑93 (MQ=255)
                                                                                     aaTGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACg    >  2:3788296/1‑93 (MQ=255)
                                                                                      aTGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGa   >  1:5808922/1‑93 (MQ=255)
                                                                                              |                                                                                      
ATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGAT  >  NC_000913/3899860‑3900038

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