Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,189,980 A→G T156T (ACT→ACC phoP ← DNA‑binding transcriptional dual regulator PhoP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,189,9800AG97.7% 150.4 / NA 44T156T (ACT→ACCphoPDNA‑binding transcriptional dual regulator PhoP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/20);  minor base T (1/0);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACC  >  NC_000913/1189739‑1190220
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
aGCGAGAGCGGTACAAAACGACGCAGTAATTTTTTCATCAGCGCAATTCGCACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGCGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTCTCAACGTTTCCATAATTGGGTATTc                                                                                                                                                                                                                                           >  2:347984/1‑250 (MQ=255)
        cGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGt                                                                                                                                                                                                                                  >  1:3340/1‑251 (MQ=255)
          ggAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTc                                                                                                                                                                                                                                 >  1:489767/1‑250 (MQ=255)
                       cAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTc                                                                                                                                                                                                                   >  2:120687/1‑251 (MQ=255)
                            aTTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCAt                                                                                                                                                                                                              >  1:389853/1‑251 (MQ=255)
                                 ttCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAAt                                                                                                                                                                                                          <  1:436897/250‑1 (MQ=255)
                                      caGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATa                                                                                                                                                                                                     <  1:70042/250‑1 (MQ=255)
                                      caGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGGACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGAGTATCAACGTTTCCTTAATGGTGTATTCGAACGCGGTCAGTTTGCTCACTTCGTCATTAATAGATa                                                                                                                                                                                                     >  2:54916/1‑250 (MQ=255)
                                       aGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACTAATCTTTGCTGACAGCTTTGCCATTATTGCGAATCAACGTATCCATAATGGTGTGTTCGAACGCGGTCAGTTTGGTCACTTCGCCATTAATAGATAAt                                                                                                                                                                                                   >  2:103417/1‑251 (MQ=255)
                                         cgcAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTc                                                                                                                                                                                                 <  1:169905/251‑1 (MQ=255)
                                          gcAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCAGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAAttac                                                                                                                                                                                                >  2:15361/1‑249 (MQ=255)
                                                  aaCAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCgg                                                                                                                                                                                                                                   <  1:151741/208‑1 (MQ=255)
                                                  aaCAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCgg                                                                                                                                                                                                                                   >  2:151741/1‑208 (MQ=255)
                                                  aaCAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCgaga                                                                                                                                                                                        <  2:464056/251‑1 (MQ=255)
                                                        atgcccgtgcctcgaccGGTGGTAATCACTTCTTGGGGCTATTGCGCTTGTATTTGTTTGCGCATACGTCCAAACAGTACATCAAGGGTATTGCTTTCCAGCAGCTCCGCATCCGGAGAGAGTTGGAGCATTAACGAATCTTTGCTTACCACTTTGCCATTATTGCGTATCAACGTTTCCTTAATGGTTTATTCGAACTCGGTCATTTTGATCACTTCGTCATTACTAGATATTTCACGGTGCGAGAGATc                                                                                                                                                                                  <  2:487731/235‑1 (MQ=255)
                                                                 ccGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGCGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCCCGGCGAGAGAGATCAACCTGac                                                                                                                                                                          >  2:222587/1‑249 (MQ=255)
                                                                 ccGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGaaa                                                                                                                                                                         <  1:38376/251‑1 (MQ=255)
                                                                      aaCGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGAATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACgggg                                                                                                                                                                    >  2:292973/1‑251 (MQ=255)
                                                                                   aCTTCTTGGGGATATTGTGCCTGAATTTTTTTGAGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGAc                                                                                                                                                       >  1:194367/1‑251 (MQ=255)
                                                                                              atatTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACAGCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGATAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCa                                                                                                                                            <  2:157656/251‑1 (MQ=255)
                                                                                                                           gTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAACTGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGc                                                                                                               >  2:164777/1‑251 (MQ=255)
                                                                                                                                       cATCACTGGTATGGCTTTACCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATc                                                                                                   >  1:159549/1‑251 (MQ=255)
                                                                                                                                       cATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGTGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATc                                                                                                   <  2:49795/251‑1 (MQ=255)
                                                                                                                                              ggTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCATTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCtt                                                                                           <  2:3340/251‑1 (MQ=255)
                                                                                                                                                 aTGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCaa                                                                                        <  1:49686/251‑1 (MQ=255)
                                                                                                                                                      tttCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATg                                                                                    >  2:126429/1‑250 (MQ=255)
                                                                                                                                                          ccGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACg                                                                               <  2:480052/251‑1 (MQ=255)
                                                                                                                                                                         cGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGTGAGCGAAATGACCTTTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAat                                                                >  1:381598/1‑251 (MQ=255)
                                                                                                                                                                         cGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGAGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGTGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAat                                                                >  1:489975/1‑251 (MQ=255)
                                                                                                                                                                                                aaTCTTTGCTGACCTCTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACt                                         <  1:347984/251‑1 (MQ=255)
                                                                                                                                                                                                aaTCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACt                                         <  1:292973/251‑1 (MQ=255)
                                                                                                                                                                                                aaTCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCAATTACGCCGCCTTATTGCCTGCATTCGCACCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTacaact                                         >  2:462189/1‑246 (MQ=255)
                                                                                                                                                                                                     ttGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACtt                                        <  1:126429/247‑1 (MQ=255)
                                                                                                                                                                                                        cTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACttt                                 >  1:71876/1‑251 (MQ=255)
                                                                                                                                                                                                          gACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACttt                                 <  1:292361/249‑1 (MQ=255)
                                                                                                                                                                                                          gACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACttg                                 <  2:475028/249‑2 (MQ=255)
                                                                                                                                                                                                               cTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGc                          >  1:423685/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGcc                         <  1:27485/250‑1 (MQ=255)
                                                                                                                                                                                                                  tGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGcc                         >  1:411107/1‑249 (MQ=255)
                                                                                                                                                                                                                    ccATTATTGCGTATCAACTTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCt                     <  1:25870/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCt                     <  1:393136/251‑1 (MQ=255)
                                                                                                                                                                                                                                aTCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGAAAATTCACGGCGAGAGCGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGCCCGCTATTAAGCCGCATTAATGCCTGCATTCGCGCCATACACTCTTCAATATGAAACGGCTTAGTCACATAATCATCAGCCCCGGCACTTAATACTACGACTTTGTCCTTCCAGCTTTCACGGGCGGt         >  2:158773/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aaCGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTa       >  1:286320/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ccATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATCATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATAcc  >  1:226088/1‑248 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
AGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACC  >  NC_000913/1189739‑1190220

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