Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,708 A→G 7.6% intergenic (+25/+20) 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,7080AG 7.6% 194.1 / 5.5 66intergenic (+25/+20)fadH/higA2,4‑dienoyl‑CoA reductase/antitoxin/DNA‑binding transcriptional repressor HigA
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (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 = 5.72e-01

GCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATC  >  NC_000913/3233560‑3233856
                                                                                                                                                    |                                                                                                                                                    
gcTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCaa                                                                                                                                                     >  2:1352700/1‑150 (MQ=255)
  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:442339/1‑150 (MQ=255)
                                                                        ggACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGc                                                                             >  1:910666/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                                                         <  2:232514/150‑1 (MQ=255)
                                                                                            cccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc                                                         >  1:105363/1‑150 (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                                      >  2:257197/1‑150 (MQ=255)
                                                                                                               gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                      <  2:1061311/150‑1 (MQ=255)
                                                                                                               gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                      <  1:1352700/150‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
GCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATC  >  NC_000913/3233560‑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: