Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 229,683 A→G I173V (ATC→GTC)  yafB → putative aldose reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913229,6830AG99.3% 960.4 / ‑3.9 304I173V (ATC→GTC) yafBputative aldose reductase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (145/157);  total (147/157)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 310 total aligned reads displayed.

GTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACT  >  NC_000913/229593‑229772
                                                                                          |                                                                                         
 ttGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTcc                                                                                        >  2:491657/1‑93 (MQ=255)
  tGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCa                                                                                       >  2:1663063/1‑93 (MQ=255)
  tGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCa                                                                                       <  1:5753333/93‑1 (MQ=255)
   ggTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCat                                                                                      >  2:4261906/1‑93 (MQ=255)
       cTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTa                                                                                  >  2:6567630/1‑93 (MQ=255)
       cTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTa                                                                                  <  2:7434863/93‑1 (MQ=255)
        tGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTAc                                                                                 <  2:7641280/93‑1 (MQ=255)
         gAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACt                                                                                >  1:5969473/1‑93 (MQ=255)
         gAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACt                                                                                <  1:395417/93‑1 (MQ=255)
            aaCATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCTTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTcc                                                                             <  1:1347977/93‑1 (MQ=255)
            aaCATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTcc                                                                             >  1:5835511/1‑93 (MQ=255)
            aaCATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTcc                                                                             <  1:3229987/93‑1 (MQ=255)
             aCATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCt                                                                            >  1:2521871/1‑93 (MQ=255)
               aTCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCtat                                                                          >  1:7105759/1‑93 (MQ=255)
               aTCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCtat                                                                          >  1:4312903/1‑93 (MQ=255)
                tCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCtata                                                                         >  2:3317986/1‑93 (MQ=255)
                tCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCtata                                                                         <  2:2593290/93‑1 (MQ=255)
                tCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCtata                                                                         >  2:1492268/1‑93 (MQ=255)
                tCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCtata                                                                         >  1:2995938/1‑93 (MQ=255)
                 cGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTatat                                                                        <  2:4055419/93‑1 (MQ=255)
                  gCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGGCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATg                                                                       >  1:5735040/1‑93 (MQ=255)
                  gCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATg                                                                       <  1:5324243/93‑1 (MQ=255)
                    tactaACCAGTTTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAAAAGACCGGCGTCCATATTACTTCCTATATGAc                                                                     >  2:6383582/1‑93 (MQ=255)
                    tactaACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGAc                                                                     >  2:6237051/1‑93 (MQ=255)
                    tactaACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGAc                                                                     >  2:2200975/1‑93 (MQ=255)
                     actaACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACg                                                                    >  1:1794721/1‑93 (MQ=255)
                     actaACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACg                                                                    >  2:6065181/1‑93 (MQ=255)
                     actaACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACg                                                                    >  2:5447514/1‑93 (MQ=255)
                      ctaACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGc                                                                   >  1:723418/1‑93 (MQ=255)
                       taACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCt                                                                  >  1:4101818/1‑93 (MQ=255)
                        aaCCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTg                                                                 <  2:1116568/93‑1 (MQ=255)
                         aCCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTgg                                                                <  2:1046661/93‑1 (MQ=255)
                         aCCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTgg                                                                <  2:1857126/93‑1 (MQ=255)
                          ccAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGc                                                               >  2:2418483/1‑93 (MQ=255)
                          ccAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGc                                                               <  1:123733/93‑1 (MQ=255)
                           cAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCg                                                              <  1:6269323/93‑1 (MQ=255)
                             gATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTa                                                            >  2:4167824/1‑93 (MQ=255)
                               ttGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATg                                                          <  2:6455497/93‑1 (MQ=255)
                               ttGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATg                                                          <  2:2916538/93‑1 (MQ=255)
                                   aCTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTaa                                                      >  2:225288/1‑93 (MQ=255)
                                   aCTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTaa                                                      <  1:6408854/93‑1 (MQ=255)
                                   aCTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTaa                                                      <  2:5582930/93‑1 (MQ=255)
                                    ctctctCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAg                                                     <  2:3046598/93‑1 (MQ=255)
                                    ctctctCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAg                                                     <  2:4465007/93‑1 (MQ=255)
                                     tctctcCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAgg                                                    <  2:1207583/93‑1 (MQ=255)
                                        ctcCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGccc                                                 <  1:7642004/93‑1 (MQ=255)
                                         tcCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                >  1:778633/1‑93 (MQ=255)
                                         tcCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                <  2:7879492/93‑1 (MQ=255)
                                         tcCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                >  1:5888311/1‑93 (MQ=255)
                                         tcCTTAACTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                <  1:1951951/93‑1 (MQ=255)
                                            ttATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaa                                             >  2:7124564/1‑93 (MQ=255)
                                            ttATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaa                                             >  1:6782486/1‑93 (MQ=255)
                                             tATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                            >  1:352906/1‑93 (MQ=255)
                                               tCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGa                                          <  2:7797637/93‑1 (MQ=255)
                                               tCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGa                                          <  1:4047570/93‑1 (MQ=255)
                                                cTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGAt                                         >  2:5047911/1‑93 (MQ=255)
                                                 tGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATg                                        >  2:3627780/1‑93 (MQ=255)
                                                  gCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGa                                       >  1:5448057/1‑93 (MQ=255)
                                                    aaaaCCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAgg                                     >  2:4558662/1‑93 (MQ=255)
                                                     aaaCCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGt                                    <  2:6026356/93‑1 (MQ=255)
                                                     aaaCCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGt                                    >  1:1594725/1‑93 (MQ=255)
                                                      aaCCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGtt                                   >  2:3230670/1‑93 (MQ=255)
                                                       aCCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTa                                  <  1:6876453/93‑1 (MQ=255)
                                                         cGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAtt                                >  1:212478/1‑93 (MQ=255)
                                                         cGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAtt                                <  1:4141808/93‑1 (MQ=255)
                                                            aaaGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCt                             >  1:1285942/1‑93 (MQ=255)
                                                            aaaGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCt                             >  2:7357210/1‑93 (MQ=255)
                                                             aaGTGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTc                            >  1:929099/1‑93 (MQ=255)
                                                                tGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTa                         <  1:4489390/93‑1 (MQ=255)
                                                                tGGTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTa                         <  2:6945524/93‑1 (MQ=255)
                                                                 ggTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTAt                        <  2:3808158/93‑1 (MQ=255)
                                                                  gTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATc                       <  1:119508/93‑1 (MQ=255)
                                                                  gTTGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATc                       >  2:7716808/1‑93 (MQ=255)
                                                                   ttGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCg                      <  1:2478987/93‑1 (MQ=255)
                                                                   ttGCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCg                      >  1:5050899/1‑93 (MQ=255)
                                                                     gCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                    >  2:2968085/1‑93 (MQ=255)
                                                                     gCCTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                    <  1:3049510/93‑1 (MQ=255)
                                                                       cTGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGc                  >  1:1165465/1‑93 (MQ=255)
                                                                        tGGGCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCt                 >  1:6493127/1‑93 (MQ=255)
                                                                          ggCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaa               <  1:2393530/93‑1 (MQ=255)
                                                                           gCTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa              >  2:5637357/1‑93 (MQ=255)
                                                                            cTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAac             <  2:5316587/93‑1 (MQ=255)
                                                                            cTAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAac             >  1:6964317/1‑93 (MQ=255)
                                                                             tAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAaca            >  2:7062335/1‑93 (MQ=255)
                                                                             tAAACAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAaca            <  2:2939804/93‑1 (MQ=255)
                                                                              aaaCAGCACGGCGTCCATATTACTTCCTATCTGACGCTGGCGTATGGTTAGGCCCTGACAGATGAGGTTACTGCTCGTATCGCAGCTAAacac           >  2:5257441/1‑93 (MQ=255)
                                                                              aaaCAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAacac           >  1:3298571/1‑93 (MQ=255)
                                                                               aaCAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAAcaca          >  2:3493932/1‑93 (MQ=255)
                                                                                aCAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTTGTATCGCAGCTAAACACaa         <  1:4599396/93‑1 (MQ=255)
                                                                                aCAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACaa         <  2:877497/93‑1 (MQ=255)
                                                                                aCAGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACaa         >  2:5945792/1‑93 (MQ=255)
                                                                                  aGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATg       <  1:558826/93‑1 (MQ=255)
                                                                                  aGCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATg       <  1:2713707/93‑1 (MQ=255)
                                                                                   gCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc      <  1:4651312/93‑1 (MQ=255)
                                                                                   gCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc      <  1:3593062/93‑1 (MQ=255)
                                                                                   gCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc      >  1:2798347/1‑93 (MQ=255)
                                                                                   gCACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc      >  1:1679041/1‑93 (MQ=255)
                                                                                    cACGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCg     <  1:2322627/93‑1 (MQ=255)
                                                                                      cGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGAc   >  1:828928/1‑93 (MQ=255)
                                                                                      cGGCGTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGAc   <  1:4810346/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACT  >  NC_000913/229593‑229772

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