Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,281,472 A→G E498G (GAG→GGG)  yjcE → putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,281,4720AG100.0% 260.1 / NA 68E498G (GAG→GGG) yjcEputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/32);  total (36/32)

GGAGGAGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGC  >  NC_000913/4281338‑4281604
                                                                                                                                      |                                                                                                                                    
ggaggaGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAAc                                                                                                                                  >  1:713726/1‑139 (MQ=255)
  aggagCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCt                                                                                                                                >  1:676165/1‑139 (MQ=255)
    gagCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTg                                                                                                                              <  1:457706/139‑1 (MQ=255)
    gagCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTg                                                                                                                              >  2:602163/1‑139 (MQ=255)
       cGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGc                                                                                                                           >  2:309653/1‑139 (MQ=255)
                   gATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTc                                                                                                               >  2:621863/1‑139 (MQ=255)
                    aTACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCt                                                                                                              >  1:872576/1‑139 (MQ=255)
                             agaATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                                     <  2:872576/139‑1 (MQ=255)
                             agaATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                                     <  1:621863/139‑1 (MQ=255)
                              gaATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTg                                                                                                    >  2:1021092/1‑139 (MQ=255)
                              gaATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTg                                                                                                    <  1:733620/139‑1 (MQ=255)
                                    cGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCt                                                                                              >  1:991726/1‑139 (MQ=255)
                                         aCCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACg                                                                                         >  2:1033834/1‑139 (MQ=255)
                                         aCCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACg                                                                                         >  1:180810/1‑139 (MQ=255)
                                           cAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTg                                                                                       <  2:676165/139‑1 (MQ=255)
                                               tgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                   <  1:602163/139‑1 (MQ=255)
                                                 ctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCg                                                                                                                          <  1:813028/98‑1 (MQ=255)
                                                 ctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCg                                                                                                                          >  2:813028/1‑98 (MQ=255)
                                                 ctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCg                                                                                                                 <  2:1022329/107‑1 (MQ=255)
                                                 ctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCg                                                                                                                 >  1:1022329/1‑107 (MQ=255)
                                                   tACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCg                                                                                                                          <  1:728864/96‑1 (MQ=255)
                                                   tACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCg                                                                                                                          >  2:728864/1‑96 (MQ=255)
                                                          gTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTaccacc                                                                        >  2:453911/1‑139 (MQ=255)
                                                                 ctcGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                                 <  1:612132/139‑1 (MQ=255)
                                                                 ctcGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                                 <  2:991726/139‑1 (MQ=255)
                                                                 ctcGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                                 <  1:1021092/139‑1 (MQ=255)
                                                                 ctcGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                                 <  1:450603/139‑1 (MQ=255)
                                                                   cgcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                               <  2:76005/139‑1 (MQ=255)
                                                                    gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTac                                                                            <  2:32932/125‑1 (MQ=255)
                                                                    gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTac                                                                            >  1:32932/1‑125 (MQ=255)
                                                                       tCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                                    <  1:564308/114‑1 (MQ=255)
                                                                       tCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                                    >  2:564308/1‑114 (MQ=255)
                                                                        cATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACgg                                                          <  2:713726/139‑1 (MQ=255)
                                                                         aTTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                   <  2:631484/113‑1 (MQ=255)
                                                                         aTTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                   >  1:631484/1‑113 (MQ=255)
                                                                         aTTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACggg                                                         >  1:48990/1‑139 (MQ=255)
                                                                         aTTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACggg                                                         >  2:1117434/1‑139 (MQ=255)
                                                                              tAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATc                                                    <  2:516350/139‑1 (MQ=255)
                                                                                 ccTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGc                                                 >  1:859653/1‑139 (MQ=255)
                                                                                    gCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAg                                                                                                                            <  1:1099343/61‑1 (MQ=255)
                                                                                    gCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAg                                                                                                                            >  2:1099343/1‑61 (MQ=255)
                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                  <  2:699796/118‑1 (MQ=255)
                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                             >  1:918846/1‑139 (MQ=255)
                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                             >  2:988490/1‑139 (MQ=255)
                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTGCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                  >  1:699796/1‑118 (MQ=255)
                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCATGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                  >  1:549249/1‑118 (MQ=255)
                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCATGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                  <  2:549249/118‑1 (MQ=255)
                                                                                       tcgtcgCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGaa                                           >  2:1114687/1‑139 (MQ=255)
                                                                                        cgtcgCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGaaa                                          <  2:289269/139‑1 (MQ=255)
                                                                                                       cGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcg                                                                    <  2:400672/98‑1 (MQ=255)
                                                                                                       cGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcg                                                                    >  1:400672/1‑98 (MQ=255)
                                                                                                           acgacgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCa                       <  1:453911/139‑1 (MQ=255)
                                                                                                            cgacgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCaaa                                 <  2:966029/128‑1 (MQ=255)
                                                                                                            cgacgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCaaa                                 >  1:966029/1‑128 (MQ=255)
                                                                                                            cgacgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCAc                      >  1:390463/1‑139 (MQ=255)
                                                                                                               cgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcg                                                                    >  1:998538/1‑90 (MQ=255)
                                                                                                               cgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcg                                                                    <  2:998538/90‑1 (MQ=255)
                                                                                                               cgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGAt                   >  2:286039/1‑139 (MQ=255)
                                                                                                                gTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGAtc                  >  2:471990/1‑139 (MQ=255)
                                                                                                                 ttGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGAtct                 >  1:230458/1‑139 (MQ=255)
                                                                                                                  tGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTACAAAAATTACTGCACGAtctc                >  2:704371/1‑139 (MQ=255)
                                                                                                                     aaGTTACGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGAt             <  1:304529/139‑1 (MQ=255)
                                                                                                                            gTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCt                                      >  1:996381/1‑107 (MQ=255)
                                                                                                                            gTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCt                                      <  2:996381/107‑1 (MQ=255)
                                                                                                                              gCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGaa    <  2:48990/139‑1 (MQ=255)
                                                                                                                               cAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCa                       <  2:1105137/119‑1 (MQ=255)
                                                                                                                               cAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCa                       >  1:1105137/1‑119 (MQ=255)
                                                                                                                                aGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGc  <  1:286039/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GGAGGAGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGC  >  NC_000913/4281338‑4281604

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