Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,899,975 G→A A452V (GCA→GTA)  bglB ← phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,899,9750GA99.7% 2337.6 / NA 642A452V (GCA→GTA) bglBphospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (294/345);  minor base C (2/0);  total (297/345)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 675 total aligned reads displayed.

TCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCTTA  >  NC_000913/3899885‑3900064
                                                                                             |                                                                                         
 cgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACa                                                                                         <  1:7417526/93‑1 (MQ=255)
 cgcatcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACa                                                                                         <  2:9289844/93‑1 (MQ=255)
    atcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATAc                                                                                      <  2:1622061/93‑1 (MQ=255)
     tcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATAcc                                                                                     >  2:4706980/1‑93 (MQ=255)
        aGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATc                                                                                  <  2:3009212/93‑1 (MQ=255)
         gACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATcc                                                                                 <  1:7587738/93‑1 (MQ=255)
         gACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATcc                                                                                 <  1:6927916/93‑1 (MQ=255)
         gACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATcc                                                                                 >  2:6381173/1‑93 (MQ=255)
           cTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGaa                                                                              <  2:3372111/93‑1 (MQ=255)
            tGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAg                                                                             >  1:2943420/1‑93 (MQ=255)
            tGTTCGGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGaa                                                                              <  2:7925688/93‑1 (MQ=255)
             gTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAg                                                                             >  2:2718085/1‑93 (MQ=255)
               tCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCt                                                                           >  1:6540032/1‑93 (MQ=255)
               tCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCt                                                                           >  2:9161614/1‑93 (MQ=255)
                 tGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCttt                                                                         >  1:676109/1‑93 (MQ=255)
                 tGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCttt                                                                         >  1:2060270/1‑93 (MQ=255)
                  gAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCtttt                                                                        >  1:880594/1‑93 (MQ=255)
                    aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                                                      <  2:1072834/93‑1 (MQ=255)
                    aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTTCATGCCATCCGAAGCTTTTCtt                                                                     >  2:113053/1‑93 (MQ=255)
                     tGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCtt                                                                     <  2:8377195/93‑1 (MQ=255)
                     tGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTCCATACCATACGAAGCTTTTCTTg                                                                    >  2:7178962/1‑93 (MQ=255)
                     tGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTg                                                                    >  2:4205961/1‑93 (MQ=255)
                       cgacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGc                                                                   <  1:3498409/93‑1 (MQ=255)
                       cgacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGc                                                                   <  1:3131154/93‑1 (MQ=255)
                       cgacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCg                                                                  <  1:509676/93‑1 (MQ=255)
                         acgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGACCCCCGCGTCCTGATCACCTCTTCATACCATGCCAAGCTTTTCTTGcct                                                                 >  2:4048441/1‑91 (MQ=255)
                         acgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCgtg                                                                <  1:4381498/93‑1 (MQ=255)
                           gaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGtt                                                              >  1:8670558/1‑93 (MQ=255)
                            aTAATTAAGGTGCTTTAATGGTTATTTTTTTTTAATGCCACCCCCCCCGTTTTGACACCTTCTACATACAATCCGCGGCTTTTCTGGCGGGtt                                                              >  2:2457044/1‑93 (MQ=255)
                             tAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTc                                                             >  1:9876851/1‑93 (MQ=255)
                              aaTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCt                                                            >  1:9246662/1‑93 (MQ=255)
                              aaTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCt                                                            <  2:9598467/93‑1 (MQ=255)
                               aTTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtg                                                          <  1:8278973/93‑1 (MQ=255)
                                ttAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtg                                                          >  2:9668060/1‑93 (MQ=255)
                                  aaGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                                                       >  1:2443151/1‑93 (MQ=255)
                                   aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                                                       <  2:8217022/93‑1 (MQ=255)
                                   aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                                                       <  1:4132252/93‑1 (MQ=255)
                                   aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                                                       <  2:36290/93‑1 (MQ=255)
                                   aGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                                                       <  2:6194140/93‑1 (MQ=255)
                                      tGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGc                                                    <  2:5650767/93‑1 (MQ=255)
                                      tGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGct                                                   <  1:3668866/93‑1 (MQ=255)
                                      tGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGct                                                   >  2:2341210/1‑93 (MQ=255)
                                       gCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGctt                                                  >  1:1250852/1‑93 (MQ=255)
                                       gCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGctt                                                  >  2:4816116/1‑93 (MQ=255)
                                        cTTTAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGctt                                                  <  2:4289220/93‑1 (MQ=255)
                                         tttAATGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttc                                                 <  2:1819717/93‑1 (MQ=255)
                                          ttAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcct                                               <  2:2208114/93‑1 (MQ=255)
                                          ttAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcct                                               <  1:778333/93‑1 (MQ=255)
                                           tAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTTTTGCGTGTTCTTGTGAGGcttcctt                                              >  1:1457609/1‑93 (MQ=255)
                                           tAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                                              >  2:2007881/1‑93 (MQ=255)
                                            aaTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcttt                                              >  2:8963767/1‑93 (MQ=255)
                                            aaTGGTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                                              <  1:8510937/93‑1 (MQ=255)
                                            aaTGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttccttc                                             <  1:9451852/93‑1 (MQ=255)
                                                gTTATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGcc                                          >  1:3941284/1‑93 (MQ=255)
                                                 ttATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCa                                         <  1:2678208/93‑1 (MQ=255)
                                                  tATTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAt                                        <  2:4487264/93‑1 (MQ=255)
                                                   aTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAtt                                       <  1:1101654/93‑1 (MQ=255)
                                                   aTTTTTTTTAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAtt                                       >  1:3286060/1‑93 (MQ=255)
                                                   aTTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTa                                      >  1:10046127/1‑93 (MQ=255)
                                                   aTTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTa                                      <  1:9045461/93‑1 (MQ=255)
                                                     ttttttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTa                                      <  1:7065225/93‑1 (MQ=255)
                                                     ttttttttAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTat                                     <  1:5797933/93‑1 (MQ=255)
                                                      tttttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTat                                     >  1:6321877/1‑93 (MQ=255)
                                                        tttttAATGACAGCCCCCCGCGTCTTGATCCCCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatca                                   <  1:4928457/93‑1 (MQ=255)
                                                        tttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatca                                   >  2:6084358/1‑93 (MQ=255)
                                                         ttttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatcat                                  >  1:2629705/1‑93 (MQ=255)
                                                         ttttAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatcatc                                 <  1:794302/93‑1 (MQ=255)
                                                          tttAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTAtcatca                                <  1:5971876/93‑1 (MQ=255)
                                                           ttAATGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGATTCCTTCGCCATTAtcatca                                <  1:5218160/93‑1 (MQ=255)
                                                           ttAATGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAc                               >  1:973263/1‑93 (MQ=255)
                                                              aTGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGAt                            <  2:314532/93‑1 (MQ=255)
                                                              aTGACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGAt                            >  2:1241622/1‑93 (MQ=255)
                                                               tGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGAt                            <  1:9136028/93‑1 (MQ=255)
                                                               tGACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGAt                            >  2:5020554/1‑93 (MQ=255)
                                                                gACAGCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATc                           <  2:7559755/93‑1 (MQ=255)
                                                                gACAGCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATcc                          >  2:4881875/1‑93 (MQ=255)
                                                                  cagccccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCca                         >  2:532538/5‑93 (MQ=255)
                                                                   atccccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCcac                        <  1:614299/90‑1 (MQ=255)
                                                                   agccccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCcac                        <  1:1708854/90‑1 (MQ=255)
                                                                    gccccccGCGGATTCATCACCCTTACATACCCTCCCAAACTTTTCTTGCGTGTTCTTGTGAGGCTTACCTTGTCCTTTTTATCACGTACcaca                       >  2:7294139/3‑93 (MQ=255)
                                                                      cccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACata                     >  1:1606384/1‑93 (MQ=255)
                                                                      cccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACata                     <  1:2229245/93‑1 (MQ=255)
                                                                       ccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataa                    >  2:3859898/1‑93 (MQ=255)
                                                                       ccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataa                    <  1:8617018/93‑1 (MQ=255)
                                                                        cccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaa                   >  1:4707515/1‑93 (MQ=255)
                                                                         ccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaat                  <  2:3281998/93‑1 (MQ=255)
                                                                         ccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaat                  <  1:9587541/93‑1 (MQ=255)
                                                                           gcgTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaataa                >  2:8809307/1‑93 (MQ=255)
                                                                             gTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAg              <  1:3398044/93‑1 (MQ=255)
                                                                              tCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGc             <  1:1790278/93‑1 (MQ=255)
                                                                              tCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGc             >  2:589224/1‑93 (MQ=255)
                                                                              gcTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGc             <  1:5603635/92‑1 (MQ=255)
                                                                               cTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGcc            <  2:8085326/93‑1 (MQ=255)
                                                                                 tgatgaCCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATTAATAAAGCCGt          <  2:8608459/93‑1 (MQ=255)
                                                                                 tGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGt          <  1:5451562/93‑1 (MQ=255)
                                                                                 tGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGt          <  1:8073761/93‑1 (MQ=255)
                                                                                  gATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTa         <  2:6787566/93‑1 (MQ=255)
                                                                                       ccTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCt    >  1:3847658/1‑93 (MQ=255)
                                                                                          tctACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCTTa  >  1:7236830/1‑93 (MQ=255)
                                                                                          tctACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCTTa  >  2:6913363/1‑93 (MQ=255)
                                                                                             |                                                                                         
TCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCTTA  >  NC_000913/3899885‑3900064

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