Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,708 A→G 6.9% intergenic (+25/+20) fadH → / ← higA 2,4‑dienoyl‑CoA reductase, NADH and FMN‑linked/antitoxinof the HigB‑HigA toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,233,7080AG 6.9% 253.9 / 7.7 87intergenic (+25/+20)fadH/higA2,4‑dienoyl‑CoA reductase, NADH and FMN‑linked/antitoxinof the HigB‑HigA toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (50/31);  new base G (2/4);  total (52/35)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCG  >  NC_000913/3233561‑3233857
                                                                                                                                                   |                                                                                                                                                     
cTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAAc                                                                                                                                                     <  2:72827/150‑1 (MQ=255)
      cAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGCCGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGACGCGCCGGATGCAACGCATcc                                                                                                                                               >  2:1621333/1‑150 (MQ=255)
      cAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATcc                                                                                                                                               <  1:1323858/150‑1 (MQ=255)
         ccGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGc                                                                                                                                            <  2:1318480/150‑1 (MQ=255)
         ccGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCGTCCGGc                                                                                                                                            >  1:229587/1‑150 (MQ=255)
          cGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCa                                                                                                                                           <  1:22674/150‑1 (MQ=255)
           gCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCATTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGcgc                                                                                                                                          <  1:1416215/150‑3 (MQ=255)
           gCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGACCACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCACCGCATCCGGCAc                                                                                                                                          >  2:103475/1‑150 (MQ=255)
           gCTGACTGATAGCCGGAAAACGGTGCATTTAATTGGCGGCGGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTGGCATCCGGcgc                                                                                                                                          <  2:454763/150‑3 (MQ=255)
            cTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGcgct                                                                                                                                         <  1:62083/150‑4 (MQ=255)
                  tatacggggAAAACGGGGCATTTAATTGGCGGCTGCTATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCTGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGt                                                                                                                                   <  1:972631/144‑1 (MQ=255)
                  gatAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGt                                                                                                                                   >  1:1085332/1‑150 (MQ=255)
                      gCGGTAAAACGGGGCATTTAATTGGCGGCTGCGATGTGGCTATGGCGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAAt                                                                                                                               <  1:46814/150‑1 (MQ=255)
                       cGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATc                                                                                                                              <  2:257517/150‑1 (MQ=255)
                       cGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCAATGGAGCTGGACGCACGACGGGCAATTGCCGAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACCCATACGGCACGTCGCGtttaga                                                                                                                              >  1:1683418/1‑146 (MQ=255)
                              aCGGGGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCGGGACGCACGACGGGCAATTGCCCAGGTACCACGGCTGCCGCTGGAGATTTAATGCGTTATTTTAAGCGCCGGATGCAAGGCATCGGGCACGTCGCGTTAATCAATAAAc                                                                                                                       <  2:85871/150‑1 (MQ=255)
                               cGGTGCATTTAATTGGCGGCTGCGCTGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAAGAAACa                                                                                                                      >  1:113748/1‑150 (MQ=255)
                               cGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACa                                                                                                                      >  1:501524/1‑150 (MQ=255)
                                ggTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACaa                                                                                                                     >  2:918607/1‑150 (MQ=255)
                                   gCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGc                                                                                                                  >  2:602521/1‑150 (MQ=255)
                                    cATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCg                                                                                                                 >  2:1230978/1‑150 (MQ=255)
                                      tttAATTGGCGGCTCCGATGTGCCAATGTAGCGGGACGCTCGCCGGGCAATTGCCCAGGGAACCCGGCTGCCGCTGGAGATTTAAATCGTTATTTGAAGGGCCGTATGAAAGGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCggg                                                                                                               <  2:1261007/150‑1 (MQ=255)
                                       ttAATGTGCGGCAGCGACGTGCCTCTCGAGCCGGCCCCACCCCGCGCAATGCCCCAGTTAATATGGCTGGCGCTGCTGATATAAATCTTTATTCCATGCCCCGTATGCACCCCAGCCGGCACTTCGCGTTAATCAATAAACAAGGCGGga                                                                                                              <  2:158545/150‑1 (MQ=255)
                                        tatttGGGGGCTGCTTTGGGGCTTGGGCGCGGGCCGCACGACGGGAATTTCCCCAGGTACCACGGGGGGCGCTGGCGATTTAAAGCGTTATTTGACGCGCCGGAGGCACCGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGgag                                                                                                             <  1:1228517/147‑1 (MQ=255)
                                        tAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGgag                                                                                                             >  1:501734/1‑150 (MQ=255)
                                        tAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGgag                                                                                                             >  1:1475016/1‑150 (MQ=255)
                                        aaaTTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACCCGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGgag                                                                                                             <  2:45403/149‑1 (MQ=255)
                                             ggcggcTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATg                                                                                                        >  1:1171492/1‑150 (MQ=255)
                                                     cGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGAATTAAATCGTTATTTGAAGCGCCGCATCCAACGCATCCTGCACGTCGCGTTAATCAATAAACAAGGCGGGACAAATGCCGATTcg                                                                                                >  2:142371/1‑150 (MQ=255)
                                                        tgtgGCTATGGAACTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATcgcgg                                                                                             >  1:697053/1‑149 (MQ=255)
                                                           ggCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGGTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGc                                                                                          >  1:475786/1‑150 (MQ=255)
                                                             cTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCa                                                                                        >  1:1072279/1‑150 (MQ=255)
                                                              tATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCaa                                                                                       <  1:48012/150‑1 (MQ=255)
                                                                tGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGCAATGCCGAATCGCGTTGCCAAtt                                                                                     >  2:1225819/1‑150 (MQ=255)
                                                                tGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAtt                                                                                     >  1:488750/1‑150 (MQ=255)
                                                                 ggAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAttt                                                                                    >  1:577555/1‑150 (MQ=255)
                                                                 ggAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAAAAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAttt                                                                                    >  1:1466657/1‑150 (MQ=255)
                                                                       ggACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGc                                                                              >  1:86247/1‑150 (MQ=255)
                                                                       ggACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGc                                                                              >  1:519948/1‑150 (MQ=255)
                                                                            cACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGtt                                                                         <  1:948513/150‑1 (MQ=255)
                                                                             aggacgGGCATTTCCCCAGGGAACACGCCTTGCGGGGTAGATTTAAATCGTTATTTTAAGCGCCGGATGCAACGCTTCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTc                                                                        <  2:969173/148‑1 (MQ=255)
                                                                              cgacgGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCACCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTcc                                                                       <  2:872871/150‑1 (MQ=255)
                                                                               gacgGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCa                                                                      >  2:498371/1‑150 (MQ=255)
                                                                                 cgGGAAATTGCCCTGGCAACCCGGCGGGCGCTGTATATTTAAATCGTTATTTGACGCGCCGGATCCAACCAATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGc                                                                    <  1:732366/150‑1 (MQ=255)
                                                                                   ggCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCCTCCGGCACGGCGCGGTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGt                                                                  >  1:1394506/1‑150 (MQ=255)
                                                                                   ggCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTACAGCGt                                                                  >  1:1139949/1‑150 (MQ=255)
                                                                                      aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTAGCGGTTATCAATAAACAAGGCGGGGGAAATGGCGAATCGCGTTGCCAATTTTTTTGCTGGTGGCAGCGTTat                                                               >  1:1771611/1‑149 (MQ=255)
                                                                                         tGCCCAGGTAACACGGCGGCGGCGGTAGATTAAAATCGATATTTGATGCGCGGGATGCACCGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAAttt                                                            <  2:339774/150‑1 (MQ=255)
                                                                                         tGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAAttt                                                            <  2:1072279/150‑1 (MQ=255)
                                                                                           cccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc                                                          >  1:960823/1‑150 (MQ=255)
                                                                                            ccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATTCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATAGCGTTGCCAATTTTTTATCGTGTTCCAGCGTTAATTTCct                                                         >  2:1120122/1‑150 (MQ=255)
                                                                                            ccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCGTCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCct                                                         <  2:229587/150‑1 (MQ=255)
                                                                                             cAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCCCGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATAGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTGcta                                                        >  1:1136617/1‑149 (MQ=255)
                                                                                             cAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCctc                                                        >  2:564049/1‑150 (MQ=255)
                                                                                               gggAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTTATCAATAAACCAGGCGGGAGCAATGCCGCATCGCGTTGCCAATTTTTTAGCGTGTTCAAGCGTTAATTTACTCtt                                                      >  1:1372285/1‑150 (MQ=255)
                                                                                               ggcaaCACGGCTGGCGCGGTATATTAAAACGTTTTTTTCACGCGCGGTATGCACCGAATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCtt                                                      <  1:1707325/147‑1 (MQ=255)
                                                                                                    cacGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCCTACGGAACCTAGCGGTAATCAATAAACAAGGCGGGAGAAATGCCCAATCGCGTTTCCAAATTTTTATGGTGTTCCAGCGTGAATTTCCTCTTccccc                                                 >  1:1376028/1‑148 (MQ=255)
                                                                                                      gggcgggcgGTGGTTATTTAATTGGTTATTTTACGCGCCGTATGCACCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTc                                               <  1:1699429/149‑1 (MQ=255)
                                                                                                       ggCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCa                                              >  1:448102/1‑150 (MQ=255)
                                                                                                          tGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGTACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGGGTTCCAGCGTTAATTTCCTCTTCCCGCTCaaaa                                           >  2:1264002/1‑150 (MQ=255)
                                                                                                             cgcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgc                                        >  1:808744/1‑150 (MQ=255)
                                                                                                              gcTGGAGATTTAAATCGTTTTTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTACCGTGTTCAAGCCTTAATTTCCTCTTCCCGCTCAAAActcg                                       >  2:1001155/1‑150 (MQ=255)
                                                                                                              gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGGTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                       >  1:1065414/1‑150 (MQ=255)
                                                                                                               gtgtagaTTTAATGCTTTATTTGAAGCGCGGGTGTGACCCCTTCGGGCACGTCGCGTTATTCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGt                                      <  2:148395/146‑1 (MQ=255)
                                                                                                               cTGGAGATTTAAATCGTTATTTGAAGCGCCGGATTCAAAGCATCCGGCACGTCTCGTTAATCAATAAACAATGCGGTAGAAATGACAAAAATCGTTTAAAATTTTTTATAGTTTTCCAGCTTTAATTTACTCTTCCCGCTAAAAAAGCGt                                      >  1:408965/1‑150 (MQ=255)
                                                                                                               cTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGt                                      >  2:1205741/1‑150 (MQ=255)
                                                                                                                tGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTg                                     <  1:39049/150‑1 (MQ=255)
                                                                                                                tGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTg                                     <  2:1171492/150‑1 (MQ=255)
                                                                                                                 ggAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAAAAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGa                                    >  2:1261986/1‑150 (MQ=255)
                                                                                                                       ttAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCAAGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCAt                              >  1:1546380/1‑150 (MQ=255)
                                                                                                                        taatgctttATTTGACGCGCCGTTTCCACCGCTGCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATa                             <  2:1395576/143‑1 (MQ=255)
                                                                                                                        tAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAAGCGTGACACCATa                             >  2:855623/1‑150 (MQ=255)
                                                                                                                                 aTTTGACGCGCCGGATCCAACGCATCCGGCACGTCGCGTTAATCAATCAACAAGGCGGTAGAAATGCCGAATCGCGTTGCCAATTTATTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGGGACACCATAGATTTACTg                    <  2:113748/150‑1 (MQ=255)
                                                                                                                                 aTTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTg                    >  2:1083351/1‑150 (MQ=255)
                                                                                                                                  tttGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGc                   <  2:129995/150‑1 (MQ=255)
                                                                                                                                    tGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTACTCTTACCGCTCAAAACGCGTGACACCATAGATTTACTGCCa                 >  2:659871/1‑150 (MQ=255)
                                                                                                                                     tacgcgcCGTATCCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCaa                <  2:730534/147‑1 (MQ=255)
                                                                                                                                     gAAGCGCCGGATGCAACGCATCCGGCACGTCGCGGTAATCAATAAACAAGGCGGGCGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAACCGCGTGACACCATAGATTTACTGCCaa                <  2:296177/150‑1 (MQ=255)
                                                                                                                                       cgtgcCGGATGCGTTGCATCCGGCGCGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTAGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAAtt              >  1:1737573/4‑150 (MQ=255)
                                                                                                                                       aGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAAtt              >  2:692359/1‑150 (MQ=255)
                                                                                                                                          gcCGGATGCGTTGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTcc           >  1:697678/1‑150 (MQ=255)
                                                                                                                                          gcCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTcc           >  1:1761457/1‑150 (MQ=255)
                                                                                                                                           ccGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCg          >  1:1326887/1‑150 (MQ=255)
                                                                                                                                             ggATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc        <  2:1085538/150‑1 (MQ=255)
                                                                                                                                             ggATCCACCGCACCCGGCACGTCGCGTTAATCAATAAACAAGTCCGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATATCCGGc        <  2:519948/150‑1 (MQ=255)
                                                                                                                                              ttcgaaaCTCAGCCGCCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCa       <  2:1626273/145‑1 (MQ=255)
                                                                                                                                                  cAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATc   >  2:585968/1‑150 (MQ=255)
                                                                                                                                                  cAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATc   >  1:1758854/1‑150 (MQ=255)
                                                                                                                                                   aaCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCg  <  1:1205741/150‑1 (MQ=255)
                                                                                                                                                   aaCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCg  >  1:1458783/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCG  >  NC_000913/3233561‑3233857

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