Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,595,921 C→T 6.2% D273N (GAC→AAC)  ypfJ ← putative neutral zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,595,9210CT 6.2% 210.2 / 2.1 64D273N (GAC→AAC) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base C (34/26);  new base T (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCGGAACTCGCGCCCAAGTAAAGTTTGTAGTGCCGAGGGAGAACAGTGGTTTTCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACTTTCGGCTCGATGC  >  NC_000913/2595638‑2596205
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                            
acGGAACTCGGGCCCAAGTAAAGTTTGTAGTGCCGAGGGGGAACAGTGGTTTTCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCGTTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAAc                                                                                                                                                                                                                                                                                         <  2:292561/288‑1 (MQ=255)
                       tttGTAGTGCCGAGGGAGAACAGTGGTTTTCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGt                                                                                                                                                                                                                                                                  >  2:481146/1‑289 (MQ=255)
                        ttGTAGTGCCGAGGGAGAACAGTGGTTTTCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCg                                                                                                                                                                                                                                                                   <  1:200333/287‑1 (MQ=255)
                         tGTAGTGCCGAGGGAGAACAGTGGTTTTCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTTGAAACCACGTTTAAACCAGCTGTAGCGtt                                                                                                                                                                                                                                                                 >  1:308385/1‑288 (MQ=255)
                                   aGGGAGAACAGTGGTTTTCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTACGGCACCACGCCACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGCTTTATTTGCGCTGGTACTGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGACCCGGAAGTTAAATGCTTTTATCAAAAGTATTGCATTGTGCCGGATCGCCGCAATCGAACCTACGTTTAAACCAGCTGTAGAGTTGCCTAGAAGt                                                                                                                                                                                                                                                       >  2:31698/1‑288 (MQ=255)
                                               ggTTTTCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGaaa                                                                                                                                                                                                                                          <  1:481146/289‑1 (MQ=255)
                                                   ttCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTga                                                                                                                                                                                                                                                            >  1:436865/1‑267 (MQ=255)
                                                   ttCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTga                                                                                                                                                                                                                                                            <  2:436865/267‑1 (MQ=255)
                                                        cATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAATTGTATGGt                                                                                                                                                                                                                                 >  1:416642/1‑289 (MQ=255)
                                                                  cggcgAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGc                                                                                                                                                                                                                        >  1:91911/1‑288 (MQ=255)
                                                                                    gTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTa                                                                                                                                                                                                     >  2:754/1‑289 (MQ=255)
                                                                                              aaGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTTGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCAt                                                                                                                                                                                            <  2:308385/288‑1 (MQ=255)
                                                                                               aGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGtactac                                                                                                                                                                                                                            >  1:429446/1‑255 (MQ=255)
                                                                                               aGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGtactac                                                                                                                                                                                                                            <  2:429446/255‑1 (MQ=255)
                                                                                                ggCATCACGCAACTTAAGTGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTCTCACTTCCGCCGTTAATGTGTGAAGCGCAGTCAGTTCAGCCTTCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCCTAGTTCCGGATCGCTGCTGTCAAAACCACGTTTAAACCAGCTGTAGAGTTGCTCAGAAGTGCCATGAGAGAAACTGTATGGTACTCCTCGCTACTGACTTTCCTGCTGTAATGGGTCACCg                                                                                                                                                                                          >  2:255190/1‑288 (MQ=255)
                                                                                                               aaGGGTAGGCTGAAAACACCAACCCTCTTCCCCGGTGAACACCAGGAAGCGGCGGATCCCGTCACGTTTCATTTGCGCTGTTAATGTGTGGAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTGACCAAAAGGATGGCATTGTGCCGGATCGGGGCTGTCGAAACGACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGGGCCATGAGTGAAACTGTCTGGTACTACTCGCGCCGGACTGTGCTGTTGTAAACGGTCATGGCCGATGGCCTGcgcc                                                                                                                                                                           <  2:416642/288‑1 (MQ=255)
                                                                                                                      tGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCa                                                                                                                                                                    >  2:285263/1‑288 (MQ=255)
                                                                                                                      tGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCa                                                                                                                                                                    >  2:355207/1‑288 (MQ=255)
                                                                                                                       gCTCAAAACACCAACCCTCTACCCCGCTCAACACCAGCAAGCGGCGGAGCCCGTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTGCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACAGTCATCGCCGATGGCCTGCGCCGCGTTCAgc                                                                                                                                                                  <  1:255190/289‑1 (MQ=255)
                                                                                                                                          ttCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTTGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGCAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTTACTTTGCTGTTTTAAACGGTCATCGCCGCAGGCCTGCGCCGCGGTCCGCGCCTCCTCCAgctacccgg                                                                                                                                               >  2:27797/1‑281 (MQ=255)
                                                                                                                                               cGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTgaga                                                                                                                                                                                                                                                          >  1:247302/1‑177 (MQ=255)
                                                                                                                                               cGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTgaga                                                                                                                                                                                                                                                          <  2:247302/177‑1 (MQ=255)
                                                                                                                                                catCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTCCATTGTGCCGGATCGCCGCTGTCGAAACCACCTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTTTAAACGGTCATCGCCGATGGCCTGCCCCCCGTTTCGCGCCTCTTCCCGATTCCCGGTTTcc                                                                                                                                          >  2:237169/3‑288 (MQ=255)
                                                                                                                                                           gCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTCCGCCGCGTTCCGCGCCTCTTACAGATCGCCGGTGTCCAGACTGCCtt                                                                                                                                >  2:496779/1‑287 (MQ=255)
                                                                                                                                                           gCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTCTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTACACGGTCATCGCCGATGGCCTGCGCCGCGTTCCGCGCCTCTTCCAGATCCCCGGTTTCCAGAACGCCtt                                                                                                                                >  1:178185/1‑287 (MQ=255)
                                                                                                                                                             aaGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTTGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTtg                                                                                                                               <  1:27797/286‑1 (MQ=255)
                                                                                                                                                             aaGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTtgctg                                                                                                                            <  1:496779/289‑1 (MQ=255)
                                                                                                                                                              aGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTtgctg                                                                                                                            >  2:183468/1‑288 (MQ=255)
                                                                                                                                                                      aTCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCAt                                                                                                                      <  1:296397/286‑1 (MQ=255)
                                                                                                                                                                          cTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGc                                                                                                               <  2:91911/289‑1 (MQ=255)
                                                                                                                                                                                gTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACAAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGGGCCTCTTCCAGATCGCCGGTTTCGCGCACCCCTTGCTGCCGCATACTACGcacccc                                                                                                          >  2:296397/1‑287 (MQ=255)
                                                                                                                                                                                               ttAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAg                                                                                          >  2:362408/1‑289 (MQ=255)
                                                                                                                                                                                                         aaGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGcccc                                                                                                                                                                                                                     <  2:364906/156‑1 (MQ=255)
                                                                                                                                                                                                         aaGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGcccc                                                                                                                                                                                                                     >  1:364906/1‑156 (MQ=255)
                                                                                                                                                                                                         aaGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGa                                                                                >  1:438828/1‑289 (MQ=255)
                                                                                                                                                                                                            cgcAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGt                                                                              <  1:306711/288‑1 (MQ=255)
                                                                                                                                                                                                               agtcagtTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTcc                                                                           <  1:754/288‑1 (MQ=255)
                                                                                                                                                                                                               agtcagtTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCa                                                                          <  2:377326/289‑1 (MQ=255)
                                                                                                                                                                                                                     ttCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTGTGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGc                                                                     >  1:451221/1‑288 (MQ=255)
                                                                                                                                                                                                                      tCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACg                                                                      <  1:183468/286‑1 (MQ=255)
                                                                                                                                                                                                                          ccATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTtgct                                                                                                                             <  2:240995/227‑1 (MQ=255)
                                                                                                                                                                                                                          ccATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTtgct                                                                                                                             >  1:240995/1‑227 (MQ=255)
                                                                                                                                                                                                                               ccTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTTCATACTATGCCCCCAGCCCCCGGCAAAACAGTCGGCCTGGAGTTCCATACGCCCAGCTAAGCg                                                          >  2:119243/1‑289 (MQ=255)
                                                                                                                                                                                                                                        gAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCTGCTCGCCGGTTTCCCGAACGCCTTGCTGCTGCATCCTATGCACCCAGACACCGGCCAAACAGTCGGCCTGGAGTTCCATTAGCACAGATCAGCGGttacctt                                                  >  2:221997/1‑283 (MQ=255)
                                                                                                                                                                                                                                        gAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTtg                                                                                                                               <  2:68339/211‑1 (MQ=255)
                                                                                                                                                                                                                                        gAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTtg                                                                                                                               >  1:68339/1‑211 (MQ=255)
                                                                                                                                                                                                                                        gAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACg                                                                      >  1:277678/1‑268 (MQ=255)
                                                                                                                                                                                                                                        gAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACg                                                                      <  2:277678/268‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGt                                   >  2:11260/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         gtgCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGa                  >  1:311784/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                          tgCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACt               >  1:43223/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                             cGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACttt             <  1:119243/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACTTTCgg          >  2:2849/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                 tCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACTTTCGGCt        <  1:31698/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   gccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACTTTCGGCTCg      <  1:285263/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCAt                                                                                                                                                                                            <  2:411567/106‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCAt                                                                                                                                                                                            >  1:411567/1‑106 (MQ=255)
                                                                                                                                                                                                                                                                                    ccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTg                                                                                  <  1:290651/212‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTg                                                                                  >  2:290651/1‑212 (MQ=255)
                                                                                                                                                                                                                                                                                    ccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGc                                          >  1:267998/1‑252 (MQ=255)
                                                                                                                                                                                                                                                                                    ccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGc                                          <  2:267998/252‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACTTTCGGCTCGa     >  1:228158/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                    ccgcTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACTTTCGGCTCGa     >  1:96590/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                        tGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACTTTCGGCTCGATGc  <  1:355207/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                            
CCGGAACTCGCGCCCAAGTAAAGTTTGTAGTGCCGAGGGAGAACAGTGGTTTTCAGCATCTGGCCGCGGCGAAATCCACAGCCAGTCGCCAGGTAAGGCATCACGCAACTTAAGAGTATGCTCAAAACACCAACCCTCTTCCCCGCTCAACACCAGCAAGCGGCGGATCCCTTCACGTTTCATTTGCGCTGTTAATGTGTGAAGCGCAGTCAGTTCAGCCATCCCTGCCCCGGAAGTTAAATGCTTTTACCAAAAGTATTGCATTGTGCCGGATCGCCGCTGTCGAAACCACGTTTAAACCAGCTGTAGCGTTGCTGAGAAGTGCCATGAGTGAAACTGTCTGGTACTACTCGCCCCTGACTTTGCTGTTGTAAACGGTCATCGCCGATGGCCTGCGCCGCGTTCAGCGCCTCTTCCAGATCGCCGGTTTCCAGAACGCCTTGCTGCTGCATACTATGCCCCCAGACACCGGCAAAACAGTCGGCCTGGAGTTCCATACGCACAGATAAGCGGTTTACTTCCGCCTGCGTCGCGTTTTGTTGCAGTTGACGAACTTTCGGCTCGATGC  >  NC_000913/2595638‑2596205

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: