Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,431,142 T→C H154R (CAC→CGC)  ytfG ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,431,1420TC99.7% 933.2 / ‑5.0 295H154R (CAC→CGC) ytfGputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (144/150);  total (144/151)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 319 total aligned reads displayed.

TCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGC  >  NC_000913/4431052‑4431231
                                                                                          |                                                                                         
tCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCgcgt                                                                                         >  1:1340129/1‑93 (MQ=255)
tCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCgcgt                                                                                         >  1:9172438/1‑93 (MQ=255)
    tAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCa                                                                                     <  2:865755/93‑1 (MQ=255)
    tAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCa                                                                                     >  1:834970/1‑93 (MQ=255)
     aaTCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAg                                                                                    >  1:7763414/1‑93 (MQ=255)
     aaTCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAg                                                                                    >  1:5865945/1‑93 (MQ=255)
      aTCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGt                                                                                   <  2:9653461/93‑1 (MQ=255)
        cACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgc                                                                                 >  2:234557/1‑93 (MQ=255)
         ccgcgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgct                                                                                <  2:4766356/92‑1 (MQ=255)
          cgcgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgctg                                                                               <  1:239518/93‑1 (MQ=255)
          cgcgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgctg                                                                               >  1:10017976/1‑93 (MQ=255)
          cgcgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgctg                                                                               >  2:515055/1‑93 (MQ=255)
           gggTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgctgc                                                                              <  1:5202654/91‑1 (MQ=255)
           gcgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgctgc                                                                              <  1:6128011/93‑1 (MQ=255)
           gcgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgctgc                                                                              >  2:565507/1‑93 (MQ=255)
           ccgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgctgc                                                                              <  1:2624995/92‑1 (MQ=255)
            cgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGcc                                                                             >  2:9902363/1‑93 (MQ=255)
             gTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCg                                                                            >  2:6059457/1‑93 (MQ=255)
              tGCCGCAGCGGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCgg                                                                           <  1:1519551/93‑1 (MQ=255)
               gCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCggg                                                                          >  2:210536/1‑93 (MQ=255)
               gCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCggg                                                                          <  1:910521/93‑1 (MQ=255)
                 cGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgc                                                                        >  2:3786875/1‑93 (MQ=255)
                  gCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgcg                                                                       >  2:2525246/1‑93 (MQ=255)
                  gCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgcg                                                                       >  1:4815235/1‑93 (MQ=255)
                  gCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgcg                                                                       >  1:8110395/1‑93 (MQ=255)
                   cAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgcgc                                                                      >  2:4639566/1‑93 (MQ=255)
                   cAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgcgc                                                                      >  1:2784725/1‑93 (MQ=255)
                    aGCTGCCGCATAAGCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCt                                                                     <  1:4118961/93‑1 (MQ=255)
                      cTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTgg                                                                   >  2:8913411/1‑93 (MQ=255)
                      cTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTgg                                                                   <  1:6411354/93‑1 (MQ=255)
                        gCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCg                                                                 >  1:9553223/1‑93 (MQ=255)
                         ccGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGa                                                                <  2:1138126/93‑1 (MQ=255)
                         ccGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGa                                                                >  2:9130610/1‑93 (MQ=255)
                         ccGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGa                                                                <  1:2078488/93‑1 (MQ=255)
                             aTAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTa                                                            >  2:2131893/1‑93 (MQ=255)
                              tAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAg                                                           >  1:9776775/1‑93 (MQ=255)
                              tAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAg                                                           >  1:4935148/1‑93 (MQ=255)
                               aaTCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGt                                                          >  2:1304389/1‑93 (MQ=255)
                                aTCTGCCCGCGTTGCTGAGGCGATTTTGCCATAGCCCGCCGCACCGATAAATAGGCCGCGTTCCAGTGCTGCCGGGGCGGCGGCGGGGTGGtt                                                         >  1:2958981/1‑93 (MQ=255)
                                 tCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGttt                                                        >  2:303626/1‑93 (MQ=255)
                                  cTGCCCGCGGTGCGGAGGCGGTTGTGCCATCGCCCGCCGCAGCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGtttt                                                       <  1:8647874/93‑1 (MQ=255)
                                      ccGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCt                                                   <  2:3283022/93‑1 (MQ=255)
                                       cgcgTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTg                                                  <  2:1623827/93‑1 (MQ=255)
                                       cgcgTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTg                                                  >  2:7864462/1‑93 (MQ=255)
                                       cgcgTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTg                                                  <  2:7464407/93‑1 (MQ=255)
                                       cgcgTTCCTGAGGCGATTTTGCCATGCCCCGCCGCACCTAGAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTg                                                  <  2:9549796/93‑1 (MQ=255)
                                        gcgTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGt                                                 >  2:7288345/1‑93 (MQ=255)
                                         cgTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTa                                                >  1:5539836/1‑93 (MQ=255)
                                          gTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAc                                               <  1:2440725/93‑1 (MQ=255)
                                          gTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAc                                               >  2:9439593/1‑93 (MQ=255)
                                           ttGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAcc                                              <  2:356836/93‑1 (MQ=255)
                                             gCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGGAGTTTTCGCTGTACCAg                                            >  1:9917375/1‑93 (MQ=255)
                                              cTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGc                                           <  2:4582122/93‑1 (MQ=255)
                                               tGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGcc                                          <  2:8459962/93‑1 (MQ=255)
                                               tGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGcc                                          >  1:1974106/1‑93 (MQ=255)
                                                gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCg                                         >  2:5066869/1‑93 (MQ=255)
                                                gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCg                                         >  1:5354673/1‑93 (MQ=255)
                                                gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCg                                         >  1:3722295/1‑93 (MQ=255)
                                                gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCg                                         >  2:9790130/1‑93 (MQ=255)
                                                 gggCGATTTTGCCATCGCCCGCCGCACCGATAAATCCGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGt                                        <  1:7340163/92‑1 (MQ=255)
                                                 aGGCGCTTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGGCCTGCCGGGGCGCTGGCGAGGAGGTTTTGGATGTGCCAGCCGt                                        >  1:7957626/1‑93 (MQ=255)
                                                 aGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGTCGAGGTAGTTTTCGCTGTACCAGCCGt                                        <  1:6831717/93‑1 (MQ=255)
                                                  ggCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGtt                                       >  1:506341/1‑93 (MQ=255)
                                                   gCGTTTTTGCCACCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTg                                      <  2:7999440/93‑1 (MQ=255)
                                                   gCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTg                                      >  1:5227586/1‑93 (MQ=255)
                                                    cgtttttGCCTCCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgc                                     <  2:3835019/90‑1 (MQ=255)
                                                      aTTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgcgc                                   <  1:3941085/93‑1 (MQ=255)
                                                       ttttGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgca                                  >  1:8753486/1‑93 (MQ=255)
                                                         ttGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagc                                >  2:3971032/1‑93 (MQ=255)
                                                         ttGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagc                                >  1:1717401/1‑93 (MQ=255)
                                                          tGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagca                               <  1:3433085/93‑1 (MQ=255)
                                                          tGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagca                               <  2:5493294/93‑1 (MQ=255)
                                                           gCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGcagcag                              >  1:6671979/1‑93 (MQ=255)
                                                           gCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGcagcag                              >  2:10025416/1‑93 (MQ=255)
                                                             cATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAggg                            <  1:6631681/93‑1 (MQ=255)
                                                               tCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGgtg                          <  1:193197/93‑1 (MQ=255)
                                                                 gcccgccGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTa                        >  1:1127078/1‑93 (MQ=255)
                                                                  cccgccGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTaa                       >  2:6620776/1‑93 (MQ=255)
                                                                  cccgccGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTaa                       >  1:6271728/1‑93 (MQ=255)
                                                                   ccgccgCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTaaa                      <  2:7702218/93‑1 (MQ=255)
                                                                   ccgcccCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTaaa                      <  2:3176721/93‑1 (MQ=255)
                                                                     gccgcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACg                    >  2:2245054/1‑93 (MQ=255)
                                                                      ccgcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGa                   <  2:6923989/93‑1 (MQ=255)
                                                                       cgcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGAt                  <  1:5806962/93‑1 (MQ=255)
                                                                         cACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGc                <  2:4008598/93‑1 (MQ=255)
                                                                          aCCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGcc               >  2:8184610/1‑93 (MQ=255)
                                                                          aCCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGcc               >  2:2077938/1‑93 (MQ=255)
                                                                           ccGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCa              >  2:8757175/1‑93 (MQ=255)
                                                                           ccGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCa              <  2:1555024/93‑1 (MQ=255)
                                                                            cGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAg             <  1:3268407/93‑1 (MQ=255)
                                                                             caTAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGa            <  2:3426123/92‑1 (MQ=255)
                                                                              aTAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGaa           <  2:8838504/93‑1 (MQ=255)
                                                                              aTAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGaa           <  2:6250358/93‑1 (MQ=255)
                                                                               tAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAAt          <  2:7062906/93‑1 (MQ=255)
                                                                                aaaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATc         <  2:6756119/93‑1 (MQ=255)
                                                                                aaaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATc         <  1:1247722/93‑1 (MQ=255)
                                                                                aaaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATc         <  1:7138179/93‑1 (MQ=255)
                                                                                aaaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATc         <  2:4072939/93‑1 (MQ=255)
                                                                                  aTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcag       >  2:5326303/1‑93 (MQ=255)
                                                                                  aTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcag       >  1:4379258/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGC  >  NC_000913/4431052‑4431231

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