Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,478,998 T→G H537P (CAT→CCT)  kefB ← potassium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,478,9980TG100.0% 217.5 / NA 58H537P (CAT→CCT) kefBpotassium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (31/27);  total (31/27)

CACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGATGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAATATGCAAATGCGGAAAGTGCT  >  NC_000913/3478852‑3479148
                                                                                                                                                  |                                                                                                                                                      
cACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGaggca                                                                                                                                                     <  1:630751/150‑1 (MQ=255)
  cGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGaggcatgc                                                                                                                                                  >  1:573568/1‑151 (MQ=255)
        aaaaTCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCg                                                                                                                                            <  1:690658/151‑1 (MQ=255)
         aaaTCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCg                                                                                                                                            >  2:453118/1‑150 (MQ=255)
         aaaTCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGt                                                                                                                                           >  1:718404/1‑151 (MQ=255)
               ctTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGAcc                                                                                                                                       >  1:619677/1‑149 (MQ=255)
               ctTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGAcc                                                                                                                                       >  1:410944/1‑149 (MQ=255)
                ttCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGc                                                                                                                                    >  2:702426/1‑151 (MQ=255)
                 tCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCa                                                                                                                                      <  1:824564/148‑1 (MQ=255)
                 tCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCg                                                                                                                                   >  2:847056/1‑151 (MQ=255)
                       tCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTg                                                                                                                              >  1:273407/1‑150 (MQ=255)
                          cgTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGcc                                                                                                                          >  1:514028/1‑151 (MQ=255)
                           gTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGccc                                                                                                                         >  2:49262/1‑151 (MQ=255)
                           gTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGccc                                                                                                                         >  2:851313/1‑151 (MQ=255)
                             cGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCa                                                                                                                       <  2:787965/151‑1 (MQ=255)
                               ggCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGc                                                                                                                     <  1:669193/151‑1 (MQ=255)
                               ggCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGc                                                                                                                     >  1:171797/1‑151 (MQ=255)
                               ggCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGc                                                                                                                     <  1:456066/151‑1 (MQ=255)
                                     ccTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAAc                                                                                                               >  1:223637/1‑151 (MQ=255)
                                        ggCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACg                                                                                                              <  1:476514/149‑1 (MQ=255)
                                         gcgcGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCAc                                                                                                           >  1:653287/1‑151 (MQ=255)
                                            cgAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTgg                                                                                                        <  2:194099/151‑1 (MQ=255)
                                             gagaAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGga                                                                                                       <  1:614168/151‑1 (MQ=255)
                                                     tGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGtt                                                                                                <  2:373704/150‑1 (MQ=255)
                                                       tACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTc                                                                                              >  2:289914/1‑150 (MQ=255)
                                                         cGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACg                                                                                           >  2:890337/1‑151 (MQ=255)
                                                          ggTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACgg                                                                                          <  2:491794/151‑1 (MQ=255)
                                                           gTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACggg                                                                                         >  1:427121/1‑151 (MQ=255)
                                                             aTCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGaa                                                                                       <  1:634710/151‑1 (MQ=255)
                                                               cGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGaaaa                                                                                     <  1:635225/151‑1 (MQ=255)
                                                               cGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGaaaa                                                                                     <  2:671712/151‑1 (MQ=255)
                                                                 gcatgcatTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACt                                                                                   >  1:44027/1‑151 (MQ=255)
                                                                   atgcatTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGc                                                                                 >  1:872991/1‑151 (MQ=255)
                                                                     gcatTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCg                                                                                <  1:360554/150‑1 (MQ=255)
                                                                      catTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTc                                                                              >  2:903509/1‑151 (MQ=255)
                                                                       atTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCa                                                                             <  1:896433/151‑1 (MQ=255)
                                                                              tGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCcctgc                                                                      >  2:613483/1‑151 (MQ=255)
                                                                                    ctcGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGt                                                                  <  1:435334/149‑1 (MQ=255)
                                                                                          gcaTCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTc                                                          <  1:275694/151‑1 (MQ=255)
                                                                                                 cATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCt                                                   >  1:571674/1‑151 (MQ=255)
                                                                                                  atatCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCtt                                                  >  2:458812/1‑151 (MQ=255)
                                                                                                   tatCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTc                                                 >  2:171776/1‑151 (MQ=255)
                                                                                                     tCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAATTCATGCGCTTcc                                                >  1:43159/1‑150 (MQ=255)
                                                                                                      ccAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCcac                                              <  1:66877/151‑1 (MQ=255)
                                                                                                       cAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCcac                                              >  2:309931/1‑150 (MQ=255)
                                                                                                         ggcggcGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACg                                           <  1:251715/151‑1 (MQ=255)
                                                                                                           cggcgAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTc                                         >  1:525819/1‑151 (MQ=255)
                                                                                                           cggcgAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTc                                         <  1:833600/151‑1 (MQ=255)
                                                                                                             gcgAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCcg                                       <  2:527691/151‑1 (MQ=255)
                                                                                                             gcgAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCcg                                       >  2:616287/1‑151 (MQ=255)
                                                                                                              cgAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTcc                                        >  2:77901/1‑149 (MQ=255)
                                                                                                               gAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCcgcg                                     >  1:483393/1‑151 (MQ=255)
                                                                                                                   aTGCAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCt                                 <  1:110905/151‑1 (MQ=255)
                                                                                                                      cAGTTGCGCTCGCTGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTcg                               >  2:896993/1‑150 (MQ=255)
                                                                                                                                cgctGTGCCTGATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAatat                    <  1:209349/151‑1 (MQ=255)
                                                                                                                                          gATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAATATGCAAATGc            <  2:515744/149‑1 (MQ=255)
                                                                                                                                          gATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAATATGCAAATGc            <  2:791498/149‑1 (MQ=255)
                                                                                                                                          gATGCGGAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAATATGCAAATGCgg          <  1:637592/151‑1 (MQ=255)
                                                                                                                                               ggAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAATATGCAAATGCGGAAAGt     <  2:538456/151‑1 (MQ=255)
                                                                                                                                               ggAGGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAATATGCAAATGCGGAAAGt     >  2:126270/1‑151 (MQ=255)
                                                                                                                                                  tgcatgcCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAATATGCAAATGCGGAAAGTGCt  >  2:529806/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCTCTCGCAGCATCCGCATATCCAGGCGGCGAAAATGCAGTTGCGCTCGCTGTGCCTGATGCGGATGCATGCCAAGCGTGACCAGCGTCTTGCGCCCCAGCTCTAACGCACTGGAGAATGTTTCACGGGAAAACTGCGTCACCCCTGCCTGTAATAACTCATGCGCTTCCACACGTCCGCGCGCTCGCGCAAGAATATGCAAATGCGGAAAGTGCT  >  NC_000913/3478852‑3479148

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