Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 598,492 A→G E338E (GAA→GAG)  cusB → copper/silver efflux system, membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913598,4920AG99.2% 447.6 / NA 125E338E (GAA→GAG) cusBcopper/silver efflux system, membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (65/59);  minor base C (0/1);  total (65/60)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 126 total aligned reads displayed.

GCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTCC  >  NC_000913/598346‑598639
                                                                                                                                                  |                                                                                                                                                   
gcgACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGagcagc                                                                                                                                                 <  2:297616/151‑1 (MQ=255)
  gACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCgg                                                                                                                                               >  4:731276/1‑151 (MQ=255)
    cccGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCgg                                                                                                                                               <  1:1341578/149‑1 (MQ=255)
    cccGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGt                                                                                                                                             <  1:883158/151‑1 (MQ=255)
     ccGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTg                                                                                                                                            >  4:1125595/1‑151 (MQ=255)
     ccGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTg                                                                                                                                            >  1:155638/1‑151 (MQ=255)
       gCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGAt                                                                                                                                          >  4:303202/1‑151 (MQ=255)
        cgcGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGAt                                                                                                                                          <  4:408216/148‑1 (MQ=255)
          cGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTAc                                                                                                                                       >  3:496261/1‑151 (MQ=255)
          cGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTAc                                                                                                                                       <  3:136177/151‑1 (MQ=255)
            cTGCAGCTGCGTCTGGAAGTCGACAACGCCGCCGAGGCCCTAAAACCGGGCATGAACCCCTGGTTGCAACTAAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTAcc                                                                                                                                      <  4:871255/150‑1 (MQ=255)
               cAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGtt                                                                                                                                   <  1:1377697/150‑1 (MQ=255)
                 gctggctcgggAGGTCGACAACGTCAACGAGGCGCTAAAGCCGGCCATGAGCGCCTGGTTTCAATTCAACTCCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATGCCGGCATCGACCAGCGGGTGATTACCGTTGAt                                                                                                                                <  3:315575/142‑1 (MQ=255)
                    gCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGc                                                                                                                              <  1:414525/150‑1 (MQ=255)
                       tCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGAc                                                                                                                          >  3:158616/1‑151 (MQ=255)
                        cTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACg                                                                                                                         >  1:1604623/1‑151 (MQ=255)
                        cTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACg                                                                                                                         >  1:634208/1‑151 (MQ=255)
                                cGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCttt                                                                                                                 >  4:1166806/1‑151 (MQ=255)
                                cGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCttt                                                                                                                 >  2:10392/1‑151 (MQ=255)
                                   cAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGt                                                                                                               >  1:492130/1‑150 (MQ=255)
                                      cGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCg                                                                                                                            >  2:161045/1‑134 (MQ=255)
                                       gCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCg                                                                                                           >  1:1321602/1‑150 (MQ=255)
                                       gCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCg                                                                                                           >  1:1287611/1‑150 (MQ=255)
                                       gCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCg                                                                                                           >  1:1552586/1‑150 (MQ=255)
                                       gCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCg                                                                                                           >  1:1849614/1‑150 (MQ=255)
                                          gacgaGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGaaa                                                                                                        >  4:397829/1‑150 (MQ=255)
                                          gacgaGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAAc                                                                                                       <  3:39862/151‑1 (MQ=255)
                                          gacgaGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAAc                                                                                                       <  4:1543828/151‑1 (MQ=255)
                                          gacgaGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAAc                                                                                                       <  4:339378/151‑1 (MQ=255)
                                            cgaGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAAcgc                                                                                                     <  2:310438/151‑1 (MQ=255)
                                            cgaGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAAcg                                                                                                      <  4:464537/150‑1 (MQ=255)
                                              aGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGt                                                                                                   >  1:802567/1‑151 (MQ=255)
                                              aGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGt                                                                                                   >  4:1480378/1‑151 (MQ=255)
                                               ggCGCTAAACCCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGtt                                                                                                  <  4:1817178/151‑1 (MQ=255)
                                                gcgcTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTg                                                                                                 >  4:1210291/1‑151 (MQ=255)
                                                  gcTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCt                                                                                               >  2:201594/1‑151 (MQ=255)
                                                    tAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGt                                                                                             >  4:32033/1‑151 (MQ=255)
                                                        aCCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTc                                                                                         <  4:1126275/151‑1 (MQ=255)
                                                         ccGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTcc                                                                                        <  1:1166307/151‑1 (MQ=255)
                                                         ccGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTc                                                                                         <  4:352709/150‑1 (MQ=255)
                                                          cGGGCATGAACTCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCa                                                                                       <  3:378320/151‑1 (MQ=255)
                                                          cGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCa                                                                                       <  4:885413/151‑1 (MQ=255)
                                                              cATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGc                                                                                    >  1:252826/1‑150 (MQ=255)
                                                               aTGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCAt                                                                                  >  4:1740111/1‑151 (MQ=255)
                                                                 gAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCg                                                                                <  1:610687/151‑1 (MQ=255)
                                                                    cGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCa                                                                              >  4:283221/1‑150 (MQ=255)
                                                                     gCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCa                                                                              >  4:427021/1‑149 (MQ=255)
                                                                      ccTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAg                                                                            <  4:1885663/150‑1 (MQ=255)
                                                                         ggTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCg                                                                         <  1:852836/150‑1 (MQ=255)
                                                                         ggTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGt                                                                        >  3:194206/1‑151 (MQ=255)
                                                                          gTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTc                                                                       >  4:319155/1‑151 (MQ=255)
                                                                            tGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCAc                                                                     >  1:106513/1‑151 (MQ=255)
                                                                                aCTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGc                                                                  >  4:468676/1‑150 (MQ=255)
                                                                                 cTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCa                                                                 <  2:393695/150‑1 (MQ=255)
                                                                                   cAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCAtt                                                               >  4:1016403/1‑150 (MQ=255)
                                                                                   cAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCAtt                                                               >  1:1281700/1‑150 (MQ=255)
                                                                                   cAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCAt                                                                <  1:893556/149‑1 (MQ=255)
                                                                                   cAACACCGCCAGCGAACCGATGCTGCGCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCAtt                                                               <  4:1410562/150‑1 (MQ=255)
                                                                                    aaCACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTAc                                                             <  2:532121/151‑1 (MQ=255)
                                                                                     acacCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACg                                                            >  1:490240/1‑151 (MQ=255)
                                                                                         cGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGctct                                                        <  2:514911/151‑1 (MQ=255)
                                                                                         cGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGctct                                                        <  2:216455/151‑1 (MQ=255)
                                                                                         cGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGctct                                                        <  1:159708/151‑1 (MQ=255)
                                                                                          gCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCtctg                                                       <  4:487263/151‑1 (MQ=255)
                                                                                           ccAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCtctgg                                                      >  4:1760248/1‑151 (MQ=255)
                                                                                             aGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCtctggtc                                                    >  1:1391728/1‑151 (MQ=255)
                                                                                               cGAACCGTTGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCtctggtctg                                                  >  3:356130/1‑151 (MQ=255)
                                                                                               cGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCtctggtctg                                                  <  4:1392295/151‑1 (MQ=255)
                                                                                                gAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCtctggtctgg                                                 <  2:215324/151‑1 (MQ=255)
                                                                                                gAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCtctggtctgg                                                 >  1:1873443/1‑151 (MQ=255)
                                                                                                 aaCCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGc                                                >  2:288474/1‑151 (MQ=255)
                                                                                                 aaCCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGc                                                <  4:1355508/151‑1 (MQ=255)
                                                                                                  aCCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCg                                               <  1:1116682/151‑1 (MQ=255)
                                                                                                      aTGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGa                                             >  4:1890384/1‑149 (MQ=255)
                                                                                                      aTGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAg                                           <  4:841512/151‑1 (MQ=255)
                                                                                                        gctgctCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAgg                                          <  1:780381/150‑1 (MQ=255)
                                                                                                             tCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGaaa                                     >  2:129978/1‑150 (MQ=255)
                                                                                                              cATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCGTTACGCTCTGGTCTGGCGGAAGGTGaaa                                     <  4:740384/149‑1 (MQ=255)
                                                                                                               aTTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAgg                                  >  3:414298/1‑151 (MQ=255)
                                                                                                               aTTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAgg                                  >  4:109335/1‑151 (MQ=255)
                                                                                                               aTTCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAg                                   <  4:895328/150‑1 (MQ=255)
                                                                                                                ttCCGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAggt                                 >  4:806922/1‑151 (MQ=255)
                                                                                                                  ccGTCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAAGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAggtgg                               <  3:36336/151‑1 (MQ=255)
                                                                                                                     tCACAAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGtt                             <  4:132655/150‑1 (MQ=255)
                                                                                                                       acaAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTcc                          <  3:218732/151‑1 (MQ=255)
                                                                                                                        caAGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCa                         <  1:1780065/151‑1 (MQ=255)
                                                                                                                         aaGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAg                        <  4:1605254/151‑1 (MQ=255)
                                                                                                                         aaGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAg                        <  4:393458/151‑1 (MQ=255)
                                                                                                                         aaGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAg                        <  3:489358/151‑1 (MQ=255)
                                                                                                                          aGCGCTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGCTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGCTTCCAGc                       >  1:1076711/1‑151 (MQ=255)
                                                                                                                            cgcTGATTGATACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCgg                     <  4:1553198/151‑1 (MQ=255)
                                                                                                                               tgattgatACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCt                  >  4:1267549/1‑151 (MQ=255)
                                                                                                                                   tgatACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGtt               <  4:568908/150‑1 (MQ=255)
                                                                                                                                   tgatACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGtt               <  1:730958/150‑1 (MQ=255)
                                                                                                                                   tgatACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTc              >  1:1012536/1‑151 (MQ=255)
                                                                                                                                    gatACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTcc             >  1:1048631/1‑151 (MQ=255)
                                                                                                                                     atACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCt            <  3:207073/151‑1 (MQ=255)
                                                                                                                                      tACCGGCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCtg           >  2:250738/1‑151 (MQ=255)
                                                                                                                                           gCAGCGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCtgattga      <  2:108201/151‑1 (MQ=255)
                                                                                                                                               cGAGCAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTcc  >  4:1356134/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTTGTACCGAAACGCGTTGCTGTTTTCCAGGCATCGCAAGGCGTCACCGCATTACGCTCTGGTCTGGCGGAAGGTGAAAAGGTGGTTTCCAGCGGCCTGTTCCTGATTGATTCC  >  NC_000913/598346‑598639

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