Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 160,649 A→G intergenic (‑45/+133) dksA ← / ← sfsA DksA/putative regulator for maltose metabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913160,6490AG99.4% 509.9 / ‑5.6 160intergenic (‑45/+133)dksA/sfsADksA/putative regulator for maltose metabolism
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (85/74);  total (86/74)
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 181 total aligned reads displayed.

GAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTT  >  NC_000913/160560‑160739
                                                                                         |                                                                                          
gagaATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGgaaa                                                                                         >  2:2402415/1‑93 (MQ=255)
gagaATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGgaaa                                                                                         >  1:2380056/1‑93 (MQ=255)
 agaATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGaaaaa                                                                                        >  2:787271/1‑93 (MQ=255)
  gaATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAt                                                                                       >  2:2092330/1‑93 (MQ=255)
    aTACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGa                                                                                     <  1:424660/93‑1 (MQ=255)
    aTACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGa                                                                                     <  1:4487248/93‑1 (MQ=255)
    aTACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGa                                                                                     <  2:1943851/93‑1 (MQ=255)
     tACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAg                                                                                    >  2:4517240/1‑93 (MQ=255)
     tACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAg                                                                                    <  2:2364179/93‑1 (MQ=255)
     tACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACGATCGATCCCCCTGTTCGGGGGGAAAATGAg                                                                                    >  2:3946486/1‑93 (MQ=255)
      aCTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAgg                                                                                   <  1:2171101/93‑1 (MQ=255)
      aCTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAgg                                                                                   <  2:849534/93‑1 (MQ=255)
      aCTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAgg                                                                                   <  2:3918843/93‑1 (MQ=255)
       cTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGc                                                                                  >  1:4741143/1‑93 (MQ=255)
       cTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGc                                                                                  >  2:2118332/1‑93 (MQ=255)
        tCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGcc                                                                                 <  2:1705656/93‑1 (MQ=255)
        tCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGcc                                                                                 <  1:1264585/93‑1 (MQ=255)
          aGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGc                                                                               >  2:1823513/1‑93 (MQ=255)
           gggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGTCGCt                                                                              >  2:2732730/1‑93 (MQ=255)
             gACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCtat                                                                            >  2:1325876/1‑93 (MQ=255)
               cGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATaa                                                                          >  2:3706013/1‑93 (MQ=255)
                gATGTTTTACTGTTTTGCCCTTCTTTCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATaaa                                                                         <  1:1614733/93‑1 (MQ=255)
                gATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCGCATGTTCGGGGGGAAAATGAGGCCGCTATaaa                                                                         >  1:745823/1‑93 (MQ=255)
                gATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATaaa                                                                         >  1:3777583/1‑93 (MQ=255)
                gATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATaaa                                                                         <  1:1036332/93‑1 (MQ=255)
                  tGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATa                                                                       <  1:3625407/93‑1 (MQ=255)
                  tGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATa                                                                       <  2:2755240/93‑1 (MQ=255)
                   gTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAg                                                                      <  1:1307036/93‑1 (MQ=255)
                   gTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAg                                                                      >  2:590095/1‑93 (MQ=255)
                    ttttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGc                                                                     >  2:3368142/1‑93 (MQ=255)
                     tttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCa                                                                    <  2:2912137/93‑1 (MQ=255)
                      ttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAg                                                                   <  1:3878639/93‑1 (MQ=255)
                       tACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGa                                                                  <  1:748126/93‑1 (MQ=255)
                       tACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGa                                                                  <  1:2598831/93‑1 (MQ=255)
                        tcGGTTTTGCCCTTCTTTCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGAt                                                                 <  1:1030769/92‑1 (MQ=255)
                         cGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATg                                                                >  1:3352958/1‑93 (MQ=255)
                          ggTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCCTGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGc                                                               >  1:3217174/1‑93 (MQ=255)
                          ggTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGc                                                               >  2:2948967/1‑93 (MQ=255)
                          ggTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGc                                                               >  2:3477883/1‑93 (MQ=255)
                            ttttGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCtt                                                             >  2:261975/1‑93 (MQ=255)
                             tttGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCttt                                                            <  1:718111/93‑1 (MQ=255)
                              ttGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCtttt                                                           <  1:1782921/93‑1 (MQ=255)
                               tGCCCTGCTTGCATGTTGATTCTCCTTATCGCGCACTATCGATCCCGATGTTCGGGGCGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTc                                                          <  2:2810963/93‑1 (MQ=255)
                                gCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTcc                                                         >  1:2905745/1‑93 (MQ=255)
                                 cccTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCg                                                        <  2:152983/93‑1 (MQ=255)
                                 cccTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCg                                                        >  1:3446815/1‑93 (MQ=255)
                                  ccTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCgg                                                       >  1:2495347/1‑93 (MQ=255)
                                  ccTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCgg                                                       >  2:1695020/1‑93 (MQ=255)
                                   cttcttGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCggg                                                      >  1:2255532/1‑92 (MQ=255)
                                   cttcttGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGa                                                      <  2:4244323/93‑1 (MQ=255)
                                   cttcttGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGa                                                      >  2:2994717/1‑93 (MQ=255)
                                       ttGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCTGATAGc                                                  >  2:935934/1‑93 (MQ=255)
                                          cATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAAt                                               >  2:773371/1‑93 (MQ=255)
                                          cATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAAt                                               >  1:4344148/1‑93 (MQ=255)
                                           ctGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAAtt                                              <  1:3598876/92‑1 (MQ=255)
                                           aTGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGGTAGCAAtt                                              >  1:2543725/1‑93 (MQ=255)
                                           aTGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAAtt                                              >  1:2031796/1‑93 (MQ=255)
                                            tGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTa                                             >  1:1497543/1‑93 (MQ=255)
                                            tGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTa                                             <  2:4575985/93‑1 (MQ=255)
                                             gTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAGAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTAt                                            >  2:2689088/1‑93 (MQ=255)
                                             gTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTAt                                            >  2:3166632/1‑93 (MQ=255)
                                              ttGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATc                                           >  2:2036421/1‑93 (MQ=255)
                                               tGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCt                                          <  2:1085463/93‑1 (MQ=255)
                                               tGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCt                                          <  2:2331733/93‑1 (MQ=255)
                                                  ttCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTaaa                                       <  2:4431039/93‑1 (MQ=255)
                                                      ccTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTa                                   <  2:4289860/93‑1 (MQ=255)
                                                       cTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATGGCAATTATCTAAACGTaa                                  >  1:952577/1‑93 (MQ=255)
                                                       cTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTaa                                  <  1:2102525/93‑1 (MQ=255)
                                                        ttAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAac                                 >  1:4072175/1‑93 (MQ=255)
                                                        ttAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAac                                 >  2:3500810/1‑93 (MQ=255)
                                                         tAACACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAaca                                <  1:3038286/93‑1 (MQ=255)
                                                          aaCACGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAacac                               <  2:3074358/93‑1 (MQ=255)
                                                            cacGCACTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACtt                             <  1:3686934/93‑1 (MQ=255)
                                                                 aCTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                        >  2:4701627/1‑93 (MQ=255)
                                                                 aCTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                        >  1:2449627/1‑93 (MQ=255)
                                                                 aCTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                        <  2:519164/93‑1 (MQ=255)
                                                                 aCTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                        >  1:4386893/1‑93 (MQ=255)
                                                                  cTATCGATCCCCATGTTCGGGTGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAAc                       >  1:339235/1‑93 (MQ=255)
                                                                   tATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACt                      >  2:1185230/1‑93 (MQ=255)
                                                                    aTCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtg                     >  1:1402769/1‑93 (MQ=255)
                                                                    aTCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtg                     >  2:3633401/1‑93 (MQ=255)
                                                                     tCGATCCCCATGTTCGGGGGGCAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtgt                    >  1:904191/1‑93 (MQ=255)
                                                                     tCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtgt                    >  2:3957333/1‑93 (MQ=255)
                                                                       gATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtgtgt                  >  2:1442495/1‑93 (MQ=255)
                                                                        aTCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtgtgtg                 <  1:658980/93‑1 (MQ=255)
                                                                         tCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGa                <  1:4533120/93‑1 (MQ=255)
                                                                         tCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGa                >  2:1728613/1‑93 (MQ=255)
                                                                          ccccATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAg               >  1:640122/1‑93 (MQ=255)
                                                                           cccATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAgg              <  1:570054/93‑1 (MQ=255)
                                                                            ccATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGa             <  2:1234205/93‑1 (MQ=255)
                                                                            ccATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGa             >  2:1991891/1‑93 (MQ=255)
                                                                            ccATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGa             >  2:3405381/1‑93 (MQ=255)
                                                                            ccATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGa             >  2:507368/1‑93 (MQ=255)
                                                                               tGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGaaaa          >  1:2207369/1‑93 (MQ=255)
                                                                               tGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGaaaa          <  1:4725435/93‑1 (MQ=255)
                                                                                gTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAg         >  2:2382660/1‑93 (MQ=255)
                                                                                gTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAg         >  2:2867401/1‑93 (MQ=255)
                                                                                    ggggggAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTa     >  2:434262/1‑93 (MQ=255)
                                                                                      ggggAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                        >  2:2249055/1‑72 (MQ=255)
                                                                                      ggggAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTAtt   <  1:4153256/93‑1 (MQ=255)
                                                                                         |                                                                                          
GAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTT  >  NC_000913/160560‑160739

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