Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,051,650 A→G intergenic (‑300/‑15) b1980 ← / → shiA hypothetical protein/putative transport protein, shikimate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,051,6500AG83.1% 507.0 / 93.2 207intergenic (‑300/‑15)b1980/shiAhypothetical protein/putative transport protein, shikimate
Reads supporting (aligned to +/- strand):  ref base A (22/13);  new base G (92/80);  total (114/93)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 209 total aligned reads displayed.

AGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAG  >  NC_000913/2051558‑2051740
                                                                                            |                                                                                          
  aTGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCa                                                                                          >  1:4058108/1‑93 (MQ=255)
  aTGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACa                                                                                          <  1:1924617/93‑1 (MQ=255)
   tGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCaa                                                                                         <  2:225018/93‑1 (MQ=255)
    gcgcTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCaaa                                                                                        >  2:6319208/1‑93 (MQ=255)
      gcTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACa                                                                                      <  2:3360650/93‑1 (MQ=255)
      gcTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACa                                                                                      <  1:5609112/93‑1 (MQ=255)
       ctctGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAg                                                                                     >  2:2151000/1‑93 (MQ=255)
       ctctGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAg                                                                                     <  2:2066295/93‑1 (MQ=255)
        tctGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAgg                                                                                    <  1:6347441/93‑1 (MQ=255)
        tctGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAga                                                                                    >  2:2417413/1‑92 (MQ=255)
         ctGTCGAGGGAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCACACTAATTATAAAAAGACAGTGCAAACAGGa                                                                                   <  2:3734245/93‑1 (MQ=255)
         ctGTAGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGACAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGa                                                                                   <  2:3452632/93‑1 (MQ=255)
           gTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAATCAGGATc                                                                                 <  1:2262810/93‑1 (MQ=255)
            tCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCa                                                                                <  2:2992467/93‑1 (MQ=255)
            tCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCa                                                                                >  1:4355712/1‑93 (MQ=255)
             cGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCAt                                                                               <  1:4725531/93‑1 (MQ=255)
              gAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCAtt                                                                              >  2:4251014/1‑93 (MQ=255)
              gAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCAtt                                                                              <  2:6565391/93‑1 (MQ=255)
               aGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTa                                                                             >  2:2595079/1‑93 (MQ=255)
               aGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTa                                                                             >  2:2799733/1‑93 (MQ=255)
               aGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTa                                                                             >  1:6610515/1‑93 (MQ=255)
                ggTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTAt                                                                            >  2:3844577/1‑93 (MQ=255)
                ggTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTAt                                                                            <  1:2353850/93‑1 (MQ=255)
                 gTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATg                                                                           >  2:6055204/1‑93 (MQ=255)
                  tAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATgg                                                                          <  2:4487370/93‑1 (MQ=255)
                    gTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGAc                                                                        >  2:6171567/1‑93 (MQ=255)
                     tctcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACt                                                                       >  1:4943219/1‑93 (MQ=255)
                     tctcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACt                                                                       >  2:2529962/1‑93 (MQ=255)
                     tctcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACt                                                                       >  2:1517602/1‑93 (MQ=255)
                       tcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTcc                                                                     >  1:5809995/1‑93 (MQ=255)
                       tcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTcc                                                                     <  2:1135514/93‑1 (MQ=255)
                       tcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTcc                                                                     >  2:642686/1‑93 (MQ=255)
                         aTAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCAc                                                                   >  1:1067661/1‑93 (MQ=255)
                          tAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACg                                                                  >  2:4582735/1‑93 (MQ=255)
                             ggTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTc                                                               <  2:5869723/93‑1 (MQ=255)
                             ggTTCGTTTATAGATCGACGGCAATGGGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTc                                                               >  1:2084059/1‑93 (MQ=255)
                              gTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCa                                                              <  2:4187427/93‑1 (MQ=255)
                                   tttATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTcc                                                         <  2:2682153/93‑1 (MQ=255)
                                    ttATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCa                                                        <  1:1015189/93‑1 (MQ=255)
                                      ataGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACt                                                      >  1:2570682/1‑93 (MQ=255)
                                          aTCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTc                                                  >  1:2595118/1‑93 (MQ=255)
                                           tCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTcc                                                 >  2:1860719/1‑93 (MQ=255)
                                            cGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTccc                                                <  1:1955640/93‑1 (MQ=255)
                                              aCGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGa                                              <  2:6283572/93‑1 (MQ=255)
                                              aCGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGa                                              <  1:6466120/93‑1 (MQ=255)
                                               cGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGAt                                             >  1:1062939/1‑93 (MQ=255)
                                                ggCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATg                                            >  1:1282450/1‑93 (MQ=255)
                                                ggCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATg                                            <  1:3583626/93‑1 (MQ=255)
                                                 gCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGa                                           >  2:3282194/1‑93 (MQ=255)
                                                  cAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGaa                                          <  1:1806346/93‑1 (MQ=255)
                                                   aaTGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAg                                         <  1:734546/93‑1 (MQ=255)
                                                   aaTGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCAACGCTCATCTCCACTCGTCCCGATGAAg                                         >  1:3956991/1‑93 (MQ=255)
                                                     tgtgAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCATACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAgcg                                       >  2:3952516/1‑91 (MQ=255)
                                                     tgtgAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAggg                                       >  2:1337577/1‑93 (MQ=255)
                                                     tgtgAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAggg                                       >  1:4919539/1‑93 (MQ=255)
                                                      gtgAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGa                                      <  2:760138/93‑1 (MQ=255)
                                                      gtgAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGa                                      >  1:5504803/1‑93 (MQ=255)
                                                      gtgAGTTACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGa                                      <  1:6266780/93‑1 (MQ=255)
                                                        gAGTTACCTTTTCCATACTAATTATAAAAAGGCAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACg                                    >  2:1441476/1‑93 (MQ=255)
                                                          gTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCt                                  >  1:1427203/1‑93 (MQ=255)
                                                           ttACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCtt                                 >  1:1661881/1‑93 (MQ=255)
                                                           ttACCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCtt                                 <  1:251452/93‑1 (MQ=255)
                                                             aCCTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTc                               >  2:3153999/1‑93 (MQ=255)
                                                              ccTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCg                              >  2:1394038/1‑93 (MQ=255)
                                                               cTTTTCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGt                             >  2:2947318/1‑93 (MQ=255)
                                                               cTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGt                             >  1:3383457/1‑93 (MQ=255)
                                                                 tttCCATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTa                           <  1:3355138/93‑1 (MQ=255)
                                                                    ccATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGt                        >  2:3634789/1‑93 (MQ=255)
                                                                    ccATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGt                        <  1:3269919/93‑1 (MQ=255)
                                                                    ccATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGt                        >  1:2159106/1‑93 (MQ=255)
                                                                     cATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTc                       >  1:2793534/1‑93 (MQ=255)
                                                                     cATACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTc                       >  2:824843/1‑93 (MQ=255)
                                                                     cATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTc                       <  1:2897660/93‑1 (MQ=255)
                                                                      aTACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTcg                      <  1:1779122/93‑1 (MQ=255)
                                                                       tACTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTcgc                     >  2:489034/1‑93 (MQ=255)
                                                                        aCTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTcgcg                    >  1:1521316/1‑93 (MQ=255)
                                                                         cTAATTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCgcgc                   >  1:5318236/1‑93 (MQ=255)
                                                                           aaTTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGccc                 <  2:744325/93‑1 (MQ=255)
                                                                           aaTTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGccc                 <  2:4454480/93‑1 (MQ=255)
                                                                            aTTATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCCCTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCcg                <  1:5674631/93‑1 (MQ=255)
                                                                             ttATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCcga               >  1:3082045/1‑93 (MQ=255)
                                                                             ttATAAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCcga               <  2:4110023/93‑1 (MQ=255)
                                                                              tataAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCcgac              >  1:4023176/1‑93 (MQ=255)
                                                                              tataAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCcgac              >  1:3181292/1‑93 (MQ=255)
                                                                               ataAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCcgacg             >  2:1479134/1‑93 (MQ=255)
                                                                               ataAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCcgacg             >  1:5255017/1‑93 (MQ=255)
                                                                               ataAAAAGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCcgacg             >  2:6239198/1‑93 (MQ=255)
                                                                                 aaaaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAg           <  1:4429187/93‑1 (MQ=255)
                                                                                 aaaaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAg           <  1:6136062/93‑1 (MQ=255)
                                                                                  aaaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGTGc          <  1:4226571/93‑1 (MQ=255)
                                                                                  aaaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGc          >  1:614750/1‑93 (MQ=255)
                                                                                   aaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCt         <  1:1337577/93‑1 (MQ=255)
                                                                                   aaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCt         >  1:3985972/1‑93 (MQ=255)
                                                                                   aaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCt         >  1:4782489/1‑93 (MQ=255)
                                                                                   aaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCt         <  2:2178178/93‑1 (MQ=255)
                                                                                   aaaGACAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCt         <  1:3818168/93‑1 (MQ=255)
                                                                                   aaaGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCt         <  1:4647933/93‑1 (MQ=255)
                                                                                       aCAGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTGAGTCGCGCCCGACGAGCTGCGt     >  2:6445784/1‑93 (MQ=255)
                                                                                         aGTGCAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTa   <  2:167109/93‑1 (MQ=255)
                                                                                          gTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAg  >  2:3893568/1‑93 (MQ=255)
                                                                                            |                                                                                          
AGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAG  >  NC_000913/2051558‑2051740

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