Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,449,045 T→G 26.8% intergenic (+3/+77) chpB → / ← ppa toxin of the ChpB‑ChpS toxin‑antitoxin system/inorganic pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,449,0450TG26.8% 111.6 / 11.6 64intergenic (+3/+77)chpB/ppatoxin of the ChpB‑ChpS toxin‑antitoxin system/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (11/34);  new base G (17/2);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 2.59e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTAATCAACTGGGAATGACGCTGGTGGCCCCCATTACGCAGGGCGGAAATTTTGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTAA  >  NC_000913/4448808‑4449282
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
cTTTAATCATCTGGGAATGCCGGGGGTGGCCCCCATTACGCAGGGGGGAAATTTTTGCCGTTTTGCCGGATGTGTCGTTCCTTTACATTGCGAAGAAGGCGATGTTCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACTCGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGgtg                                                                                                                                                                                                                                                       <  2:207289/239‑1 (MQ=255)
cTTGAATCAACGGGGAATGAACGTGGCGGCCCCCATTACGCAGGGCGGAAATTTGGCCCGTTATGCCGTATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGgt                                                                                                                                                                                                                                                        <  1:266524/238‑1 (MQ=255)
          cTGGGAATGACGCTGGTGGCCCCCATTACGCAGGGCGGAAATTTTGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGGGGGGGAATAAGGGGGGtttatttt                                                                                                                                                                                                                                            >  1:1161083/1‑238 (MQ=255)
               aaTGACGCTGGTGGCCCCCATTACGCAGGGCGGAAATTTTGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGGGGAATAAGGGGGGGTTATTTTTtgcg                                                                                                                                                                                                                                       >  1:436424/1‑239 (MQ=255)
                    cGCTGGTGGCCCCCATTACGCAGGGCGGAAATTTTGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCGGGCTGCGGATGAGGTGGTGGAAGAGGCGTTTTTACGCTTGCAGGGGGGGGGGGGAAAAGGGGGGttttttttttgttttttt                                                                                                                                                                                                                                  >  2:744453/1‑228 (MQ=255)
                    cGCTGGGGGCCCCCATTACGCAGGGCGGAAATTTGGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCAGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCa                                                                                                                                                                                                                                   <  1:1047857/238‑1 (MQ=255)
                                 cctttagCAGGGCGAAAATTTTGCCCGGTATGGCGGATGTGGCGTTCCTTTACATTGCGGAGGAGGCGATGTGGCAGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTa                                                                                                                                                                                                                    <  2:1272489/233‑1 (MQ=255)
                                   ttACGCAGGGCGGAAATTTTGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGGTTTTTTTTTGCGGGGATAAAAAAAACCCTACt                                                                                                                                                                                                                  >  1:379039/1‑239 (MQ=255)
                                           ggCGGAAATTTTGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGGGGAATAAGTGGGGTTTATTTTTCGCGGGGATAAAAAAAACCTTCCTAACCGaa                                                                                                                                                                                                           >  1:683674/1‑238 (MQ=255)
                                           ggCGGAAATTTTGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGGGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGAGGGGGTATAAGGTGTGTTTTTTTATCGCGGGGATAAAAAAACCCTTCTTAACCGaa                                                                                                                                                                                                           >  1:532110/1‑238 (MQ=255)
                                           ggCGGAAATTTTGCCCGTTATGCCGAATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGCGGTGCTTGTGAATCAGGTGCGGCGTAGGGATCTACACGCCCGGCTGGCAGAGCGTATTGATCTGGCTGCGGTTGCGGTGGTGGAAGGGGTGCTATTACGCTTCTAGGCGGGTGGGGAACAAGGGGGGTTGTTTTATCCCGGGCAAAAACAAACCCTTACTAAcagcac                                                                                                                                                                                                          >  2:344707/1‑234 (MQ=255)
                                                aaaTTTTGGCCGTTTGGGCGGATTTTGGTTTCCTTTACAGTGCGAAGAAGGGGATGTGCACGGCGGGGGTCTGGTGAATCAGGTGCGGATGATGGAACTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCg                                                                                                                                                                                                      <  2:460502/238‑1 (MQ=255)
                                                  aTTTTGCCCGTTAGGCCGGGTGTGAGGTTCCTTTACATGGCGAAGGAGGGGATGTGCAGGGCGTGGTGCTGATGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATGGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCg                                                                                                                                                                                                  <  2:831553/239‑1 (MQ=255)
                                                    tttGGCCGTTAAGCCGGATTTAGCGGTCCTTTACATTGCGAAGAAGGCGAATTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGtt                                                                                                                                                                                                <  1:147696/239‑1 (MQ=255)
                                                    tttGCCCGTTTGGCCGGATTTGGGGGTCCTTTACATTGCGAAGAAGGCGATGTGCCGGGGGTGGGGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGATGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGtt                                                                                                                                                                                                <  1:1253807/239‑1 (MQ=255)
                                                               tGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGttatta                                                                                                                                                                                    <  1:751029/239‑1 (MQ=255)
                                                                         aGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGCGGGGGAATCAGGTGTGTTTATTTATCGCGGGCATAAAAAAAACCTTACTCACCGAAGCCCCGCGTTCAGGGTTATTACGCCaacaga                                                                                                                                                                          >  1:453331/1‑234 (MQ=255)
                                                                          gggTTCCTTTATATGGGGAAAAAGGAGTTGTGCACGGGGGGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACGCGCCCGGCTGGCAAAGCGTATTGGGCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAgaagaa                                                                                                                                                                         <  2:436424/237‑1 (MQ=255)
                                                                           acagctttttCATTGCGGAGAAGGGGTTGGGCGGGGGGTGGTGCTGGTGAAGCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAgaagaa                                                                                                                                                                         <  2:417601/232‑1 (MQ=255)
                                                                            gTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCTGGCGTTCAGGGTTATTACGCCAgaagaa                                                                                                                                                                         <  1:757714/237‑1 (MQ=255)
                                                                              tCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGACTAAGGTGGGTTTATTTATCGCGGGCATAAAAAAAACCTTACTAACCGAAGCCCGGCGGTTAGGGTTATTTCGCCCGAAGGAAtta                                                                                                                                                                     >  1:37700/1‑239 (MQ=255)
                                                                                          cGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGGGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCACGGTTATTACGCCAGAAGAACTTATTTTTTcccg                                                                                                                                                          >  1:417601/1‑234 (MQ=255)
                                                                                          cGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGGGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGGTATTACGCCAGAAAAACTTTTTTATTTTTTgc                                                                                                                                                        >  1:859992/1‑239 (MQ=255)
                                                                                          cGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGCGGTGGAATAAGGGGGGTTTTTTTATCTCGGGCTTAAAAAAAACCCTTATAACGAGAGCCCGGGGTTCAGGGGGTTTAAGCCAAAAAAAAATATTTTATTTTTTg                                                                                                                                                         >  2:1177098/1‑239 (MQ=255)
                                                                                          cGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGAGGGGGGCTAAGGGGGGGTTATTTGTAGCGGGCATAAAAAAACCCCTACTAAACGGAAGCCGGCAGTTAGGGGTATTACGCCAGAAGAACATTTTTAtttttttt                                                                                                                                                        >  2:853296/1‑237 (MQ=255)
                                                                                           gaagaaGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACCCGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTgcg                                                                                                                                                       <  1:119414/239‑1 (MQ=255)
                                                                                            aagaagGCGATGTTCCCGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTgcg                                                                                                                                                       <  1:1131435/238‑1 (MQ=255)
                                                                                            aagaagGCCATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTgcg                                                                                                                                                       <  1:434477/238‑1 (MQ=255)
                                                                                                    gATTTGCACGGTGTTGGTGTGGTTATTCAGGTGCGGTTGAGGGTTCTACCCGCGCGGCTGGCAAAGCGTATTGGGCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGaa                                                                                                                                               <  2:162048/238‑1 (MQ=255)
                                                                                                    gATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGGTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCCCGAag                                                                                                                                              >  1:722263/1‑239 (MQ=255)
                                                                                                           aCGGCGGTGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCaa                                                                                                                                       <  1:796168/239‑1 (MQ=255)
                                                                                                              gCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGGCTGGCTGCGGATGAGGTGGAGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGGGGGATCAGGTGTGTTTCTTTATTGCGGGCATAAAAAAACCCTTACTAACCGAAGACCGGCGGTTGAGGGTGTTACGCCCGAAAGAATTATTTCTTTTTCTCGACGACCAAGGAGGCAAcca                                                                                                                                    >  2:410533/1‑237 (MQ=255)
                                                                                                               cGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAgaagaa                                                                                                                                                                         <  1:1103602/202‑1 (MQ=255)
                                                                                                                          gTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGGGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGGGGTGGAAAAGGGGGGGTTTTTTTTTTGCGGGGTTAAAAAAAACCTTGTTAACCGAGAACCGGCGGTTCGGGGTGTTACGACAGAAGACCTTTTTTTTTATTTGTGCGCGCGAAGGAGGCAAAGATATGAGGTTTa                                                                                                                         >  2:669742/1‑239 (MQ=255)
                                                                                                                            gAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAg                                                                                                                        <  1:669742/238‑1 (MQ=255)
                                                                                                                               tCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCCGAAGAACTTATTTATTCTTTGCGCGCTCGAAAGAGGCCACGATTTCAGCTTTAGCg                                                                                                                      >  1:301831/1‑237 (MQ=255)
                                                                                                                               tCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGGGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCCGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCg                                                                                                                      >  1:705598/1‑237 (MQ=255)
                                                                                                                                           ggtggattctcacGCCCGGTGGGCAAGCGTTTGTAGTCTGGCGGGGGGTGTGGGGGTGGAATGTGTCTTATTGCGCTTACGGGGGGTGGGGAAGAAGGTGTGTTTATTGATCGCGGTCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTGTAGCGGCTTCTGCGTTTTcc                                                                                                     <  2:239938/229‑1 (MQ=255)
                                                                                                                                             tgGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTcc                                                                                                     <  1:853296/238‑1 (MQ=255)
                                                                                                                                                aTCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCg                                                                                                                                                                                                      <  1:54323/143‑1 (MQ=255)
                                                                                                                                                   tACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCACGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGtt                                                                                                                                                                                                <  1:93038/145‑1 (MQ=255)
                                                                                                                                                   tACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCACGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGtt                                                                                                                                                                                                >  2:93038/1‑145 (MQ=255)
                                                                                                                                                    acacGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTc                                                                                            <  1:978331/239‑1 (MQ=255)
                                                                                                                                                       cGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGGGGAATAAGGGGGGTTTATTTTTCGCGGGGCTAAAAAAAACCTTACTAACCGAAGCCCGGCCTTTAGGGTTATTTACCCCAAAGAAATTCTTTTTTTTTTGGCCGCTCGAAGGGGGCCACACTTTCACATTTAGCGGCTTTTTTGTTTTTCCCAACCtgaaa                                                                                           >  2:1129216/1‑235 (MQ=255)
                                                                                                                                                                     tgggttttGGTCTGGGTGCGGATTGGGTGGTGGAAAGGGCGTTATTACGCTTTCAGGGGGTGGTGGAATGTGGTGTTGTTATTTATCGCGGGCATAAAAAAACCCTTTCTAACCGAAGCCCGGCGTTCAGTGTTATTACGCCAGTAGATCTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGc                                                                          <  2:680115/233‑1 (MQ=255)
                                                                                                                                                                             tGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCTTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTc                                                                   <  1:602695/239‑1 (MQ=255)
                                                                                                                                                                               ggcTGGCTGCGGATGAGGTGGTTGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGa                                                                 <  1:410533/237‑1 (MQ=255)
                                                                                                                                                                               gTCTTGGTGCGGGTGAGGGGGTGGAAGAGGCGTTATTTCGCTTGCAGGCGGTGGTTGAAAAAGGGGTGTTTATTTACCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGa                                                                 <  2:301831/239‑1 (MQ=255)
                                                                                                                                                                               gTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGa                                                                 <  1:340912/239‑1 (MQ=255)
                                                                                                                                                                                tCTGGCTGCGGAGTAGGTGGTGGAAGAGGCGTTATTAAGCTTGCAGGCCGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGa                                                                 <  1:1016244/238‑1 (MQ=255)
                                                                                                                                                                                  tGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGGATAAGGGGTGTTTATTTATCGCGGGCATAAAAAAAACCTTACTAACCGAAGCCCGGCGTTTAAGGGTATTACGCCCGAAAGAATTCTTTATTCTTTGTGCGCTCGAAGGAGGGAACGGTTTCAGGTTTTGCGGCTTTTGCGTTTTTCCCACCCTCCACTTTTCCCCACCTTGCCttttctgggg                                                               >  2:1084253/1‑233 (MQ=255)
                                                                                                                                                                                  tGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTTCCAACCCTCAAATTTTACCCCCTTGCCctttttgcggg                                                              >  1:837617/1‑233 (MQ=255)
                                                                                                                                                                                  tGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGGGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAAACGAAGGCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTTCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTTTGCGTTTTTCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGt                                                              >  1:1131858/1‑239 (MQ=255)
                                                                                                                                                                                    gCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGGGGGGTTATTTATTGCGGGCATAAAAAAACCCTTACTAACCGAAGACCGGCGTTTAGGGGTTTTACGCCCGAAGAACTTCATTATTCTTTGCGCGCTCGACGGAGGCAACGATTTCAGCTTTATCGGCTTTTCCGGTTTTCCCACCCTCAAAATTTCACCCCCTCCCCttttttccgtct                                                             >  2:285568/1‑232 (MQ=255)
                                                                                                                                                                                              aGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGttat                                                                                                                                                                                      <  2:229448/112‑1 (MQ=255)
                                                                                                                                                                                              aGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAAACGAAGCCCGGCGTTCAGGGttat                                                                                                                                                                                      >  1:229448/1‑112 (MQ=255)
                                                                                                                                                                                                         gagGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGGGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGGGGCAACGATTTCAGCTTAACGGCTTCTGCGGTTTTCCAACCTTCCAACTTTAACCACTTTCCCTTTTTCCGGTCTTTTTTATGCTCgtagaaagtag                                      >  1:587373/1‑234 (MQ=255)
                                                                                                                                                                                                                       cTTGCAGGCGGTGGTGGAATAAGGGGTGGTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGCAGTTTAGGGTTATTTAGCCAGAAGAAATTATTTATTCTTTGCGCGCCCGAAAGAGGAAACCATTTTCCCTTTAGTGGCTTTTCGGTTTTTCCACCCTCTACCTTTTCCCCCCTTGCCTTTTTTTCGGGTTTTGTATTGTTCGAGAAAGTGAGAGGTTTGTGCttt                         >  2:875592/1‑239 (MQ=255)
                                                                                                                                                                                                                         tggcggCGGGGGTGGAATTAGGTGTGTTTATTTTTCGCGGGCCAAAAAAACCCCTTACTAACCGAAGCCCGGCGTGTAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTca                       <  2:705598/235‑1 (MQ=255)
                                                                                                                                                                                                                         tctgggggggggtggAATAGGGGTGTTTTATTTATGCGGGGCATAAAAAACCCCTCTCTAACCGAATCCCGGCGTTCAGGGTTATTTCGCCAGAAGACTTTATTTATTCTTTGCGCGCTCGAAGGAGGGAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTGGACTTTTTCGAGGGCGTTGTAGTGCTCGAAGAAGTGCGCGATTTGCGCTTTca                       <  2:859992/229‑1 (MQ=255)
                                                                                                                                                                                                                         tGCAGGCGGTTGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTca                       <  1:620757/239‑1 (MQ=255)
                                                                                                                                                                                                                         tGCAGGCGGTGGTGGAATAAGGTTGTTTTAATTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTca                       <  1:24615/239‑1 (MQ=255)
                                                                                                                                                                                                                          gCAGGCGGTGGTGGAATAAGGGGTGTTTATTTATCGCGGGCATAAAAAAACCCTTAATAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTTCGCGCACGAAGGAGGAAACACGTTCAACTTTTGCGGCGTCTCCGGTTTTCCAACCTTCAACTTTCACCCACATGCATTTTTTCAGGGTTTTTTTGGGCTCCAAGAAGTGGGGGAGATTCGCTTTcag                      >  2:1148520/1‑239 (MQ=255)
                                                                                                                                                                                                                                ggtggtGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGtt                 <  1:1014670/238‑1 (MQ=255)
                                                                                                                                                                                                                                      ggAATAAGGTTTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCa           <  1:1153999/238‑1 (MQ=255)
                                                                                                                                                                                                                                      ggAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCa           <  1:1177098/238‑1 (MQ=255)
                                                                                                                                                                                                                                              gtgtgtTTTGTTTTCGCGTGCATAAACAACCCCTTATGAACCGAAGCCCGGCGTTCAGTGTTATTTGGCCAGATGAACTTATTTATTTGTTGCGCTCTCGAAGGAGGCAACGATTTCAGGTTTAGAGGTTTCTGCGTTTTGCCAACCGTCAACGTTCATCCAATTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTaa  <  2:236620/239‑1 (MQ=255)
                                                                                                                                                                                                                                              gtgtgtTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGGCCTCGAAGCAGTGAGCGCTTTGCGCTTTCAGCAGTTCAGGCAGATCGTTa   >  1:620811/1‑238 (MQ=255)
                                                                                                                                                                                                                                              gtgtgtTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAAATTTCACCCACTTGCCTTTTTCGCGGGCTTTTTAGTGCTCGAAGAAGTGAGCGATTTGAGCTTTCAGCAGTTCAGGCAGATCGTTa   >  1:1117493/1‑238 (MQ=255)
                                                                                                                                                                                                                                              gtgtgtTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTATTCGTTTTCCCAACGCTAAAATTTCACCCCCTTGACTTTTTCGGGGCCTTTGTAGGGCTAGCAGAAGTGGAGGCTTTTCCCTTTATACAGTTCAGGCAGAACGTTaa  >  2:1201205/1‑239 (MQ=255)
                                                                                                                                                                                                                                               tgtgtTTTTTTACCGAGGTCATAAAAAAACCCTTCCTAACCGAAGCTCGGGGTTCAGGGTTATTACGCCAGAAGAATTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGGTTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCCCCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTaa  <  2:379039/238‑1 (MQ=255)
                                                                                                                                                                                                                                               tgtgtTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGGTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTaa  <  1:240609/238‑1 (MQ=255)
                                                                                                                                                                                                                                               tgtgtTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTaa  <  1:1010663/238‑1 (MQ=255)
                                                                                                                                                                                                                                               tgtgtTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTaa  <  1:1084253/238‑1 (MQ=255)
                                                                                                                                                                                                                                               tgtgtTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTaa  <  1:420411/238‑1 (MQ=255)
                                                                                                                                                                                                                                               tgtgtTTAATTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTa   <  1:316322/237‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
CTTTAATCAACTGGGAATGACGCTGGTGGCCCCCATTACGCAGGGCGGAAATTTTGCCCGTTATGCCGGATTTAGCGTTCCTTTACATTGCGAAGAAGGCGATGTGCACGGCGTGGTGCTGGTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCACCCACTTGCCTTTTTCGAGGTCTTTGTAGTGCTCGAAGAAGTGAGCGATTTGCGCTTTCAGCAGTTCAGGCAGATCGTTAA  >  NC_000913/4448808‑4449282

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