Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,252,701 C→T P66S (CCG→TCG)  ubiC → chorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,252,7010CT100.0% 220.4 / NA 63P66S (CCG→TCG) ubiCchorismate lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/33);  total (30/33)

GATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAA  >  NC_000913/4252569‑4252829
                                                                                                                                    |                                                                                                                                
gATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCgctgct                                                                                                                           >  1:295156/1‑140 (MQ=255)
gATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCgctgct                                                                                                                           >  1:86730/1‑140 (MQ=255)
gATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCgctgct                                                                                                                           <  2:497296/140‑1 (MQ=255)
gATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCgctgct                                                                                                                           >  2:648766/1‑140 (MQ=255)
     gCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGa                                                                                                                      <  1:663896/140‑1 (MQ=255)
      cAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGaa                                                                                                                     <  2:853645/140‑1 (MQ=255)
      cAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGaa                                                                                                                     >  1:499452/1‑140 (MQ=255)
       aaCTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAACCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAGTGAAATCCCCGAAGAACTGTCGCTGCTGCCGaaa                                                                                                                    <  2:765077/140‑1 (MQ=255)
          tGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAagag                                                                                                                 >  1:731215/1‑140 (MQ=255)
           gCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGt                                                                                                                <  1:433913/140‑1 (MQ=255)
             tCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGtct                                                                                                              <  1:413485/140‑1 (MQ=255)
                 cTGGCTGTTGCTGGAGGGTTCCATGACAAAACGTTTTCAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGGTTGGCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGtt                                                                                                          >  2:747237/1‑140 (MQ=255)
                 cTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGtct                                                                                                              >  1:308381/1‑136 (MQ=255)
                 cTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGtct                                                                                                              <  2:308381/136‑1 (MQ=255)
                 cTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGtt                                                                                                          <  1:827631/140‑1 (MQ=255)
                 cTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGtt                                                                                                          >  1:172232/1‑140 (MQ=255)
                           cTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACg                                                                                                <  2:499452/140‑1 (MQ=255)
                              gaggaTTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGa                                                                                             >  2:458968/1‑140 (MQ=255)
                              gaggaTTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGa                                                                                             >  2:841334/1‑140 (MQ=255)
                               aggaTTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGaa                                                                                            <  2:444683/140‑1 (MQ=255)
                                 gaTTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAAt                                                                                          <  1:841334/140‑1 (MQ=255)
                                  aTTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAAtt                                                                                         >  1:34576/1‑140 (MQ=255)
                                   ttCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAAttt                                                                                        >  1:379965/1‑140 (MQ=255)
                                     ccATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTg                                                                                      >  2:434331/1‑140 (MQ=255)
                                           cAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTAtgt                                                                                >  2:799318/1‑140 (MQ=255)
                                             aaaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              <  2:731215/140‑1 (MQ=255)
                                             aaaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              <  2:379965/140‑1 (MQ=255)
                                             aaaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              <  2:259105/140‑1 (MQ=255)
                                             aaaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              <  2:565161/140‑1 (MQ=255)
                                              aaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGcc                                                                             >  1:656413/1‑140 (MQ=255)
                                              aaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGcc                                                                             >  1:74252/1‑140 (MQ=255)
                                              aaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGcc                                                                             >  2:684406/1‑140 (MQ=255)
                                              aaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              <  2:796951/139‑1 (MQ=255)
                                              aaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              >  1:358086/1‑139 (MQ=255)
                                              aaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              <  2:358086/139‑1 (MQ=255)
                                              aaCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              >  1:796951/1‑139 (MQ=255)
                                               aCGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCg                                                                            <  1:434331/140‑1 (MQ=255)
                                                  ttttGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATg                                                                         <  1:684406/140‑1 (MQ=255)
                                                    ttGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGt                                                                       <  1:648766/140‑1 (MQ=255)
                                                    ttGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGt                                                                       <  1:579219/140‑1 (MQ=255)
                                                     tGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTg                                                                      >  2:770942/1‑140 (MQ=255)
                                                               ggAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGACGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              <  2:46884/122‑1 (MQ=255)
                                                               ggAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGGTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGc                                                                              >  1:46884/1‑122 (MQ=255)
                                                                gAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCtt                                                           <  1:770942/140‑1 (MQ=255)
                                                                 aaaaaCGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTg                                                          <  2:797780/140‑1 (MQ=255)
                                                                     aCGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCgg                                                      <  2:161144/140‑1 (MQ=255)
                                                                        gTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCg                                                   <  2:295156/140‑1 (MQ=255)
                                                                           aGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTAc                                                <  1:799318/140‑1 (MQ=255)
                                                                             cGTGATGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACcg                                              <  1:56419/140‑1 (MQ=255)
                                                                                 aCGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGACATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTa                                                 <  2:770270/133‑1 (MQ=255)
                                                                                 aCGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTa                                                 >  1:770270/1‑133 (MQ=255)
                                                                                   gatgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACcgtcg                                           >  2:96823/1‑137 (MQ=255)
                                                                                   gatgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACcgtcg                                           <  1:96823/137‑1 (MQ=255)
                                                                                   gatgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTc                                        <  2:74252/140‑1 (MQ=255)
                                                                                       atCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCtgt                                    >  1:647882/1‑140 (MQ=255)
                                                                                           gcgAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGtt                                         >  2:71657/1‑131 (MQ=255)
                                                                                           gcgAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGtt                                         <  1:71657/131‑1 (MQ=255)
                                                                                            cgAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCaa                               >  1:787405/1‑140 (MQ=255)
                                                                                                 ggTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTa                          >  1:61206/1‑140 (MQ=255)
                                                                                                   tttGTCGAGCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGGGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAg                        <  1:334567/140‑1 (MQ=255)
                                                                                                           gCAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCgg                >  2:495332/1‑140 (MQ=255)
                                                                                                            cAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGa               >  2:536020/1‑140 (MQ=255)
                                                                                                            cAGAATGAAATCCCCGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGa               >  1:818755/1‑140 (MQ=255)
                                                                                                                         ccGAAGAACTGTCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACaa  >  2:569805/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                
GATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAA  >  NC_000913/4252569‑4252829

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: