Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 504,622 A→G G25G (GGT→GGC fsr ← putative fosmidomycin efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913504,6220AG100.0% 340.9 / NA 96G25G (GGT→GGCfsrputative fosmidomycin efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (47/49);  total (47/49)

TGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAG  >  NC_000913/504488‑504757
                                                                                                                                      |                                                                                                                                       
tGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCaa                                                                                                                                     <  2:1762157/139‑1 (MQ=255)
   aaGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACg                                                                                                                                  <  1:831426/139‑1 (MQ=255)
       tGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGt                                                                                                                              >  2:167576/1‑139 (MQ=255)
          gggTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGtt                                                                                                                             <  2:1145764/137‑1 (MQ=255)
          gggTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTCAGAGGCTGATAGCACCTAAAATGCCAAACGATGtt                                                                                                                             >  1:1145764/1‑137 (MQ=255)
                 cATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTg                                                                                                                    <  1:167576/139‑1 (MQ=255)
                   tGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCt                                                                                                                      <  2:14496/135‑1 (MQ=255)
                   tGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCt                                                                                                                      >  1:14496/1‑135 (MQ=255)
                        aTCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTg                                                                                                             >  2:1418158/1‑139 (MQ=255)
                         tCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGa                                                                                                            >  2:1644971/1‑139 (MQ=255)
                         tCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGa                                                                                                            >  2:1428390/1‑139 (MQ=255)
                         tCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGa                                                                                                            >  2:574266/1‑139 (MQ=255)
                             cATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCa                                                                                                        >  2:1239099/1‑139 (MQ=255)
                             cATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCa                                                                                                        >  1:779597/1‑139 (MQ=255)
                                  aTGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCcgc                                                                                                   >  2:1079726/1‑139 (MQ=255)
                                       agagAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGcgcccg                                                                                              >  2:954093/1‑139 (MQ=255)
                                       agagAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGcgcccg                                                                                              >  1:451605/1‑139 (MQ=255)
                                       agagAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGcgcccg                                                                                              >  1:65478/1‑139 (MQ=255)
                                        gagaAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGcgcccgc                                                                                             >  2:294692/1‑139 (MQ=255)
                                          gaAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           <  1:1449466/139‑1 (MQ=255)
                                          gaAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           <  1:1418158/139‑1 (MQ=255)
                                          gaAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           <  1:294692/139‑1 (MQ=255)
                                          gaAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGAACATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCCTTGACGCAGCCGcgcccgcaa                                                                                           <  1:1079726/139‑3 (MQ=255)
                                          aaaaaTTCTGACTGAAGCAGCGGATAAATCGCAAGAATCAGCGATTGGATCAGGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           <  1:497127/138‑1 (MQ=255)
                                             aaTTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAg                                                                                        >  2:1796435/1‑139 (MQ=255)
                                                tCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           <  2:166364/133‑1 (MQ=255)
                                                tCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           >  1:166364/1‑133 (MQ=255)
                                                       gaAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAATGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGgttg                                                                              >  2:335052/1‑139 (MQ=255)
                                                         agcagcGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGgttgtt                                                                            <  2:451605/139‑1 (MQ=255)
                                                         agcagcGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGgttgtt                                                                            <  1:148417/139‑1 (MQ=255)
                                                           cagcGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCa                                                                          <  2:65478/139‑1 (MQ=255)
                                                                gATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCAt                                                                     <  1:1007392/139‑1 (MQ=255)
                                                                gATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCAt                                                                     >  1:1224716/1‑139 (MQ=255)
                                                                  tAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTg                                                                   >  1:1617219/1‑139 (MQ=255)
                                                                         ccAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTc                                                                       >  2:181962/1‑128 (MQ=255)
                                                                         ccAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTc                                                                       <  1:181962/128‑1 (MQ=255)
                                                                           aGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGcc                                                                                                                       <  2:1257631/78‑1 (MQ=255)
                                                                           aGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGcc                                                                                                                       >  1:1257631/1‑78 (MQ=255)
                                                                            gAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGaaa                                                         <  1:954093/139‑1 (MQ=255)
                                                                             aaTCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAg                                                        >  2:191545/1‑139 (MQ=255)
                                                                                       tGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGtt                                                                                                                             >  2:1062759/1‑60 (MQ=255)
                                                                                       tGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGtt                                                                                                                             <  1:1062759/60‑1 (MQ=255)
                                                                                         gATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTgggg                                                                                 <  2:228278/102‑1 (MQ=255)
                                                                                         gATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTgggg                                                                                 >  1:228278/1‑102 (MQ=255)
                                                                                           tCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACg                                                                                                                                  >  2:972464/1‑51 (MQ=255)
                                                                                           tCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACg                                                                                                                                  <  1:972464/51‑1 (MQ=255)
                                                                                           tCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                               <  1:638027/134‑1 (MQ=255)
                                                                                           tCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                               >  2:638027/1‑134 (MQ=255)
                                                                                             aTGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGtt                                                                                                                             >  2:58090/1‑54 (MQ=255)
                                                                                             aTGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGtt                                                                                                                             <  1:58090/54‑1 (MQ=255)
                                                                                              tGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCa                                                                      >  1:837571/1‑108 (MQ=255)
                                                                                              tGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCa                                                                      <  2:837571/108‑1 (MQ=255)
                                                                                                tCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCg                                     >  1:209756/1‑139 (MQ=255)
                                                                                                   ttCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGtt                                        >  2:539194/1‑133 (MQ=255)
                                                                                                   ttCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGtt                                        <  1:539194/133‑1 (MQ=255)
                                                                                                    tCAGCAGATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGt                                 <  2:779597/139‑1 (MQ=255)
                                                                                                     caggagATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTa                                <  1:598947/139‑1 (MQ=255)
                                                                                                     cagcagATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           >  1:766030/1‑80 (MQ=255)
                                                                                                     cagcagATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           <  2:766030/80‑1 (MQ=255)
                                                                                                     cagcagATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTa                                <  2:1164100/139‑1 (MQ=255)
                                                                                                     cagcagATGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTa                                <  1:1428390/139‑1 (MQ=255)
                                                                                                           aTGTGAGGGGCTGATAGCACCTAAAATGCCAAACGATTTCCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAAt                          <  1:1796435/139‑1 (MQ=255)
                                                                                                           aTGTGAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAAt                          <  2:209756/139‑1 (MQ=255)
                                                                                                             gtgAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           <  2:1781097/72‑1 (MQ=255)
                                                                                                             gtgAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           >  1:1781097/1‑72 (MQ=255)
                                                                                                             gtgAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                               >  1:528299/1‑116 (MQ=255)
                                                                                                             gtgAGAGGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                               <  2:528299/116‑1 (MQ=255)
                                                                                                                  agGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTgggg                                                                                 <  2:587557/77‑1 (MQ=255)
                                                                                                                  agGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTgggg                                                                                 >  1:587557/1‑77 (MQ=255)
                                                                                                                  agGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGgttgtt                                                                            >  2:1750650/1‑82 (MQ=255)
                                                                                                                  agGCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGgttgtt                                                                            <  1:1750650/82‑1 (MQ=255)
                                                                                                                    gCTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAtt                 >  1:507992/1‑139 (MQ=255)
                                                                                                                     cTGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAttt                >  1:1385200/1‑139 (MQ=255)
                                                                                                                      tGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtg               <  1:574266/139‑1 (MQ=255)
                                                                                                                      tGATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtg               <  1:191545/139‑1 (MQ=255)
                                                                                                                       gATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTc                                                                       <  2:1431025/82‑1 (MQ=255)
                                                                                                                       gATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTc                                                                       >  1:1431025/1‑82 (MQ=255)
                                                                                                                       gATAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtgt              <  2:507992/139‑1 (MQ=255)
                                                                                                                        aTAGCCCCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTc                                                      <  1:518293/98‑1 (MQ=255)
                                                                                                                        aTAGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTc                                                      >  2:518293/1‑98 (MQ=255)
                                                                                                                          aGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCg                                     >  2:1266701/1‑113 (MQ=255)
                                                                                                                          aGCACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCg                                     <  1:1266701/113‑1 (MQ=255)
                                                                                                                            cACCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatta         >  1:1574697/1‑139 (MQ=255)
                                                                                                                             aCCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCa                                                                      >  2:323761/1‑77 (MQ=255)
                                                                                                                             aCCTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCa                                                                      <  1:323761/77‑1 (MQ=255)
                                                                                                                              ccTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           <  1:792738/55‑1 (MQ=255)
                                                                                                                              ccTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCca                                                                                           >  2:792738/1‑55 (MQ=255)
                                                                                                                              ccTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTgggg                                                                                 >  2:1046417/1‑65 (MQ=255)
                                                                                                                              ccTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTgggg                                                                                 <  1:1046417/65‑1 (MQ=255)
                                                                                                                               cTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGa                                    <  2:1076496/109‑1 (MQ=255)
                                                                                                                               cTAAAATGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGa                                    >  1:1076496/1‑109 (MQ=255)
                                                                                                                                  aaaTGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                               >  2:749600/1‑95 (MQ=255)
                                                                                                                                  aaaTGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                               <  2:1036020/95‑1 (MQ=255)
                                                                                                                                  aaaTGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                               <  1:749600/95‑1 (MQ=255)
                                                                                                                                  aaaTGCCAAACGATGTTCGGGCCTTGGTCCTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                               >  1:1036020/1‑95 (MQ=255)
                                                                                                                                   aaTGCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAg  <  2:1385200/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAG  >  NC_000913/504488‑504757

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