Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 374,060 G→A E323E (GAG→GAA)  mhpE → 4‑hydroxy‑2‑ketovalerate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913374,0600GA97.4% 537.5 / ‑0.1 156E323E (GAG→GAA) mhpE4‑hydroxy‑2‑ketovalerate aldolase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (77/75);  total (79/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

Only 100 of 175 total aligned reads displayed.

TCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGAC  >  NC_000913/373970‑374152
                                                                                          |                                                                                            
tCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAggaaga                                                                                            >  2:1444463/1‑93 (MQ=255)
tCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAggaaga                                                                                            >  1:1649066/1‑93 (MQ=255)
tCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAggaaga                                                                                            <  2:3496670/93‑1 (MQ=255)
tCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCAAggaaga                                                                                            >  1:706989/1‑93 (MQ=255)
 cACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGat                                                                                           >  2:1910774/1‑93 (MQ=255)
  aCTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGata                                                                                          <  2:1737857/93‑1 (MQ=255)
   cTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGatat                                                                                         <  1:17829/93‑1 (MQ=255)
   cTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGatat                                                                                         <  2:199545/93‑1 (MQ=255)
   cTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGatat                                                                                         >  1:2295752/1‑93 (MQ=255)
    tgtgAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATg                                                                                        <  1:655242/93‑1 (MQ=255)
      tgAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGAt                                                                                      <  1:3079736/93‑1 (MQ=255)
      tgAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGAt                                                                                      <  1:2086851/93‑1 (MQ=255)
       gAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATc                                                                                     <  1:3214825/93‑1 (MQ=255)
       gAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATc                                                                                     <  1:4065965/93‑1 (MQ=255)
       gAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATc                                                                                     <  1:2445371/93‑1 (MQ=255)
         aaCGGCGGCGGCGCGTTATGGCTTTAGTGCGGGGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGt                                                                                   >  2:3670567/1‑93 (MQ=255)
         aaCGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGt                                                                                   <  2:3016460/93‑1 (MQ=255)
         aaCGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGt                                                                                   <  2:2323992/93‑1 (MQ=255)
           cggcggcggCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTg                                                                                 >  1:3418432/1‑93 (MQ=255)
               ggcggcgCGTTATGGCTTATGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGt                                                                             <  1:4617707/93‑1 (MQ=255)
               ggcggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGt                                                                             >  2:2004149/1‑93 (MQ=255)
                gcggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTg                                                                            <  1:1382872/93‑1 (MQ=255)
                 cggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTgg                                                                           <  2:2487347/93‑1 (MQ=255)
                  ggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGgc                                                                          >  1:301248/1‑93 (MQ=255)
                  ggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTAACGTgcg                                                                          >  1:3719922/1‑91 (MQ=255)
                   gcgcGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGgcg                                                                         >  2:1515/1‑93 (MQ=255)
                    cgcgTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGgcgc                                                                        <  1:627823/93‑1 (MQ=255)
                     gcgTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCt                                                                       <  2:2438136/93‑1 (MQ=255)
                     gcgTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCt                                                                       <  1:3967883/93‑1 (MQ=255)
                      cgTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTg                                                                      <  2:682432/93‑1 (MQ=255)
                       gTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTgg                                                                     <  2:3104197/93‑1 (MQ=255)
                        ttATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGa                                                                    >  1:3862406/1‑93 (MQ=255)
                        ttATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGa                                                                    >  2:4115939/1‑93 (MQ=255)
                         tATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGAt                                                                   >  2:1656546/1‑93 (MQ=255)
                          aTGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATc                                                                  >  2:2137188/1‑93 (MQ=255)
                          aTGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATc                                                                  >  2:574096/1‑93 (MQ=255)
                           tGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCt                                                                 >  1:1898997/1‑93 (MQ=255)
                             gCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTgc                                                               <  1:3892792/93‑1 (MQ=255)
                             gCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTgc                                                               >  2:1741444/1‑93 (MQ=255)
                             gCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTgc                                                               >  1:4105109/1‑93 (MQ=255)
                               ttAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTgcgc                                                             >  2:4001588/1‑93 (MQ=255)
                                  aGTGCGTTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaac                                                          <  1:1285550/93‑1 (MQ=255)
                                  aGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaac                                                          <  1:184377/93‑1 (MQ=255)
                                     gCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaacaac                                                       >  1:554607/1‑93 (MQ=255)
                                       ggTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaacaacaa                                                     >  1:2667953/1‑93 (MQ=255)
                                         tGGATATTCTCGTTGAGCTGGGCAGACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaat                                                   <  1:3181183/93‑1 (MQ=255)
                                          ggATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaata                                                  >  1:492477/1‑93 (MQ=255)
                                            atatTCTCGTTGAGCTTGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaataat                                                >  2:2074046/1‑93 (MQ=255)
                                              atTCTCGTTGAGCTGGGCAAACGCCGGCTGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGGGCAACAACAAATAatga                                              >  2:3620033/1‑93 (MQ=255)
                                              atTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatga                                              >  1:791887/1‑93 (MQ=255)
                                               ttCTCGTTGAGCTGGGCAAACTCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatgat                                             <  2:2426567/93‑1 (MQ=255)
                                               ttCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaacaa                                                        <  2:1366498/82‑1 (MQ=255)
                                               ttCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatgat                                             <  2:1068051/93‑1 (MQ=255)
                                               ttCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatgat                                             <  2:3664337/93‑1 (MQ=255)
                                               ttCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatgat                                             >  2:1238517/1‑93 (MQ=255)
                                               ttCTCGTTGAGCTGGGCAAACGCCGGATGGCTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatgat                                             >  1:142184/1‑93 (MQ=255)
                                                    gTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTg                                        >  2:2654789/1‑93 (MQ=255)
                                                     ttGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGc                                       >  1:1045982/1‑93 (MQ=255)
                                                      tGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGcc                                      <  2:3204741/93‑1 (MQ=255)
                                                       gAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCg                                     >  2:514388/1‑93 (MQ=255)
                                                       gAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCg                                     >  2:2028300/1‑93 (MQ=255)
                                                         gCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCgag                                   >  2:3043714/1‑93 (MQ=255)
                                                           tGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGaa                                 >  1:259292/1‑93 (MQ=255)
                                                           tGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGaa                                 <  1:3217618/93‑1 (MQ=255)
                                                             ggCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACg                               <  1:2326105/93‑1 (MQ=255)
                                                             ggCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACg                               >  1:2426035/1‑93 (MQ=255)
                                                             ggCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACg                               <  2:3727219/93‑1 (MQ=255)
                                                               cAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTg                             >  2:4654382/1‑93 (MQ=255)
                                                                 aaCGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCa                           >  1:2949709/1‑93 (MQ=255)
                                                                    gCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAttt                        <  2:3815287/93‑1 (MQ=255)
                                                                     ccGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAtttt                       <  1:3419305/93‑1 (MQ=255)
                                                                     ccGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAtttt                       >  1:1879898/1‑93 (MQ=255)
                                                                       ggATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGt                     <  2:2013932/93‑1 (MQ=255)
                                                                       ggATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGt                     >  1:1363489/1‑93 (MQ=255)
                                                                       ggATGGTTGGCGGCCAGGAAGATATGATCGTTGAAGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGt                     >  2:186596/1‑93 (MQ=255)
                                                                         gtGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTAt                   <  1:377896/92‑1 (MQ=255)
                                                                         aTGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTAt                   <  1:3121567/93‑1 (MQ=255)
                                                                         aTGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTAt                   <  2:3224140/93‑1 (MQ=255)
                                                                          tGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATg                  <  1:2013395/93‑1 (MQ=255)
                                                                          tGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATg                  >  1:3881406/1‑93 (MQ=255)
                                                                          tGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATg                  >  1:4124566/1‑93 (MQ=255)
                                                                            gTTGTCAGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGcc                <  2:2090554/93‑1 (MQ=255)
                                                                            gTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGcc                >  1:2531360/1‑93 (MQ=255)
                                                                            gTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGcc                <  1:2918669/93‑1 (MQ=255)
                                                                              tGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCa              <  2:627131/93‑1 (MQ=255)
                                                                              tGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCa              >  1:3333709/1‑93 (MQ=255)
                                                                                   gcCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAAGGTGCATTTTGTATGCCCATCCCa         >  2:2606393/1‑93 (MQ=255)
                                                                                   gcCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCa         >  2:3572500/1‑93 (MQ=255)
                                                                                   gcCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCa         >  1:251219/1‑93 (MQ=255)
                                                                                   gcCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCa         <  2:2574336/93‑1 (MQ=255)
                                                                                    ccAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAg        >  2:4503494/1‑93 (MQ=255)
                                                                                    ccAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAg        >  2:4503493/1‑93 (MQ=255)
                                                                                     cAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGc       <  1:368869/93‑1 (MQ=255)
                                                                                     cAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGc       >  2:4087149/1‑93 (MQ=255)
                                                                                     cAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGc       >  2:2944814/1‑93 (MQ=255)
                                                                                      aggaagATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCa      >  1:2895151/1‑93 (MQ=255)
                                                                                       ggaagaTATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCAc     >  2:4184489/1‑93 (MQ=255)
                                                                                       ggaagaTATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCAc     >  2:2895655/1‑93 (MQ=255)
                                                                                       ggaagaTATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCAc     >  2:769413/1‑93 (MQ=255)
                                                                                          ggaTATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGAc  >  1:4040831/1‑93 (MQ=255)
                                                                                          |                                                                                            
TCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGAC  >  NC_000913/373970‑374152

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