Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,419,344 G→A A38V (GCT→GTT)  ydaG ← Rac prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,419,3440GA100.0% 242.7 / NA 74A38V (GCT→GTT) ydaGRac prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/41);  total (33/41)

TTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCG  >  NC_000913/1419214‑1419478
                                                                                                                                  |                                                                                                                                      
ttGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGAATCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTct                                                                                                                                <  2:257056/139‑1 (MQ=255)
   cTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcc                                                                                                                             <  1:641960/139‑1 (MQ=255)
       tGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCtgcctgcc                                                                                                                         <  1:374419/139‑1 (MQ=255)
             aCGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGAt                                                                                                                   <  1:309486/139‑1 (MQ=255)
               gTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATcc                                                                                                                 >  2:1620935/1‑139 (MQ=255)
                  aTGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcc                                                                                                                             <  2:185226/124‑1 (MQ=255)
                  aTGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcc                                                                                                                             >  1:185226/1‑124 (MQ=255)
                   tGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATccc                                                                                                                >  1:1137445/1‑136 (MQ=255)
                   tGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATccc                                                                                                                <  2:1137445/136‑1 (MQ=255)
                   tGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGt                                                                                                             <  2:33939/139‑1 (MQ=255)
                         tGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                       <  2:1206258/139‑1 (MQ=255)
                         tGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                       <  1:893843/139‑1 (MQ=255)
                           aaCGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcc                                                                                                                             >  1:1606773/1‑115 (MQ=255)
                           aaCGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcc                                                                                                                             <  2:1606773/115‑1 (MQ=255)
                           aaCGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgctg                                                                                                     >  2:1217388/1‑139 (MQ=255)
                           aaCGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgctg                                                                                                     >  1:92610/1‑139 (MQ=255)
                           aaCGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgctg                                                                                                     >  1:1153452/1‑139 (MQ=255)
                           aaCGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgctg                                                                                                     >  2:200836/1‑139 (MQ=255)
                              gTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                  <  2:220692/139‑1 (MQ=255)
                                     ccTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATc                                                                                                                  <  1:798410/116‑1 (MQ=255)
                                     ccTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATc                                                                                                                  >  2:798410/1‑116 (MQ=255)
                                     ccTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                       >  1:1661754/1‑127 (MQ=255)
                                     ccTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                       <  2:1661754/127‑1 (MQ=255)
                                        ggTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        <  1:1217388/139‑1 (MQ=255)
                                        ggTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        <  1:1620935/139‑1 (MQ=255)
                                             gggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGcc                                                                                   <  2:92610/139‑1 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATc                                                                                                                  <  1:335799/107‑1 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATc                                                                                                                  >  2:335799/1‑107 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                       >  1:1463671/1‑118 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                       <  2:1463671/118‑1 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        <  1:440377/133‑1 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        >  2:440377/1‑133 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCt                                                                                  >  2:1528416/1‑139 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCt                                                                                  >  1:1071073/1‑139 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCt                                                                                  >  2:280668/1‑139 (MQ=255)
                                                aaCACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTgg                                                                                <  1:280668/139‑1 (MQ=255)
                                                    cgccgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        >  1:935356/1‑127 (MQ=255)
                                                    cgccgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        <  2:935356/127‑1 (MQ=255)
                                                            acCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATa                                                                    >  2:239444/1‑139 (MQ=255)
                                                             ccAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATAc                                                                   <  1:200836/139‑1 (MQ=255)
                                                             ccAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATAc                                                                   <  2:1153452/139‑1 (MQ=255)
                                                                gggATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGc                                                                <  2:1042833/139‑1 (MQ=255)
                                                                   aTTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGcaac                                                             <  1:239444/139‑1 (MQ=255)
                                                                     tCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    >  1:931831/1‑114 (MQ=255)
                                                                     tCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    <  2:931831/114‑1 (MQ=255)
                                                                      cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGcaa                                                              <  1:1305483/135‑1 (MQ=255)
                                                                      cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGcaa                                                              >  2:1305483/1‑135 (MQ=255)
                                                                      cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          >  1:953926/1‑139 (MQ=255)
                                                                           ttgCTGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATc                                                     <  1:41517/139‑1 (MQ=255)
                                                                                gAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATa                                                <  2:1071073/139‑1 (MQ=255)
                                                                                gAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATa                                                <  2:1615791/139‑1 (MQ=255)
                                                                                   gTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTg                                             >  1:430409/1‑139 (MQ=255)
                                                                                        ttttttGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctg                                                                                                                               <  2:457685/52‑1 (MQ=255)
                                                                                        ttttttGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctg                                                                                                                               >  1:457685/1‑52 (MQ=255)
                                                                                        ttttttGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTg                                                                                          >  2:1024155/1‑89 (MQ=255)
                                                                                        ttttttGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTg                                                                                          <  1:1024155/89‑1 (MQ=255)
                                                                                           tttGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACtt                                     >  1:917498/1‑139 (MQ=255)
                                                                                           tttGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACtt                                     >  1:398189/1‑139 (MQ=255)
                                                                                                   ttGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGc                             <  2:953926/139‑1 (MQ=255)
                                                                                                   ttGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGc                             <  1:1528416/139‑1 (MQ=255)
                                                                                                    tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                       >  2:1211174/1‑64 (MQ=255)
                                                                                                    tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                       <  1:1211174/64‑1 (MQ=255)
                                                                                                    tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    <  2:804813/83‑1 (MQ=255)
                                                                                                    tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    >  1:804813/1‑83 (MQ=255)
                                                                                                     gCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCAc                           <  2:430409/139‑1 (MQ=255)
                                                                                                                     gttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTAtctct               >  1:1289504/1‑135 (MQ=255)
                                                                                                                     gttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTAtctct               <  2:1289504/135‑1 (MQ=255)
                                                                                                                        gCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          >  1:768333/1‑89 (MQ=255)
                                                                                                                        gCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          <  2:768333/89‑1 (MQ=255)
                                                                                                                             gTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCa                                                    <  1:1586834/90‑1 (MQ=255)
                                                                                                                             gTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCa                                                    >  2:1586834/1‑90 (MQ=255)
                                                                                                                              tataACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGtt                                                                              <  2:325019/63‑1 (MQ=255)
                                                                                                                              tataACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGtt                                                                              >  1:325019/1‑63 (MQ=255)
                                                                                                                              tataACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCg  >  1:568656/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                      
TTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCG  >  NC_000913/1419214‑1419478

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: ‑