Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,637,944 T→C H68H (CAT→CAC uspA → universal stress protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,637,9440TC100.0% 745.1 / NA 233H68H (CAT→CACuspAuniversal stress protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (126/107);  total (126/107)

Only 100 of 254 total aligned reads displayed.

CGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637854‑3638033
                                                                                          |                                                                                         
cGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGc                                                                                         <  1:4147813/93‑1 (MQ=255)
 gTAGATGTAAACTACTCTGACCTATACATCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCa                                                                                        <  1:204194/93‑1 (MQ=255)
 gTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCa                                                                                        <  1:2567557/93‑1 (MQ=255)
   aGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACt                                                                                      <  2:2799257/93‑1 (MQ=255)
   aGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACt                                                                                      <  1:3340927/93‑1 (MQ=255)
      tGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAc                                                                                   >  1:3836208/1‑93 (MQ=255)
       gTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAcc                                                                                  >  2:2224308/1‑93 (MQ=255)
          aaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAg                                                                               >  2:295109/1‑93 (MQ=255)
           actactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTgagcgagc                                                                              >  2:2579000/1‑93 (MQ=255)
           actactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGc                                                                              <  2:3818866/93‑1 (MQ=255)
            ctactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCt                                                                             <  2:3413826/93‑1 (MQ=255)
             tactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCtt                                                                            >  2:159157/1‑93 (MQ=255)
               ctctGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTc                                                                          >  2:1357422/1‑93 (MQ=255)
                tctGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTcc                                                                         >  2:3777004/1‑93 (MQ=255)
                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       >  2:2306162/1‑93 (MQ=255)
                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       >  2:395850/1‑93 (MQ=255)
                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       <  1:3443168/93‑1 (MQ=255)
                    aCCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTa                                                                     >  2:3503284/1‑93 (MQ=255)
                     ccTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTaa                                                                    <  2:1746705/93‑1 (MQ=255)
                     ccTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTaa                                                                    <  1:2460247/93‑1 (MQ=255)
                     ccTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTaa                                                                    >  2:1278828/1‑93 (MQ=255)
                      cTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAAt                                                                   >  1:1012412/1‑93 (MQ=255)
                      cTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAAt                                                                   >  2:1331757/1‑93 (MQ=255)
                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  2:1729659/1‑93 (MQ=255)
                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  2:4376684/1‑93 (MQ=255)
                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  1:1853135/1‑93 (MQ=255)
                        ataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACGTCACGCACTGACCGAGCTTTCCACTAATGc                                                                 >  1:2004632/1‑93 (MQ=255)
                        ataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGc                                                                 >  1:3633903/1‑93 (MQ=255)
                         taCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCa                                                                <  2:3340316/93‑1 (MQ=255)
                         taCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCa                                                                >  1:870152/1‑93 (MQ=255)
                            acCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGc                                                             >  2:4057922/1‑93 (MQ=255)
                             ccGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCt                                                            <  1:1485767/93‑1 (MQ=255)
                               gggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAc                                                          >  1:702376/1‑93 (MQ=255)
                                ggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAcc                                                         >  1:1338320/1‑93 (MQ=255)
                                ggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAcc                                                         <  2:4655908/93‑1 (MQ=255)
                                 gCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAccc                                                        >  2:1014512/1‑93 (MQ=255)
                                   ttATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTATCCaa                                                      >  2:2375631/1‑93 (MQ=255)
                                   ttATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCaa                                                      <  1:2940025/93‑1 (MQ=255)
                                   ttATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCaa                                                      >  1:3112136/1‑93 (MQ=255)
                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     <  2:1802191/93‑1 (MQ=255)
                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     >  2:416598/1‑93 (MQ=255)
                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     >  1:4325079/1‑93 (MQ=255)
                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     <  2:3024043/93‑1 (MQ=255)
                                     aTTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATc                                                    >  1:4687141/1‑93 (MQ=255)
                                     aTTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATc                                                    >  1:441860/1‑93 (MQ=255)
                                     aTTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATc                                                    <  2:3458288/93‑1 (MQ=255)
                                       tGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCAc                                                  <  2:3278323/93‑1 (MQ=255)
                                        gATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACt                                                 <  1:3735743/93‑1 (MQ=255)
                                         aTGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTg                                                >  1:1212459/1‑93 (MQ=255)
                                          tgtgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGa                                               <  1:4457706/93‑1 (MQ=255)
                                             gAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAc                                            >  2:2096108/1‑93 (MQ=255)
                                              aaTCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAcc                                           >  2:769333/1‑93 (MQ=255)
                                               aTCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAccc                                          >  1:327185/1‑93 (MQ=255)
                                               aTCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAccc                                          <  2:4090231/93‑1 (MQ=255)
                                                  tGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGa                                       <  2:4538360/93‑1 (MQ=255)
                                                  tGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGa                                       <  2:1140154/93‑1 (MQ=255)
                                                  tGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGa                                       >  2:2457549/1‑93 (MQ=255)
                                                   gggTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGag                                      <  1:2362568/93‑1 (MQ=255)
                                                    ggTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagc                                     <  2:4452578/93‑1 (MQ=255)
                                                        atatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggca                                 >  1:2821460/1‑93 (MQ=255)
                                                        atatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggca                                 >  1:1443231/1‑93 (MQ=255)
                                                         tatGCAGAAACGCATCTCTGAAGAGATACATCACAAGCTGACCGTGCTTTACACTAACGCAGGCTACTCAATCACTGAAATCCTGtgcagcag                                >  2:290211/1‑93 (MQ=255)
                                                         tatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcag                                >  1:4209782/1‑93 (MQ=255)
                                                         tatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcag                                >  2:1930468/1‑93 (MQ=255)
                                                         tatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcag                                <  2:2423011/93‑1 (MQ=255)
                                                         tatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcag                                >  2:2771623/1‑93 (MQ=255)
                                                                 aaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAcc                        >  1:2609437/1‑93 (MQ=255)
                                                                 aaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGAGAcc                        >  1:1183560/1‑93 (MQ=255)
                                                                   cGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTg                      >  1:4054162/1‑93 (MQ=255)
                                                                   cGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTg                      <  1:1117687/93‑1 (MQ=255)
                                                                     cATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTggg                    >  2:534131/1‑93 (MQ=255)
                                                                       tctcTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGcc                  <  2:263706/93‑1 (MQ=255)
                                                                        ctctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCa                 >  2:4249474/1‑93 (MQ=255)
                                                                        ctctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCa                 <  1:2971754/93‑1 (MQ=255)
                                                                        ctctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCa                 <  2:1392280/93‑1 (MQ=255)
                                                                         tctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAg                >  1:2876971/1‑93 (MQ=255)
                                                                          ctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAgg               >  1:604274/1‑93 (MQ=255)
                                                                           tGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAGCCCTGAGCGGCAGCGGCGACCTGGGCCAGGt              <  2:3608704/93‑1 (MQ=255)
                                                                           tGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAACCCCTGAGCGGCAGCGGCGACCTGGGCCAGGt              <  2:1962857/93‑1 (MQ=255)
                                                                           tGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGt              >  2:2143427/1‑93 (MQ=255)
                                                                            gAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAACCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGtt             >  2:3738271/1‑93 (MQ=255)
                                                                             aaGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTc            <  1:2207253/93‑1 (MQ=255)
                                                                             aaGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTc            >  2:3681014/1‑93 (MQ=255)
                                                                              agagACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCt           >  2:617166/1‑93 (MQ=255)
                                                                               gagaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTg          >  2:3164960/1‑93 (MQ=255)
                                                                               gagaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTg          >  1:818089/1‑93 (MQ=255)
                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        >  1:3517643/1‑93 (MQ=255)
                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        <  2:183465/93‑1 (MQ=255)
                                                                                   cacaTCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCg      >  1:1963291/1‑93 (MQ=255)
                                                                                    acaTCGAGCACTGGCCGAGCTTTACTCTAATGCAGGCTACCCAAACGCTGAAACCCCGAGCTGCAGCGGGGACCTGGGCCAGGTTCTGGTCGa     <  2:1553570/93‑1 (MQ=255)
                                                                                    acaTCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGa     >  1:1879363/1‑93 (MQ=255)
                                                                                    acaTCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGa     <  1:4255174/93‑1 (MQ=255)
                                                                                    acaTCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGa     >  1:2122711/1‑93 (MQ=255)
                                                                                     catcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt    <  2:4625082/93‑1 (MQ=255)
                                                                                      atcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATg   >  1:4119845/1‑93 (MQ=255)
                                                                                      atcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATg   <  1:3911251/93‑1 (MQ=255)
                                                                                      atcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATg   >  1:3056208/1‑93 (MQ=255)
                                                                                       tcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGc  <  2:131857/93‑1 (MQ=255)
                                                                                       tcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGc  <  1:552954/93‑1 (MQ=255)
                                                                                       tcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGc  >  1:1041358/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637854‑3638033

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