Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 959,388 (T)8→7 intergenic (+70/‑75) aroA → / → ycaL 5‑enolpyruvylshikimate‑3‑phosphate synthetase/putative heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913959,3810T.100.0% 610.0 / NA 136intergenic (+63/‑82)aroA/ycaL5‑enolpyruvylshikimate‑3‑phosphate synthetase/putative heat shock protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (43/93);  total (43/93)

Only 100 of 150 total aligned reads displayed.

GCTGGCGCGGATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTATATATGATTA  >  NC_000913/959291‑959469
                                                                                          |                                                                                        
  tGGCGCGGATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGttttt                                                                                      <  1:1654344/93‑1 (MQ=255)
   ggCGCGGATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGtttttt                                                                                     <  1:1267915/93‑1 (MQ=255)
    gcgcGGATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGttttttt                                                                                    <  2:383916/93‑1 (MQ=255)
    gcgcGGATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGttttttt                                                                                    <  2:2576870/93‑1 (MQ=255)
     cgtggATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGtttttttc                                                                                   <  1:6143697/90‑2 (MQ=255)
        ggATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGtttttttcgga                                                                                >  1:4805820/1‑89 (MQ=255)
        ggATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGtttttttcgga                                                                                <  2:4501092/93‑5 (MQ=255)
        ggATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGtttttttcgga                                                                                <  1:2528768/93‑5 (MQ=255)
         gATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGAc                                                                              >  2:5388711/1‑93 (MQ=255)
         gATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGAc                                                                              >  2:7762488/1‑93 (MQ=255)
           ttAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCt                                                                            <  1:5717053/93‑1 (MQ=255)
             aGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTtg                                                                          <  2:212388/93‑1 (MQ=255)
                cAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGa                                                                       >  2:2733692/1‑93 (MQ=255)
                cAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGa                                                                       >  1:5667523/1‑93 (MQ=255)
                 aGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAg                                                                      <  1:6857063/93‑1 (MQ=255)
                 aGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAg                                                                      <  2:6023024/93‑1 (MQ=255)
                  ggCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGt                                                                     >  2:1229859/1‑93 (MQ=255)
                      gCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCAtt                                                                 <  1:6177471/93‑1 (MQ=255)
                       ccTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCAttt                                                                <  1:583597/93‑1 (MQ=255)
                       ccTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCAttt                                                                >  1:4747909/1‑93 (MQ=255)
                         tgaatgaaCAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTg                                                              >  1:5870857/1‑93 (MQ=255)
                          gaatgaaCAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGa                                                             <  1:1095476/93‑1 (MQ=255)
                           aatgaaCAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGAt                                                            >  1:7294505/1‑93 (MQ=255)
                           aatgaaCAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGAt                                                            <  1:1962898/93‑1 (MQ=255)
                             tgaaCAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTa                                                          >  2:4037967/1‑93 (MQ=255)
                              gaaCAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTaa                                                         <  2:2849102/93‑1 (MQ=255)
                              gaaCAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTaa                                                         <  1:3496536/93‑1 (MQ=255)
                               aacaacGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAAt                                                        >  1:3027212/1‑93 (MQ=255)
                               aacaacGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAAt                                                        >  1:3027211/1‑93 (MQ=255)
                                  aacGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGt                                                     >  2:6500544/1‑93 (MQ=255)
                                   acGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTa                                                    <  2:5903933/93‑1 (MQ=255)
                                    cGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAg                                                   <  1:7804854/93‑1 (MQ=255)
                                    cGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAg                                                   <  1:2650426/93‑1 (MQ=255)
                                     gggCAATAAATGGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATTGTAGc                                                  >  2:4568949/1‑93 (MQ=255)
                                      ggCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCg                                                 <  2:3135463/93‑1 (MQ=255)
                                      ggCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCg                                                 <  2:3485773/93‑1 (MQ=255)
                                         aataaataGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCg                                              >  2:2688510/1‑93 (MQ=255)
                                          ataaataGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGc                                             >  1:505836/1‑93 (MQ=255)
                                          ataaataGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGc                                             >  1:6839950/1‑93 (MQ=255)
                                            aaataGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCtt                                           <  2:7718968/93‑1 (MQ=255)
                                               taGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTc                                        <  1:668250/93‑1 (MQ=255)
                                               taGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTc                                        >  2:5569446/1‑93 (MQ=255)
                                               taGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTc                                        >  2:1508701/1‑93 (MQ=255)
                                                 gCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCaa                                      <  1:4679514/93‑1 (MQ=255)
                                                 gCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCaa                                      <  1:5451074/93‑1 (MQ=255)
                                                 gCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCaa                                      <  1:453385/93‑1 (MQ=255)
                                                 gCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCaa                                      <  1:190252/93‑1 (MQ=255)
                                                  ccAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAAt                                     <  1:1249152/93‑1 (MQ=255)
                                                  ccAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAAt                                     <  2:7290882/93‑1 (MQ=255)
                                                   cAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATg                                    <  2:561120/93‑1 (MQ=255)
                                                   cAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATg                                    >  2:7011796/1‑93 (MQ=255)
                                                   cAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATg                                    <  1:1332058/93‑1 (MQ=255)
                                                   cAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATg                                    <  1:1332059/93‑1 (MQ=255)
                                                    aaaTCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGt                                   <  1:188289/93‑1 (MQ=255)
                                                    aaaTCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGt                                   <  2:2994598/93‑1 (MQ=255)
                                                     aaTCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTa                                  >  1:4805011/1‑93 (MQ=255)
                                                     aaTCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTa                                  <  1:7495365/93‑1 (MQ=255)
                                                      aTCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTaa                                 <  1:4519981/93‑1 (MQ=255)
                                                      aTCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTaa                                 <  1:2099523/93‑1 (MQ=255)
                                                      aTCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTaa                                 >  1:14417/1‑93 (MQ=255)
                                                       tctttcttTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAg                                <  1:4172174/93‑1 (MQ=255)
                                                       tctttcttTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAg                                <  1:5059417/93‑1 (MQ=255)
                                                       tctttcttTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAg                                >  2:3214951/1‑93 (MQ=255)
                                                       tctttcttTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAg                                <  1:964911/93‑1 (MQ=255)
                                                          ttctttATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAgttg                             <  2:1028419/93‑1 (MQ=255)
                                                            ctttATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAgttgtt                           <  2:5829721/93‑1 (MQ=255)
                                                            ctttATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAgttgtt                           >  2:7436149/1‑93 (MQ=255)
                                                            ctttATCAAAACGTCGGCACATTGCCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAgttgtt                           <  2:711523/93‑1 (MQ=255)
                                                             tttATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGttgttg                          <  2:1626964/93‑1 (MQ=255)
                                                              ttATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGa                         <  2:5387191/93‑1 (MQ=255)
                                                              ttATCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGa                         <  2:4451246/93‑1 (MQ=255)
                                                                aTCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATa                       <  2:7337044/93‑1 (MQ=255)
                                                                aTCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATa                       <  1:3949246/93‑1 (MQ=255)
                                                                aTCAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATa                       <  1:2679305/93‑1 (MQ=255)
                                                                  cAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACa                     <  1:823082/93‑1 (MQ=255)
                                                                  cAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACa                     <  1:2310283/93‑1 (MQ=255)
                                                                  cAAAACGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACa                     <  1:6977790/93‑1 (MQ=255)
                                                                  cAAAACGGCGGCACATTGTAGGCTGTTTTTGCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCCCTTGTCAATGTAAGTTGTTGATACa                     <  1:823084/93‑1 (MQ=255)
                                                                   aaaaCGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGCTGTTGATACat                    >  2:7506351/1‑93 (MQ=255)
                                                                    aaaCGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACata                   >  2:7107387/1‑93 (MQ=255)
                                                                       cGTCGGCTCATTGTCGGCGTGTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACataata                <  2:6413380/93‑1 (MQ=255)
                                                                       cGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACataata                >  2:2568067/1‑93 (MQ=255)
                                                                       cGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACataata                <  1:1738978/93‑1 (MQ=255)
                                                                       cGTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACataata                <  2:6413381/93‑1 (MQ=255)
                                                                        gTCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAatat               >  2:4908056/1‑93 (MQ=255)
                                                                         tCGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATAtt              <  2:2008406/93‑1 (MQ=255)
                                                                          cGGCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATAttt             <  1:3724156/93‑1 (MQ=255)
                                                                           ggCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTta            <  1:4057539/93‑1 (MQ=255)
                                                                           ggCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTta            <  2:5229588/93‑1 (MQ=255)
                                                                           ggCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTta            <  1:4157607/93‑1 (MQ=255)
                                                                           ggCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTta            <  1:1448091/93‑1 (MQ=255)
                                                                           ggCACATTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGATTGTCAATGTAAGTTGTTGATACATAATATTta            >  1:1095751/1‑93 (MQ=255)
                                                                             cacaTTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTtata          >  2:1662015/1‑93 (MQ=255)
                                                                              acaTTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTtatat         >  1:7203028/1‑93 (MQ=255)
                                                                                aTTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTatatat       >  2:2285085/1‑93 (MQ=255)
                                                                                aTTGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTatatat       <  1:5456964/93‑1 (MQ=255)
                                                                                 ttGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTATATATg      <  2:2429263/93‑1 (MQ=255)
                                                                                  tGTCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTATATATGa     <  1:7762489/93‑1 (MQ=255)
                                                                                    tCGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTATATATGAtt   >  1:652870/1‑93 (MQ=255)
                                                                                     cGGCGTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTATATATGATTa  <  2:3387345/93‑1 (MQ=255)
                                                                                          |                                                                                        
GCTGGCGCGGATTAGCCAGGCAGCCTGAATGAACAACGGGCAATAAATAGCCAAATCTTTCTTTATCAAAACGTCGGCACATTGTCGGCGTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTATATATGATTA  >  NC_000913/959291‑959469

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