Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,595,919 C→T 6.0% 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,9190CT 6.0% 217.8 / 2.0 66D273N (GAC→AAC) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base C (35/27);  new base T (2/2);  total (37/29)
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.69e-01

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

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