Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,631,245 C→T S171S (AGC→AGT slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,2450CT96.5% 176.3 / NA 57S171S (AGC→AGTsltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/24);  minor base G (0/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGC  >  NC_000913/4631003‑4631490
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
gcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAACCCGGGACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCCGAGGGAAGAAGCCTGGCAAGGGGCGAAAAAGCTGGGGCTAACCGGCAAAAGCCAGCCTAACGCCCGGGCCAAGTTATTTAGTGTCTGgcg                                                                                                                                                                                                                                                >  2:311163/1‑251 (MQ=255)
   aaCCCCACGCTTCCTACCGCTCGCGCGCTGCAATCGCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCTAAAGAGCTGTTGCTATCCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGgcgtg                                                                                                                                                                                                                                              <  1:213627/250‑1 (MQ=255)
   aaCCACACGCTGCCACCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAATACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTGCTACTATGCAAAATGGAACACCGGGCGGAGTGAAGAAGCCTGGCAAGGGGCGAAAGATCTGTGGCTAGCCGGCAAGCGCCAGCCTAACGCCTGTGGCAAGTTATTTAGTGGCTGgcgtgc                                                                                                                                                                                                                                             <  2:222914/251‑1 (MQ=255)
         ccGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAgg                                                                                                                                                                                                                                       <  1:454108/250‑1 (MQ=255)
                 tccCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACa                                                                                                                                                                                                                                 >  1:65185/1‑249 (MQ=255)
                   ccgctaccACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGACTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGAt                                                                                                                                                                                                                             <  1:121352/244‑1 (MQ=255)
                   ccGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCATCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGAt                                                                                                                                                                                                                             <  1:40548/251‑1 (MQ=255)
                    cGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATc                                                                                                                                                                                                                            <  2:364418/251‑1 (MQ=255)
                                               cAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGt                                                                                                                                                                                                  <  1:250964/250‑1 (MQ=255)
                                                  tGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGGAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGGAAACAAGATCCGCTGGCGTATTTAGAGCGGATCCGTCTgg                                                                                                                                                                                              >  1:324833/1‑251 (MQ=255)
                                                       gggCGCGGCGTGAAGACAGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACGACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATg                                                                                                                                                                                         <  2:222181/250‑1 (MQ=255)
                                                         gcgcGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGaa                                                                                                                                                                                       <  1:88515/251‑1 (MQ=255)
                                                                     gACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGGAGCGCAATGTAATTAATACTATGCAAAATGGACAACCGGGCAGAGAGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACCAGATATTTAGTGTCTGGCGTGCGTCAGGGAAACAAGATCCGCTGGGGTATTTAGAGAGGATCCGGCCGGCGATGGAAACGGGtgac                                                                                                                                                                             >  2:138277/1‑246 (MQ=255)
                                                                     gACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAac                                                                                                                                                                             <  1:138277/249‑1 (MQ=255)
                                                                           cGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCt                                                                                                                                                                     <  1:342286/251‑1 (MQ=255)
                                                                                  tGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTaa                                                                                                                                                                <  1:117585/249‑1 (MQ=255)
                                                                                      aGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGc                                                                                                                                                          <  1:311163/251‑1 (MQ=255)
                                                                                            tAGCCCGTAAAAGCCCTGAACTACCGAAGCGCAATGTAATTACTTCTATGCAAAATGGTACACCGGGCAGAGTGAAGATGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTAAGGGTCTGTCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggca                                                                                                                                                     <  2:206177/250‑1 (MQ=255)
                                                                                                ccGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggcagggca                                                                                                                                                <  1:173658/251‑1 (MQ=255)
                                                                                                         cccGGAACTACCGAAGCGCAATGTCATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGAGCGTCAGGTAAACAAGATCCGCGGGCGTATTTAGAGCGAAACCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGGAAAAGAGCTGGCAGGGCAGAtgcctac                                                                                                                                       >  2:345897/1‑249 (MQ=255)
                                                                                                            ggAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCg                                                                                                                                     >  1:193070/1‑250 (MQ=255)
                                                                                                                  aCCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTAcc                                                                                                                               <  1:473292/250‑1 (MQ=255)
                                                                                                                   ccGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTAcc                                                                                                                               >  1:210673/1‑249 (MQ=255)
                                                                                                                    cGAAGCGCTATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTg                                                                                                                                                                                                                       <  2:42785/160‑1 (MQ=255)
                                                                                                                    cGAAGCGCTATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTg                                                                                                                                                                                                                       >  1:42785/1‑160 (MQ=255)
                                                                                                                                 aaTTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGc                                                                                                               >  1:102747/1‑251 (MQ=255)
                                                                                                                                  aTTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGACTGTGACAAGATATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGAAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATTCCTGCCGATTACCAGACTATCGCCTCGGc                                                                                                               >  1:235606/1‑250 (MQ=255)
                                                                                                                                     actactATTCAAAATGGAAAACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGGTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGAAAAGTTATTTAGTGGCTGGCGTGCGTCAGGTAAACAAGATCCGCGGGCGTATGTAGAGCGTATCCGTCGGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTTCCAGACTATCGCCTCGGCAAt                                                                                                            <  2:38623/250‑1 (MQ=255)
                                                                                                                                      ctactaTGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTACAGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTACCACTGGACTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACCATCGCCTCGGCAATc                                                                                                           >  1:443669/1‑250 (MQ=255)
                                                                                                                                       tactaTGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCAt                                                                                                         >  1:277638/1‑251 (MQ=255)
                                                                                                                                          tatgcaacctGGAACACCGGGCAGAGTGAAGAAGCCTGGCCAGGGGCGAAAGAGATGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCGGATGCCTGCCGATTACCAGACTTTCGCCTCGGCAAtactttt                                                                                                      >  2:39711/10‑245 (MQ=255)
                                                                                                                                          taTGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTc                                                                                                      >  1:465578/1‑251 (MQ=255)
                                                                                                                                               aaaaTGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCCATCATTTCCCTgg                                                                                                 >  1:460553/2‑251 (MQ=255)
                                                                                                                                                     gAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGaac                                                                                            >  1:366460/1‑250 (MQ=255)
                                                                                                                                                       acacCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGaacaac                                                                                         <  2:65185/251‑1 (MQ=255)
                                                                                                                                                           cGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCt                                                                                      >  1:91169/1‑250 (MQ=255)
                                                                                                                                                              gCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAAt                                                                                   >  1:255012/1‑250 (MQ=255)
                                                                                                                                                               cAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAc                                                                                 <  2:193070/251‑1 (MQ=255)
                                                                                                                                                               cAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAc                                                                                 >  1:280920/1‑251 (MQ=255)
                                                                                                                                                                     gaagaaGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGGACt                                                                           >  2:164942/1‑251 (MQ=255)
                                                                                                                                                                         aagCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCa                                                                                                          >  1:429267/1‑216 (MQ=255)
                                                                                                                                                                         aagCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGGTATTTAGTGTCTGGCGGGCGGCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCa                                                                                                          <  2:429267/216‑1 (MQ=255)
                                                                                                                                                                                    ggggCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCGGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTAc                                                            <  2:210673/251‑1 (MQ=255)
                                                                                                                                                                                         gAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGCTGCCTGCCGATTACCAGACTATCGCCTCGGCAACCATTTCCCTGGAGAACAACCCTAATACGGTACTGACCTTCGAGCGaacaac                                                         >  2:315138/1‑249 (MQ=255)
                                                                                                                                                                                         gAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTg                                                      <  1:25933/251‑1 (MQ=255)
                                                                                                                                                                                                cggtgGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGAc                                               <  2:277638/249‑1 (MQ=255)
                                                                                                                                                                                                    ggCTAACCGGCAGGAGCCAGCCAAACGCCTGTGACAAGTTATTTAGGGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTTGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACGGTCGCTACCGAt                                           <  2:460553/251‑1 (MQ=255)
                                                                                                                                                                                                      cTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTCCTGACCTTCGCGCGTACAACTGGCGCGCCCGAttt                                         >  2:237990/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccagccTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAACTgg                           >  1:219526/1‑251 (MQ=255)
                                                                                                                                                                                                                         cTAACGCCTGTGACCAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGAAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCCCTGGCGAACAACCCTATTACGGCACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGgcggcg                      >  1:366279/1‑251 (MQ=255)
                                                                                                                                                                                                                         atAACGCCTGTGACCAGTTATTTAGTGTCAGGAGTGCGTCAGGTAACCAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACAATGGCCTGGTAACAGTGATGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGCCCGATTTTTCCCTTCAAATGgcggcg                      >  1:366284/2‑251 (MQ=255)
                                                                                                                                                                                                                          tAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGcggcgg                     >  1:266288/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaCGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGcggcggg                    >  2:219009/1‑250 (MQ=255)
                                                                                                                                                                                                                           aaCGCCTGTGACAAGTTATTTAGTGTCTGGCGGGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGgcggcg                      <  2:428999/249‑1 (MQ=255)
                                                                                                                                                                                                                             agCCTATGAAAAGTTATTTCGTGCCTGGGTTGCGTCAGGTAAACAAGATACGCTGGCGTTTTTAGAGCGTACCCGTCTGGCGATGCAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGACGCCTGACGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTCATAAGGTACTGACCTTCGCGCGTCAACACGGCGCGAACGGTTTTACCCGTCAAAAGGaggcggggg                  >  2:483482/2‑248 (MQ=255)
                                                                                                                                                                                                                               ccTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACCACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTgggg                >  2:349661/1‑249 (MQ=255)
                                                                                                                                                                                                                                             tttAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGgc  >  1:270997/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGC  >  NC_000913/4631003‑4631490

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