Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,689,996 A→C L89* (TTA→TGA)  bcsZ ← endo‑1,4‑D‑glucanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,689,9960AC97.0% 546.8 / 3.2 167L89* (TTA→TGA) bcsZendo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (85/74);  total (89/78)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 172 total aligned reads displayed.

CATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTAAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACATGCCGTAACTTT  >  NC_000913/3689899‑3690092
                                                                                                 |                                                                                                
cATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTcaa                                                                                                <  1:1273520/100‑1 (MQ=255)
  tCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGa                                                                                              >  1:2877689/1‑100 (MQ=255)
  tCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGa                                                                                              >  1:1756084/1‑100 (MQ=255)
    cAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAAc                                                                                            >  2:2620956/1‑100 (MQ=255)
       aCATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCt                                                                                         >  1:2225297/1‑100 (MQ=255)
       aCATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTGAAAAGAAcgtt                                                                                         <  1:3651840/100‑4 (MQ=255)
         aTCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGa                                                                                       >  2:817317/1‑100 (MQ=255)
          tCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAg                                                                                      >  1:3577911/1‑100 (MQ=255)
           cACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGc                                                                                     <  2:270979/100‑1 (MQ=255)
            aCCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCg                                                                                    <  2:1268287/100‑1 (MQ=255)
            aCCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCg                                                                                    >  2:1751759/1‑100 (MQ=255)
             ccATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCga                                                                                   >  1:3835164/1‑100 (MQ=255)
                tCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGAt                                                                                <  1:874523/100‑1 (MQ=255)
                tCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGAt                                                                                <  1:701866/100‑1 (MQ=255)
                   gAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTAAAGAACCCTGAGCGAGATTGt                                                                             >  1:2794558/1‑100 (MQ=255)
                    aGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGtt                                                                            <  1:147936/100‑1 (MQ=255)
                     ggCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTc                                                                           <  1:3364632/100‑1 (MQ=255)
                     ggCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTc                                                                           <  2:48653/100‑1 (MQ=255)
                      gCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCt                                                                          >  2:980180/1‑100 (MQ=255)
                       cccaaTTGCTGTCCAGCACACCCACCTTACTGTACTCTTTCTTCCACCACAGCCAGGCGAACAAACGTTCTTTTGAAGAACCCTGAGCGAGGTTGTTCTg                                                                         <  2:556190/97‑1 (MQ=255)
                       ccGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTg                                                                         >  2:1433055/1‑100 (MQ=255)
                       ccGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTg                                                                         >  2:3027669/1‑100 (MQ=255)
                       ccGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTg                                                                         >  2:3463577/1‑100 (MQ=255)
                       ccGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTg                                                                         >  1:3996192/1‑100 (MQ=255)
                         gAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCg                                                                       >  2:2550936/1‑100 (MQ=255)
                             tGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCa                                                                   >  1:594426/1‑100 (MQ=255)
                              gCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTATTTTCCAGAACCCTGAGCGCGATTGTTCTGCGTCCAg                                                                  >  2:777447/1‑100 (MQ=255)
                               cTGTCCAGCACTTCCCGCTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAATAACCCTGAGCGAGATTGTTCTGCGTCCAGt                                                                 <  1:3354352/100‑1 (MQ=255)
                               cTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGt                                                                 >  1:3128784/1‑100 (MQ=255)
                                tGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGtn                                                                <  1:3237756/100‑2 (MQ=255)
                                 gTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCg                                                               >  2:3629430/1‑100 (MQ=255)
                                  tCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGCTCTGCGTCCAGTCga                                                              <  2:3034700/100‑1 (MQ=255)
                                   ccAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCgag                                                             >  1:3495550/1‑100 (MQ=255)
                                   ccAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCgag                                                             <  1:403303/100‑1 (MQ=255)
                                   ccAGCACTTCCCACTTACTGTTCTCTTACTTGCCCCACAGCCAGGAGGGCATACTTTCTTTGCAAGAACCCTGTGCGGGATTGTTCTGCGGCCCGTCgag                                                             >  2:2681568/1‑100 (MQ=255)
                                    cAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCgaga                                                            >  1:25194/1‑100 (MQ=255)
                                    cAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCgaga                                                            >  1:1453034/1‑100 (MQ=255)
                                       cACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAata                                                         <  1:2590533/100‑1 (MQ=255)
                                       cACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAata                                                         <  2:805072/100‑1 (MQ=255)
                                         cTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAAtatt                                                       >  1:112188/1‑100 (MQ=255)
                                          ttCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAAtatta                                                      >  2:1586015/1‑100 (MQ=255)
                                           tCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAAtattat                                                     >  2:1083355/1‑100 (MQ=255)
                                            cccACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATc                                                    <  1:1372307/100‑1 (MQ=255)
                                             ccACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCg                                                   <  1:1780678/100‑1 (MQ=255)
                                                 ttACTGTACTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGACATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAg                                               <  1:1070999/100‑1 (MQ=255)
                                                    cTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTg                                            <  2:919554/100‑1 (MQ=255)
                                                    cTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGGTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAACGCTg                                            >  1:587491/1‑100 (MQ=255)
                                                     atctctcTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCATGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGc                                           <  2:3572146/97‑1 (MQ=255)
                                                      gatctcTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGAGTGTTCTGCGTCCAGTCGGGAATATTATCGAAAGCTGCa                                          >  1:582129/3‑100 (MQ=255)
                                                      gTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCa                                          <  1:208389/100‑1 (MQ=255)
                                                       ttCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGCTTGTTCTGCGTCCAGTCGCGACTATTATCGAAAGCTGCAc                                         >  2:1875407/1‑100 (MQ=255)
                                                       ttCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCAc                                         <  2:736941/100‑1 (MQ=255)
                                                         ctctTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACgg                                       <  1:1060435/100‑1 (MQ=255)
                                                         ctctTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACgg                                       >  2:2876796/1‑100 (MQ=255)
                                                         ctctTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACgg                                       <  2:2788767/100‑1 (MQ=255)
                                                         ctctTTCTTGCCCCACAGCCAGGCGGGCAAACGTACTTTTCAAGAACCCTGAGAGAGATTGTGCTGCGTCCAGTGGAGAATATTATCGAAAGCTGCACgg                                       >  1:3794528/1‑100 (MQ=255)
                                                           ctttcttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTc                                     >  1:3643372/1‑100 (MQ=255)
                                                           ctttcttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTc                                     >  2:3662726/1‑100 (MQ=255)
                                                           ctttcttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTc                                     >  2:1121944/1‑100 (MQ=255)
                                                           ctttcttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTc                                     >  2:191398/1‑100 (MQ=255)
                                                           ctttcttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTACGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTc                                     >  1:3844528/1‑100 (MQ=255)
                                                            tttcttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACTGTCg                                    <  2:3381028/100‑1 (MQ=255)
                                                             ttcttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCgg                                   >  2:724716/1‑99 (MQ=255)
                                                              tcttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGGGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGtt                                  <  2:324055/100‑1 (MQ=255)
                                                               cttGCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTa                                 <  2:2684800/100‑1 (MQ=255)
                                                                  gCCCCACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAgcc                              >  2:2179047/1‑100 (MQ=255)
                                                                   ccccACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAgccg                             <  1:3448986/100‑1 (MQ=255)
                                                                    cccACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAgccgc                            >  1:3103155/1‑100 (MQ=255)
                                                                     ccACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAATCTGCACGGTCGTTAgccgcc                           <  1:3188878/100‑1 (MQ=255)
                                                                     ccACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAgccgcc                           >  1:74745/1‑100 (MQ=255)
                                                                     ccACAGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAgccgcc                           <  1:2114486/100‑1 (MQ=255)
                                                                       acaGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAg                         <  1:3317304/100‑1 (MQ=255)
                                                                        caGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAgg                        <  1:3853141/100‑1 (MQ=255)
                                                                        caGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAgg                        >  1:630712/1‑100 (MQ=255)
                                                                         gccctggCGGGCAAACGTTCTTTTCAAGAACGCTGAGCGGGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAggg                       <  2:2775957/95‑1 (MQ=255)
                                                                         aGCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAggg                       <  2:3421775/100‑1 (MQ=255)
                                                                          gCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTTCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGc                      >  2:1672794/1‑100 (MQ=255)
                                                                          gCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGc                      <  1:2661496/100‑1 (MQ=255)
                                                                          gCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGc                      >  2:1069149/1‑100 (MQ=255)
                                                                          gCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGc                      >  2:3783183/1‑100 (MQ=255)
                                                                          gCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGGCCAGTCGCGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGc                      >  1:3004593/1‑100 (MQ=255)
                                                                          gCCAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGACTGTTCTGCGTCCAGTCGAGAATATTATCGAAGGCTGCACGGTCGTTAGCCGCCAGGGc                      >  2:3801986/1‑100 (MQ=255)
                                                                            cAGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCag                    >  1:3887817/1‑100 (MQ=255)
                                                                             aGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCaga                   <  2:1426536/100‑1 (MQ=255)
                                                                             aGGCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCaga                   <  2:1119205/100‑1 (MQ=255)
                                                                              ggCGGGTAAACGTTCTTTTCAAGAACCCGGAGCGAGATTGTTCTGCGTCCACTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCagag                  <  2:2981998/100‑1 (MQ=255)
                                                                              ggCGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCagag                  <  2:2316785/100‑1 (MQ=255)
                                                                               gCGGGCAAACGTTCCTTTCAAGAACCCTGAGCGCGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAgaga                 <  2:2405323/100‑1 (MQ=255)
                                                                                cGGGCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGaa                >  1:1833337/1‑100 (MQ=255)
                                                                                 gggCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAAc               <  2:3914629/100‑1 (MQ=255)
                                                                                 gggCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAAc               <  2:212609/100‑1 (MQ=255)
                                                                                  ggCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACa              >  1:2835797/1‑100 (MQ=255)
                                                                                  ggCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACa              >  1:42816/1‑100 (MQ=255)
                                                                                   gCAAACGTTCTTTTCAAGTACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACAt             <  1:960094/100‑1 (MQ=255)
                                                                                   gCAAACGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACAt             <  2:2039472/100‑1 (MQ=255)
                                                                                     cagggTTCTTGAAAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACATGc           >  2:1132243/5‑100 (MQ=255)
                                                                                     aaaCGTTCTTTTCAAGGACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACATGc           <  1:3682052/100‑1 (MQ=255)
                                                                                     aaaCGTTCTTTTCAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACATGc           <  2:1362152/100‑1 (MQ=255)
                                                                                       gggTTCTTGAAAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGGATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACATGcgg         >  2:3981602/3‑98 (MQ=255)
                                                                                           tCTTTTCAAGAACCCTGAGCGAGATTGTTCCGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACATGCCGTAAc     <  1:2058253/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACAGCCAGGCGGGCAAACGTTCTTTTAAAGAACCCTGAGCGAGATTGTTCTGCGTCCAGTCGAGAATATTATCGAAAGCTGCACGGTCGTTAGCCGCCAGGGCAGAGAACATGCCGTAACTTT  >  NC_000913/3689899‑3690092

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