Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,710 C→T 7.6% intergenic (+27/+18) fadH → / ← higA 2,4‑dienoyl‑CoA reductase/antitoxin/DNA‑binding transcriptional repressor HigA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,233,7100CT 7.6% 193.8 / 4.8 66intergenic (+27/+18)fadH/higA2,4‑dienoyl‑CoA reductase/antitoxin/DNA‑binding transcriptional repressor HigA
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base T (3/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

TGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATC  >  NC_000913/3233562‑3233856
                                                                                                                                                    |                                                                                                                                                  
gggCGCAACCGCTGATTGCTAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGAATGCAACg                                                                                                                                                   <  2:1385123/149‑1 (MQ=255)
   cgcAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                                                                                                >  1:298539/1‑150 (MQ=255)
     cAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATcc                                                                                                                                              <  1:794171/150‑1 (MQ=255)
        ccGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGc                                                                                                                                           <  2:893467/150‑1 (MQ=255)
        ccGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGc                                                                                                                                           >  1:1192785/1‑150 (MQ=255)
          gCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCGGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGcgc                                                                                                                                         <  1:347955/150‑3 (MQ=255)
            tgattgatAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGt                                                                                                                                       >  1:52894/1‑150 (MQ=255)
                 gatAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCTGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGt                                                                                                                                  >  2:206737/1‑150 (MQ=255)
                 gatAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGt                                                                                                                                  >  2:289130/1‑150 (MQ=255)
                    aGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                                                                                                               >  2:343281/1‑150 (MQ=255)
                     gCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAAt                                                                                                                              >  2:47081/1‑150 (MQ=255)
                                 tGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAgg                                                                                                                  >  2:1291941/1‑150 (MQ=255)
                                    aTTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCgg                                                                                                               >  2:450131/1‑150 (MQ=255)
                                      ttAATTGGCGGCTGCGATGTGGCTATGGAGCGGGACGCACGACGGGCAATTGCCCAGGGAACAGGGGTGGCGCTGGAGATTTAATTCGTTATTTGACGCGCGGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGga                                                                                                             <  1:715459/150‑1 (MQ=255)
                                       tAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGATGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGgag                                                                                                            >  2:42617/1‑150 (MQ=255)
                                            ggcggcTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATg                                                                                                       <  2:977255/150‑1 (MQ=255)
                                              cggcTGCGATGTGGCTAGGGAGCTGGCCGCACGACGGGCATTTGCCCAGGGAACACGGCTGGCGCTGGTGATTTAAATCGTTATTTGATGCGCGGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGcc                                                                                                     <  2:409253/150‑1 (MQ=255)
                                               ggcTGCTATGTGGCTTTGGAGCTGGACGCACGAGGGGAAATTGCCCAGGGAACAGGGGTGGGGCTGGATATTTAAATCGTTATTTCAAGCGCCGGATGCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCg                                                                                                    <  1:1423745/150‑1 (MQ=255)
                                               ggcTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTTACGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCg                                                                                                    <  1:550995/150‑1 (MQ=255)
                                               ggcTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCg                                                                                                    >  1:1072196/1‑150 (MQ=255)
                                                    cGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATcg                                                                                               >  2:648434/1‑150 (MQ=255)
                                                      aTGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATcgcg                                                                                             >  1:713831/1‑150 (MQ=255)
                                                       tgtgGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATcc                                                                                                                                              <  1:822723/100‑1 (MQ=255)
                                                       tgtgGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATcc                                                                                                                                              >  2:822723/1‑100 (MQ=255)
                                                        gtgGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGtt                                                                                           <  2:713831/150‑1 (MQ=255)
                                                        gtgGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGtt                                                                                           <  2:1415595/150‑1 (MQ=255)
                                                           gCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGcc                                                                                        >  2:747569/1‑150 (MQ=255)
                                                            cTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCa                                                                                       >  2:581858/1‑150 (MQ=255)
                                                                ggAGCTGGACGCACGAGGGGCAATTGCCCTGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGTATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAttt                                                                                   <  2:1142283/150‑1 (MQ=255)
                                                                gcaGCTGCACGCACGACGGGCAATGGCCCCGGCAACACGGCGGGCGGTGGAGATTTAAATCGTTATTTTACGCGCGGTATGCAACCCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAttt                                                                                   <  2:1312750/148‑1 (MQ=255)
                                                                   gCTGGACGCACGACGGGCAATTGCCCAGGCAACACGGCTGGCGCGGTACATTTAAATCTTTATTTGAAGCGCCGTATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAtttttt                                                                                <  2:1234484/150‑1 (MQ=255)
                                                                     gggACGCACGCCGGCCAATTGCCCAGGCAACAGGGCTGGCGCTGGATATTAAACTCGTTTTTTCAAGCGCGGGATGCAACCCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAg                                                                              <  2:1392043/149‑1 (MQ=255)
                                                                      ggACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGc                                                                             >  1:910666/1‑150 (MQ=255)
                                                                      ggACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGc                                                                             >  1:442339/1‑150 (MQ=255)
                                                                       gACGCACGACGGGCAATTGCCCAGGGAACAGGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCg                                                                            <  2:776039/150‑1 (MQ=255)
                                                                       gACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCGGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAGCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCg                                                                            <  2:411521/150‑1 (MQ=255)
                                                                        aCGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCgt                                                                           >  2:1269888/1‑150 (MQ=255)
                                                                            acgacgGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTc                                                                       >  1:1264338/1‑150 (MQ=255)
                                                                             cgacgGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAACGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGttac                                                                      >  2:1207296/1‑148 (MQ=255)
                                                                               aaggtcAATTTCCCAGGAAACACGTCTCGCTCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAg                                                                    >  2:1104869/6‑150 (MQ=255)
                                                                                cgGGCAATTTCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCAGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGc                                                                   >  2:1316634/1‑150 (MQ=255)
                                                                                 gggCAATTGCCCTGGGAACAGGGCGGGCGGGGGATTTTTAAATCTTTATTTGACGCGCGGTATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCg                                                                  <  1:289130/150‑1 (MQ=255)
                                                                                          cccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc                                                         >  1:105363/1‑150 (MQ=255)
                                                                                          cccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc                                                         <  2:232514/150‑1 (MQ=255)
                                                                                           ccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCct                                                        >  1:1035977/1‑150 (MQ=255)
                                                                                                gAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTcc                                                   >  1:829046/1‑150 (MQ=255)
                                                                                                     cGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTc                                              >  1:98028/1‑150 (MQ=255)
                                                                                                        cTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaaa                                           <  1:581858/150‑1 (MQ=255)
                                                                                                          ggCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAc                                         >  2:868589/1‑150 (MQ=255)
                                                                                                             gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                      <  1:1352700/150‑1 (MQ=255)
                                                                                                             gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                      <  2:1061311/150‑1 (MQ=255)
                                                                                                             gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                      >  2:257197/1‑150 (MQ=255)
                                                                                                              cTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGt                                     <  1:206737/150‑1 (MQ=255)
                                                                                                                 gagaTTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGac                                  <  2:1072196/150‑1 (MQ=255)
                                                                                                                          aTCGTTATTTGGAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGAt                         >  2:2740/1‑150 (MQ=255)
                                                                                                                              ttATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTAc                     <  2:442339/150‑1 (MQ=255)
                                                                                                                              ttATTTCACGCGCCGGAGGCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTAc                     <  1:1291941/150‑1 (MQ=255)
                                                                                                                               tATTTCAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACt                    <  2:290725/150‑1 (MQ=255)
                                                                                                                                aTTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTg                   <  1:42617/150‑1 (MQ=255)
                                                                                                                                  ttGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGcc                 >  1:878749/1‑150 (MQ=255)
                                                                                                                                       gtgcCGGATGCGTTGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAAttt            >  1:601902/3‑150 (MQ=255)
                                                                                                                                       gtgcCGGATGCGTTGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAAttt            >  1:940877/3‑150 (MQ=255)
                                                                                                                                        tgcCGGATGCGTTGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTc           >  1:790201/2‑150 (MQ=255)
                                                                                                                                          ccGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCg         <  1:450131/150‑1 (MQ=255)
                                                                                                                                            ggATGCAACGCACCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc       <  1:1104869/150‑1 (MQ=255)
                                                                                                                                                 cAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATc  >  2:653750/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATC  >  NC_000913/3233562‑3233856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: