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,9800AG89.7% 129.2 / NA 39T156T (ACT→ACC) phoPDNA‑binding transcriptional dual regulator PhoP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/13);  minor base T (0/4);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 2.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAA  >  NC_000913/1189740‑1190224
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
gCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCg                                                                                                                                                                                                                                             >  1:360934/1‑250 (MQ=255)
            aaaaaaGACGCAGGAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGATGGTAATCACTTCTTGGGGATCTTGCGCCTGAATTTTTTTGCGCAGACGTCACCTCAGTACATCAATGGTATGGTTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGACATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGACCGCGGTCAGtt                                                                                                                                                                                                                                >  2:48232/1‑251 (MQ=255)
             aaaaaGACGCAGTACTTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGttt                                                                                                                                                                                                                               >  2:80095/1‑251 (MQ=255)
                               tttCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAAt                                                                                                                                                                                                             >  1:546406/1‑251 (MQ=255)
                                           aaTTCGAACAGATAGACCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGAGAAGTCACg                                                                                                                                                                                                 <  1:192790/251‑1 (MQ=255)
                                               cGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCg                                                                                                                                                                                              <  1:235885/250‑1 (MQ=255)
                                                     gATAGACCTGCCCGAGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGCACATCCATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTACCGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCgagag                                                                                                                                                                                          >  2:253116/1‑248 (MQ=255)
                                                      aTAGCCCTGGCCGCGAACGGAGGTAATCACTTCTTGGGGAGATTGTGCCTGAATTTTTTTGCGCAGACGTCACATCAGTACATCAATGGTATGGCTTTCCAGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATATTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCgagag                                                                                                                                                                                          >  2:291564/1‑247 (MQ=255)
                                                             tGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGACTTTTTTTGCGCAGACGTCCCATCCGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTg                                                                                                                                                                               >  1:245218/1‑251 (MQ=255)
                                                             tGGCCGCGAACGGCGGTCCTCACTTCTTGAGGATCTTGCGCCTGAATTTTTTTGCGCCGACGTCCAATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTg                                                                                                                                                                               >  2:76622/1‑251 (MQ=255)
                                                                           ggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGc                                                                                                                                                                  <  2:199390/250‑1 (MQ=255)
                                                                                                       tGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTAtt                                                                                                                                     >  2:10575/1‑251 (MQ=255)
                                                                                                        gAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTAtt                                                                                                                                     <  1:299541/250‑1 (MQ=255)
                                                                                                        gAATTTTTTTGCGCAGACGTCACATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTACATAATGGTGTATTCGAACGCGGTAAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTAtt                                                                                                                                     >  2:240984/1‑250 (MQ=255)
                                                                                                         aaTTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAAt                                                                                                                                                                                                      <  2:427895/184‑1 (MQ=255)
                                                                                                              ttttGCGCAGACGTCCCATCAGTTCATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGATCATTAACGAATCTTTGCTGACCACTTTGCCATGATGGCGTATCAACGTTTCCATAATTGTGTATTCGAACGCGGTCAGTTTGATCTCTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGGGAAGCCAGACCGCTATTAcgccg                                                                                                                               <  2:295893/250‑1 (MQ=255)
                                                                                                              ttttGAGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCAGGAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTCcgccgc                                                                                                                              >  1:369143/1‑251 (MQ=255)
                                                                                                                    gctttctgcctaTCAGTTCATTAATGTTATTGCTTTCCCGCAGCGCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCTCGTTGCCATTATTTCGTCTCAACTTTTCCATAATGGTGTATTCGCACGCGGTCAGTTTGATCGCCTCGTCATTAATAGATGATTCACGGCGAGAGAGAGCAACCTGAAATGTTGGGAGCGAAATGACCTATGTCGCCATTCTGCTATTCCGTCGCATTAAt                                                                                                                        <  2:492574/240‑1 (MQ=255)
                                                                                                                          gTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGCCCGCTATTACGCCGCATTAATGCCTGc                                                                                                                  >  1:53821/1‑251 (MQ=255)
                                                                                                                                      cATCAAGGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATTGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATc                                                                                                      <  1:268715/251‑1 (MQ=255)
                                                                                                                                              gTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTc                                                                                              >  2:499321/1‑251 (MQ=255)
                                                                                                                                                     tttCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATg                                                                                        >  1:200036/1‑250 (MQ=255)
                                                                                                                                                         ccGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACAACTTTGCCATTATTCCGTATCAAAGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACg                                                                                   >  1:356101/1‑251 (MQ=255)
                                                                                                                                                                 ccGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCTCTTTGCCATTATTGCGTATCAACGTGTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTGCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATTCCTGCAGTCGCGCCATCTCCTCTTCAATATGAAACGGTTTAg                                                                             <  2:356101/249‑1 (MQ=255)
                                                                                                                                                                  cGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTca                                                                          <  2:200036/251‑1 (MQ=255)
                                                                                                                                                                     aTCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGATTAGTCACAt                                                                       >  1:435915/1‑251 (MQ=255)
                                                                                                                                                                         ggATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACGTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAatc                                                                   <  1:80095/251‑1 (MQ=255)
                                                                                                                                                                                           aaCGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTa                                                  >  2:169271/1‑250 (MQ=255)
                                                                                                                                                                                             cGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAAt                                                >  1:558324/1‑250 (MQ=255)
                                                                                                                                                                                              gAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGCACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATAc                                              >  2:204526/1‑251 (MQ=255)
                                                                                                                                                                                                         gACCATTTTGGCGTTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTAACGGCGAGAGAGCTCAACCTGACACGGGGGGAGCGAAATGACCTTTGAAGCCAGACCGCTATTAAGCCGCATTAATACCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGt                                   >  2:273420/1‑251 (MQ=255)
                                                                                                                                                                                                                  gCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAg                           <  1:169271/250‑1 (MQ=255)
                                                                                                                                                                                                                    cATTATTGCGTATCAACGTTTCAATAATTGTGTATTCGAACGCGGTCAGTTTGAGCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGAGCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCAATCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCTCATAATCATCAGCTCCGACACGTAATACTTCGACTTTGTCCTGCCAGCtt                        <  2:405413/251‑1 (MQ=255)
                                                                                                                                                                                                                        attGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCAc                    <  1:10575/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTa           <  1:150305/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaCGTTTCCAGAATTGTGTATTCGAACGCGCTCAGTTTGATCACTTCGTCATTAATAGATAATTCACTGCGAGATAGATCAACCTGAAACGGGGGGAGCGAAATGACCTTTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAAGCATCAGCACCGGCACTTAATACTTCGTCTTTGTCCTGCCAGCTTTCACTGGCGGTTaa          <  1:436325/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aCGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGCTTAAt         >  2:43372/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGTTTCCATAATGGAGTATTCGACCGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCAATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAAt         >  1:163539/1‑250 (MQ=255)
                                                                                                                                                                                                                                          cATAATGGTGTATTCGAACGCGGTCAGTTTGAACACTTCTTCATTAATAGATAATTCACGGCGAGAGAGATCACCCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGaa  <  2:19314/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
GCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAA  >  NC_000913/1189740‑1190224

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