Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,190,053 A→G F132S (TTT→TCT)  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,190,0530AG94.9% 121.7 / NA 39F132S (TTT→TCT) phoPDNA‑binding transcriptional dual regulator PhoP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (19/18);  minor base T (0/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATCGGCAGTGAAACATCGTTGCTACGCCAGCGGCGAATCAGTGACAGACCGTCCTCGTCTGGCA  >  NC_000913/1189810‑1190287
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
aCGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATGAGTGCGTATCAACGTTTCCATAAGAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACTTGATACggggg                                                                                                                                                                                                                                     <  2:697102/251‑1 (MQ=255)
aCGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGAAGACGTCCCATCAGTACATCAATGGTATGGCATTCCCGCAGCTACGCATCCGGATAGAGTTGGCGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTAAATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACggggg                                                                                                                                                                                                                                     >  1:705845/1‑251 (MQ=255)
   gtggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGAATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCCTTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGAATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACggggggcg                                                                                                                                                                                                                                  >  1:98009/1‑249 (MQ=255)
        aaTCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATGAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGaa                                                                                                                                                                                                                              <  1:38001/250‑1 (MQ=255)
           cACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGAAGGGGGGAGCGAAATg                                                                                                                                                                                                                           >  1:378243/1‑250 (MQ=255)
              ttcttGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCt                                                                                                                                                                                                                       >  1:239523/1‑251 (MQ=255)
                   ggggATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGa                                                                                                                                                                                                                   <  1:700549/250‑1 (MQ=255)
                      gATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGcc                                                                                                                                                                                                               <  1:585933/251‑1 (MQ=255)
                              gCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATGCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCt                                                                                                                                                                                                       <  1:346049/251‑1 (MQ=255)
                                 tGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTAtt                                                                                                                                                                                                    >  1:23952/1‑251 (MQ=255)
                                              gcAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAAt                                                                                                                                                                                       <  2:396823/251‑1 (MQ=255)
                                                 gACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGcc                                                                                                                                                                                    <  1:546223/251‑1 (MQ=255)
                                                   cGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTg                                                                                                                                                                                  >  2:282997/1‑251 (MQ=255)
                                                    gTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGc                                                                                                                                                                                 >  1:277320/1‑251 (MQ=255)
                                                     tCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGCTAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGc                                                                                                                                                                                 >  1:438562/1‑250 (MQ=255)
                                                                  tCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACCTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTAGCGCCATCAc                                                                                                                                                                   >  1:109656/1‑251 (MQ=255)
                                                                        gTATGGCTTTCCCGCAGCTCCGCATCCGGTTAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTc                                                                                                                                                             >  2:583798/1‑251 (MQ=255)
                                                                          aTGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACCTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACTGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTTTGAAGCCAGACCGCTACTACGCCGCATTAATGCCTGCATTCGCGCCATCCCCTCTTCaa                                                                                                                                                           <  2:298764/251‑1 (MQ=255)
                                                                                          tcagcagccGGATAGAGTTGGAGCATTAACGACTCTTTGTTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCCTTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGATATTACGCCGCATTAATGCCTGCCTTCGCGCCCTCACCTCTTCAATATGAGACAGTTTatt                                                                                                                                           >  2:399132/8‑249 (MQ=255)
                                                                                                     aTAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAatcat                                                                                                                                <  1:252659/251‑1 (MQ=255)
                                                                                                         agTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGc                                                                                                                            >  1:140678/1‑251 (MQ=255)
                                                                                                             ggAGCAGTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCTGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAAGGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCAcc                                                                                                                         <  1:380156/250‑1 (MQ=255)
                                                                                                                       cGAATCTTTGCTGACCACTGTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTAAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAAt                                                                                                               >  1:119039/1‑250 (MQ=255)
                                                                                                                                          ttGCCCTTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGCGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCACTATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGcc                                                                                            >  2:589310/1‑250 (MQ=255)
                                                                                                                                           tGCCATTATTGCGTATCAACTTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATGCACGGCGAGAGAGATCAACCTGATACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCTCCTCGTCAATATGAAACGGTTGAGTCACATAAGCAGCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCa                                                                                           <  2:698844/250‑1 (MQ=255)
                                                                                                                                               attattGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATAAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCttt                                                                                      >  1:496475/1‑251 (MQ=255)
                                                                                                                                                ttattGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTc                                                                                     >  2:393772/1‑251 (MQ=255)
                                                                                                                                                  attGCGTATCAACGTTTACATAAAAGGTTCTTCAAACGCGGTCCGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGACATGACCTACGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGACACGGTTTAGGAACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCa                                                                                    >  2:374497/1‑250 (MQ=255)
                                                                                                                                                     gCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATGCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCAGTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTGCGACGTTGTCCTGCCAGCTTTCACggg                                                                                <  2:438562/251‑1 (MQ=255)
                                                                                                                                                        tATCACCGCTTCCATACTCGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGCTAATTCACGGCGAGAGAGCTCAACCTGAGACGGGGGGAGCGAAATGAACTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTACTGCCAGCTTTCACGGGCgg                                                                             >  2:327659/1‑251 (MQ=255)
                                                                                                                                                         aTCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACGTCGACGTTGTCCTGCCAGCTTTCACGGGCGGt                                                                            <  2:140678/251‑1 (MQ=255)
                                                                                                                                                                   ccATAATAGTGTATTCGAACGCGGGCAGTTTGATCACGTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTGGTCCTGCCAGCTTTCACGGGCGGTTAATACCAg                                                                   <  2:510834/250‑1 (MQ=255)
                                                                                                                                                                       aataGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATc                                                               <  1:194393/250‑1 (MQ=255)
                                                                                                                                                                                             aGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCATCCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATCGTCAGTGAAACATCGTTGCTACg                                        <  2:586521/251‑1 (MQ=255)
                                                                                                                                                                                                     cACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATCGGCAGTGAAACATCGTTGCTACGCCAgcggc                                <  1:583798/251‑1 (MQ=255)
                                                                                                                                                                                                          cGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATCGGCAGTGAAACATCGTTGCTACGCCAGCGGCGAAt                            <  1:336051/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCGCGGGCGGTTAATACCAGAATCGGCAGTGAAACATCGTTGCTACGCCa                                     >  1:251579/1‑218 (MQ=255)
                                                                                                                                                                                                                                 cGAGAGAGATCAACCTGAGACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATCGGCAGTGAAACATCGTTGCTACGCCa                                     <  2:251579/218‑1 (MQ=255)
                                                                                                                                                                                                                                   agagagagCAACCTGAGACGGGGGGAGCGAAATTACCTGTGAAGCCAGACCGCTATTACGCCGCATTAAGGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATCGGCAGTGAAACATCGTTGCTACGCCAGCGGCGAATCAGTGACAGACCGTCCTCGTCTGGCa  <  1:117000/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
ACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATCGGCAGTGAAACATCGTTGCTACGCCAGCGGCGAATCAGTGACAGACCGTCCTCGTCTGGCA  >  NC_000913/1189810‑1190287

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