Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 601,318 A→C N869H (AAC→CAC)  cusA → copper/silver efflux system, membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913601,3180AC100.0% 223.5 / NA 59N869H (AAC→CAC) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (34/25);  total (34/25)

CCCGACGAGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGG  >  NC_000913/601170‑601464
                                                                                                                                                    |                                                                                                                                                  
cccGACGAGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCAc                                                                                                                                                  <  2:496029/151‑1 (MQ=255)
       aGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGc                                                                                                                                           <  1:317956/151‑1 (MQ=255)
        gCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCt                                                                                                                                          <  2:189548/151‑1 (MQ=255)
                   atCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAGCCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATg                                                                                                                                <  2:783289/150‑1 (MQ=255)
                         gCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCg                                                                                                                          >  1:832164/1‑150 (MQ=255)
                             gcgATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGAc                                                                                                                     >  2:432494/1‑151 (MQ=255)
                             gcgATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGAc                                                                                                                     <  1:652716/151‑1 (MQ=255)
                               gATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGt                                                                                                                   <  1:79779/151‑1 (MQ=255)
                                  cGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTtg                                                                                                                 >  2:259411/1‑150 (MQ=255)
                                   gTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTtgat                                                                                                               >  1:815080/1‑151 (MQ=255)
                                             gtCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTc                                                                                                     <  2:449631/151‑1 (MQ=255)
                                              tCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCg                                                                                                    >  2:128054/1‑151 (MQ=255)
                                                 gtggtTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGc                                                                                                 >  2:633695/1‑151 (MQ=255)
                                                    gtTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCtgt                                                                                              >  1:5703/1‑151 (MQ=255)
                                                        aCGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTa                                                                                          <  2:735728/151‑1 (MQ=255)
                                                         cGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTAt                                                                                         >  1:123935/1‑151 (MQ=255)
                                                           aTTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCt                                                                                       >  2:729082/1‑151 (MQ=255)
                                                             ttGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTgg                                                                                     >  1:779553/1‑151 (MQ=255)
                                                                cAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGt                                                                                  >  2:515012/1‑151 (MQ=255)
                                                                 aaaaaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGtt                                                                                 >  2:764068/1‑151 (MQ=255)
                                                                   aaaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTcc                                                                               >  1:427410/1‑151 (MQ=255)
                                                                    aaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCg                                                                              <  2:799225/151‑1 (MQ=255)
                                                                      gCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTc                                                                            >  1:137846/1‑151 (MQ=255)
                                                                         aTAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCggg                                                                         >  1:448605/1‑151 (MQ=255)
                                                                           aGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGt                                                                        <  1:858091/150‑1 (MQ=255)
                                                                              tGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGc                                                                    >  1:811403/1‑151 (MQ=255)
                                                                                aaaaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTgg                                                                      >  2:477265/1‑147 (MQ=255)
                                                                                 aaaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGaa                                                                 <  1:255383/151‑1 (MQ=255)
                                                                                 aaaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGaa                                                                 <  2:767367/151‑1 (MQ=255)
                                                                                 aaaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGa                                                                  >  1:775830/1‑150 (MQ=255)
                                                                                 aaaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGa                                                                  >  2:352007/1‑150 (MQ=255)
                                                                                  aaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGa                                                                  >  1:844820/1‑149 (MQ=255)
                                                                                  aaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGa                                                                  >  1:885475/1‑149 (MQ=255)
                                                                                   aaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGc                                                               >  1:130806/1‑151 (MQ=255)
                                                                                       gCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGtt                                                            <  1:877946/150‑1 (MQ=255)
                                                                                                 ccGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATc                                                  >  1:772176/1‑150 (MQ=255)
                                                                                                 ccGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATc                                                  >  1:530391/1‑150 (MQ=255)
                                                                                                 ccGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATc                                                  >  2:421194/1‑150 (MQ=255)
                                                                                                  cGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATcag                                                <  1:522136/151‑1 (MQ=255)
                                                                                                     gCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATcagcag                                             <  1:223619/151‑1 (MQ=255)
                                                                                                     gCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATcagcag                                             <  1:551920/151‑1 (MQ=255)
                                                                                                      cACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCagcagc                                            <  2:435656/151‑1 (MQ=255)
                                                                                                       aCCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCg                                           <  1:96321/151‑1 (MQ=255)
                                                                                                       aCCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCg                                           >  1:618659/1‑151 (MQ=255)
                                                                                                           gCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTAcc                                       >  2:565886/1‑151 (MQ=255)
                                                                                                            cGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCg                                      >  2:282272/1‑151 (MQ=255)
                                                                                                             gTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGt                                     >  2:852820/1‑151 (MQ=255)
                                                                                                                gCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTg                                  <  1:767798/151‑1 (MQ=255)
                                                                                                                        cGGGCAGTTCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCtggtg                          <  2:399141/151‑1 (MQ=255)
                                                                                                                               ttCGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGggcggc                    >  2:632650/1‑150 (MQ=255)
                                                                                                                                 cGAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATc                 <  1:622166/151‑1 (MQ=255)
                                                                                                                                  gAGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCt                <  2:386051/151‑1 (MQ=255)
                                                                                                                                   aGCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTg               >  1:490155/1‑151 (MQ=255)
                                                                                                                                    gCTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTgg              <  1:487639/151‑1 (MQ=255)
                                                                                                                                     cTACTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGt             >  1:647393/1‑151 (MQ=255)
                                                                                                                                       aCTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGtt            <  1:297381/150‑1 (MQ=255)
                                                                                                                                       aCTGGAGCGCGCCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGCCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGt             >  2:844701/1‑149 (MQ=255)
                                                                                                                                             gcgcgcCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTgtgg     <  1:432251/151‑1 (MQ=255)
                                                                                                                                                 gcCCACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGtggtgg  >  2:179079/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCGACGAGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGG  >  NC_000913/601170‑601464

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