Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,901,321 T→G A3A (GCA→GCC bglB ← phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,901,3210TG99.6% 787.4 / NA 248A3A (GCA→GCCbglBphospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base G (129/117);  minor base A (1/0);  total (131/117)
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

Only 100 of 251 total aligned reads displayed.

TGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAC  >  NC_000913/3901231‑3901410
                                                                                          |                                                                                         
tGAGTTCGAGATCCCTTCGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGc                                                                                         <  2:3075750/93‑1 (MQ=255)
  aGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGGTTCTGGAAAGGCtt                                                                                       >  1:3828800/1‑93 (MQ=255)
    gTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTc                                                                                     <  1:4101799/93‑1 (MQ=255)
     tCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCa                                                                                    <  1:3001855/93‑1 (MQ=255)
       gagaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATa                                                                                  >  1:1140777/1‑93 (MQ=255)
       gagaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATa                                                                                  >  1:3693037/1‑93 (MQ=255)
       gagaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATa                                                                                  <  1:3227788/93‑1 (MQ=255)
       gagaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATa                                                                                  >  1:2713458/1‑93 (MQ=255)
       gagaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGATTTCATa                                                                                  >  1:3575618/1‑93 (MQ=255)
        agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATaa                                                                                 >  1:1952161/1‑93 (MQ=255)
        agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAACAAATGTTTCTGGAAAGGCTTTCATat                                                                                 >  1:1690460/1‑92 (MQ=255)
        acaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATaa                                                                                 <  1:2439729/91‑1 (MQ=255)
         gaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAt                                                                                >  2:93885/1‑93 (MQ=255)
         gaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAt                                                                                <  1:1009515/93‑1 (MQ=255)
         gaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAt                                                                                >  2:2574728/1‑93 (MQ=255)
          aTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAtt                                                                               >  2:2108411/1‑93 (MQ=255)
           tCCCTTTGCCATCTTCCTGGCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTa                                                                              >  2:4584573/1‑93 (MQ=255)
           tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTa                                                                              >  2:1639233/1‑93 (MQ=255)
            cccTTTGCCATCTTCCTGCCAGGCACCTTCCTCCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTaa                                                                             <  1:1390686/93‑1 (MQ=255)
            cccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTaa                                                                             <  1:1760144/93‑1 (MQ=255)
            cccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTaa                                                                             >  1:3131362/1‑93 (MQ=255)
             ccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAAc                                                                            <  2:2627042/93‑1 (MQ=255)
             ccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAAc                                                                            <  2:3349558/93‑1 (MQ=255)
              cTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACt                                                                           >  2:3885774/1‑93 (MQ=255)
               tttGCCATTTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAACTGTTTCTGGAAAGGCTTTTATAATTAACtt                                                                          >  2:2885733/1‑92 (MQ=255)
               tttGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTc                                                                          <  2:194259/93‑1 (MQ=255)
                  gCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCtt                                                                       <  2:4528739/93‑1 (MQ=255)
                   ccATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCttt                                                                      >  2:1441225/1‑93 (MQ=255)
                    cATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCtttt                                                                     >  1:3996874/1‑93 (MQ=255)
                     aTCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTa                                                                    >  1:4648197/1‑93 (MQ=255)
                       cTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATc                                                                  <  2:603015/93‑1 (MQ=255)
                       cTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATc                                                                  <  2:1604540/93‑1 (MQ=255)
                       cTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATc                                                                  <  2:962375/93‑1 (MQ=255)
                       cTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATc                                                                  <  1:1881541/93‑1 (MQ=255)
                         tCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGt                                                                >  2:352372/1‑93 (MQ=255)
                          ccTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGtt                                                               <  2:4575671/93‑1 (MQ=255)
                           cTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTa                                                              <  2:2452576/93‑1 (MQ=255)
                            tGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAg                                                             <  1:3419739/93‑1 (MQ=255)
                            tGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAg                                                             <  1:3961425/93‑1 (MQ=255)
                             gCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGc                                                            <  1:977008/93‑1 (MQ=255)
                             gCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGc                                                            <  1:1581635/93‑1 (MQ=255)
                                 ggCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAat                                                        >  1:1823456/1‑93 (MQ=255)
                                    aCCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatgat                                                     <  1:2046124/93‑1 (MQ=255)
                                     ccTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatgatg                                                    <  1:2604328/93‑1 (MQ=255)
                                     ccTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatgatg                                                    <  2:1324685/93‑1 (MQ=255)
                                      cTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATgg                                                   <  2:538393/93‑1 (MQ=255)
                                          cACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATaa                                               <  1:4003243/93‑1 (MQ=255)
                                          cACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATaa                                               >  1:78005/1‑93 (MQ=255)
                                            ccTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACa                                             <  2:3803119/93‑1 (MQ=255)
                                             cTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAg                                            >  1:606710/1‑93 (MQ=255)
                                              tGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGc                                           >  1:4464964/1‑93 (MQ=255)
                                              tGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGc                                           <  2:1019154/93‑1 (MQ=255)
                                               ggttGGGAGTTGGTGCCCCCCCCCAAAGAAATTGTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCg                                          <  2:361297/91‑1 (MQ=255)
                                               gATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCg                                          >  2:945766/1‑93 (MQ=255)
                                                 ttGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGt                                        >  2:2926534/1‑93 (MQ=255)
                                                    gCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGttta                                     >  2:4395927/1‑91 (MQ=255)
                                                    gCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCa                                     >  1:712817/1‑93 (MQ=255)
                                                     cAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCAc                                    <  1:558671/93‑1 (MQ=255)
                                                      aGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCAcc                                   >  2:1907843/1‑93 (MQ=255)
                                                      aGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCAcc                                   <  2:737094/93‑1 (MQ=255)
                                                       gCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCt                                  <  2:2819071/93‑1 (MQ=255)
                                                        cTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTg                                 <  2:3573614/93‑1 (MQ=255)
                                                           ttGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcgc                              <  2:638953/93‑1 (MQ=255)
                                                              cgccgcGCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTa                           <  1:2474782/93‑1 (MQ=255)
                                                              cgcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTGTCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTa                           >  2:805142/1‑93 (MQ=255)
                                                               gcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTAt                          >  2:294018/1‑93 (MQ=255)
                                                               gcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTAt                          >  2:3310343/1‑93 (MQ=255)
                                                               gcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTAt                          >  1:3090300/1‑93 (MQ=255)
                                                               gcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTAt                          >  2:3722990/1‑93 (MQ=255)
                                                                ccACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATc                         >  1:2913714/1‑93 (MQ=255)
                                                                  aCCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTg                       >  1:2299287/1‑93 (MQ=255)
                                                                  aCCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTg                       >  1:2324819/1‑93 (MQ=255)
                                                                  aCCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTg                       >  2:4347467/1‑93 (MQ=255)
                                                                  aCCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTg                       >  2:792977/1‑93 (MQ=255)
                                                                     ccAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTgcgc                    >  2:2525139/1‑93 (MQ=255)
                                                                     ccAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTgcgc                    >  1:1897827/1‑93 (MQ=255)
                                                                     ccAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTgcgc                    >  2:1298197/1‑93 (MQ=255)
                                                                      cAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgcc                   >  1:2118997/1‑93 (MQ=255)
                                                                       aaaGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccg                  >  2:2208481/1‑93 (MQ=255)
                                                                        aaGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgt                 <  2:1508617/93‑1 (MQ=255)
                                                                         aGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtg                >  1:2234766/1‑93 (MQ=255)
                                                                         aGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtg                >  2:3986143/1‑93 (MQ=255)
                                                                          gAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtgc               >  1:1288972/1‑93 (MQ=255)
                                                                          gAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtgc               >  2:4178691/1‑93 (MQ=255)
                                                                            aaTGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtgccg             >  2:3444579/1‑93 (MQ=255)
                                                                            aaTGTTCCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtgccg             <  1:2155180/93‑1 (MQ=255)
                                                                              tGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGccgtgccgtg           <  1:3750016/93‑1 (MQ=255)
                                                                              tGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGccgtgccgtg           >  1:1643095/1‑93 (MQ=255)
                                                                              tGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGccgtgccgtg           <  2:3159735/93‑1 (MQ=255)
                                                                               gTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTgg          >  1:46922/1‑93 (MQ=255)
                                                                                tttCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTggg         <  1:2837411/93‑1 (MQ=255)
                                                                                 ttCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTgggg        >  1:3145119/1‑93 (MQ=255)
                                                                                  tCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGt       <  2:122864/93‑1 (MQ=255)
                                                                                    tGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTaa     >  2:3401715/1‑93 (MQ=255)
                                                                                    tGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTaa     >  2:3423264/1‑93 (MQ=255)
                                                                                    tGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTaa     >  1:4517146/1‑93 (MQ=255)
                                                                                      gAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATa   >  2:4481811/1‑93 (MQ=255)
                                                                                      gAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATa   <  1:3263282/93‑1 (MQ=255)
                                                                                       aaaGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAc  <  1:4241376/93‑1 (MQ=255)
                                                                                       aaaGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAc  <  2:4301750/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAC  >  NC_000913/3901231‑3901410

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