Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 399,746 T→C 6.0% F52L (TTC→CTC)  yaiZ → DUF2754 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913399,7460TC 6.0% 215.8 / 2.6 66F52L (TTC→CTC) yaiZDUF2754 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base C (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TTCGTTGTGAGTCACAAAATGGCAACCTTAAGTTAAATTTATTAACAATTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACA  >  NC_000913/399468‑400024
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                      
ttCGTTGTGAGTCACAAAATGGCAACCTTAAGTTAAATTTATTAACAATTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCtggctggc                                                                                                                                                                                                                                                                              <  2:375728/289‑1 (MQ=255)
  cGTTGTGAGTCACAAAATGGCAACCTTAAGTTAAATTTATTAACAATTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGGTTTGAAGCAACAGGTTGGCAAACCTATGCCGTGGGTCTGGTGCCGTGGGAGATTAGTTACTGGGTGgggg                                                                                                                                                                                                                                                                            >  2:7096/1‑289 (MQ=255)
  cGTTGTGAGTCACAAAATGGCAACCTTAAGTTAAATTTATTAACAATTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCgg                                                                                                                                                                                                                                                                            >  2:350878/1‑289 (MQ=255)
         aGTCACAAAATGGCAACCTTAAGTTAAATTTATTAACAATTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTCATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGat                                                                                                                                                                                                                                                                     >  2:391142/1‑289 (MQ=255)
             acaAAATGGCAACCTTAAGTTAAATTTATTAACAATTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGattatt                                                                                                                                                                                                                                                                 <  1:227283/289‑1 (MQ=255)
                       aaCCTTAAGTTAAATTTATTAACAATTACGAAAATGTTGTCTACAACACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTcgcg                                                                                                                                                                                                                                                       >  1:38435/1‑289 (MQ=255)
                                   aaTTTATTAACAATTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGCTTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGCAATTg                                                                                                                                                                                                                                           >  1:158679/1‑289 (MQ=255)
                                               aTTACGCAAAGGTTGTCTACACCACAGAGTGTTTTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTTCAATGATACTGCCTGTCAATATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTTGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTTGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAAt                                                                                                                                                                                                                               <  1:201468/289‑1 (MQ=255)
                                               aTTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGAGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAAt                                                                                                                                                                                                                               <  1:6473/289‑1 (MQ=255)
                                                                   gtcacaGACGGTTGTTGCAGAATATGCAAGGATGTTGTTGTTCGTTAACTGAGGTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGATGGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGTTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCAGGCGTTa                                                                                                                                                                                                           <  1:12554/287‑1 (MQ=255)
                                                                   accacAGACGGTTGTTGTAGAATATGCAAGGATGTTGTTTTTCTTTAACGGAGCTGCCATTAATCTGCCTGTAAAAATCCGCCGGGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTa                                                                                                                                                                                                           <  2:38435/289‑1 (MQ=255)
                                                                   accacAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTa                                                                                                                                                                                                           <  2:193754/289‑1 (MQ=255)
                                                                   accacAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTCTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGa                                                                                                                                                                                                                           >  1:246900/1‑273 (MQ=255)
                                                                   accacAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTCTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGa                                                                                                                                                                                                                           <  2:246900/273‑1 (MQ=255)
                                                                     cacaGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTaaa                                                                                                                                                                                                         <  1:350878/289‑1 (MQ=255)
                                                                        aGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGCACACGCCCAATAAGAATATTATCGGGCGTTAAAata                                                                                                                                                                                                      >  2:356757/1‑289 (MQ=255)
                                                                               tgttgCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACAtt                                                                                                                                                                                               >  1:200183/1‑289 (MQ=255)
                                                                                      gAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGttt                                                                                                                                                                                        >  1:173513/1‑289 (MQ=255)
                                                                                        atatGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTc                                                                                                                                                                                      >  1:180565/1‑289 (MQ=255)
                                                                                           tGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGATTTCAAt                                                                                                                                                                                   <  1:10537/289‑1 (MQ=255)
                                                                                                     ttgttTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCggg                                                                                                                                                                                                                                                                           <  1:370366/191‑1 (MQ=255)
                                                                                                     ttgttTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCggg                                                                                                                                                                                                                                                                           >  2:370366/1‑191 (MQ=255)
                                                                                                       ggttttCGTGAACGGAGCGGACAGGAATCTGCCGGTAAAAATCCGCCGTGACGGGAACTACTATGCGTTCGTCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGTGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTGAAAATATTACATTGTGGTTTGAAATGCGTTATCCGCa                                                                                                                                                                       <  2:158679/287‑1 (MQ=255)
                                                                                                            tCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTTTGGTTTTCAATGCCTTATCCGACGCGTg                                                                                                                                                                  >  2:217670/1‑289 (MQ=255)
                                                                                                             cGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTg                                                                                                                                                                  >  2:115144/1‑288 (MQ=255)
                                                                                                             cGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGa                                                                                                                                                                 <  1:261300/289‑1 (MQ=255)
                                                                                                               ttAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGAc                                                                                                                                                                >  1:317714/1‑288 (MQ=255)
                                                                                                                                aaTCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGa                                                                                                                                                                                                                           >  2:83184/1‑212 (MQ=255)
                                                                                                                                aaTCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGa                                                                                                                                                                                                                           <  1:83184/212‑1 (MQ=255)
                                                                                                                                  tCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTa                                                                                                                                                                                                           >  1:355809/1‑226 (MQ=255)
                                                                                                                                  tCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTa                                                                                                                                                                                                           <  2:355809/226‑1 (MQ=255)
                                                                                                                                     gCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATTGCTTATCCGCAACGTTAAGCTCCAGCGCCCGTTCAATCaaa                                                                                                                                         >  2:4010/1‑289 (MQ=255)
                                                                                                                                        tGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACgt                                                                                                                                      <  2:173513/289‑1 (MQ=255)
                                                                                                                                          tAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACgtgt                                                                                                                                    >  2:127235/1‑289 (MQ=255)
                                                                                                                                                tCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGCTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTCCAATCCAACGTTTGATCAg                                                                                                                              >  2:87722/1‑289 (MQ=255)
                                                                                                                                                  cgccgTGACTGGCACTACTATGCGTTCTCCATTGGCCTTATATTCATTCTTAAGGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTAATCAGAt                                                                                                                            <  1:4010/289‑1 (MQ=255)
                                                                                                                                                  cgccgTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGAt                                                                                                                            <  1:115144/289‑1 (MQ=255)
                                                                                                                                                        gACTGGCACTACTTTGCGTTTGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgtg                                                                                                                       <  2:95287/288‑1 (MQ=255)
                                                                                                                                                        gACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgtg                                                                                                                       <  2:317714/288‑1 (MQ=255)
                                                                                                                                                               actactATGCGTTCGCCATTGGCCTTGTATTCATTCTTAATGGCGTGGTGGGGTGACTGGGATTTGAAGCAAAAGGGTGGCAGGCCTATGCCGTCGGTCTGGGGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGCGCGTTAAAATATTACATTGTGGTTTGCAATGCGTTATCCGAAGCGTGACGCTCCAGCGCCAGTGCAATCAAACGTGTGATCAGATCGTTGTAACCCa                                                                                                                <  1:356757/288‑1 (MQ=255)
                                                                                                                                                                           tCGCCATTGGCATTATATTCATTCTTAATGGCGTGGTGGGGTTCCTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCGGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTCTCGGGCGTTAAAATATTACATTGTGGTTTTCTATGCGTTATCCGCAGCTTGACGCTCCAGCTCAAGGTCAATCACACGTGCTATCAGATCGGTGTAAGCCAGACCGCTGGCtt                                                                                                    >  1:99380/1‑288 (MQ=255)
                                                                                                                                                                             gCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTAAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGcc                                                                                                 >  1:165406/1‑289 (MQ=255)
                                                                                                                                                                              ccATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGATGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCca                                                                                                <  2:181038/289‑1 (MQ=255)
                                                                                                                                                                                     ccgtataTTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTTGCAGACCTATGCCGTCGGTCTGGTGACGTGGCTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCtt                                                                                          <  1:7096/285‑1 (MQ=255)
                                                                                                                                                                                     ccTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCtt                                                                                          <  2:200183/288‑1 (MQ=255)
                                                                                                                                                                                         atatTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACAtt                                                                                                                                                                                               <  2:23754/183‑1 (MQ=255)
                                                                                                                                                                                         atatTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACAtt                                                                                                                                                                                               >  1:23754/1‑183 (MQ=255)
                                                                                                                                                                                         atatTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGat                                                                                     <  2:165406/289‑1 (MQ=255)
                                                                                                                                                                                          tatTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGGGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGata                                                                                    <  1:217670/289‑1 (MQ=255)
                                                                                                                                                                                               attcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCca                                                                                                <  2:91004/272‑1 (MQ=255)
                                                                                                                                                                                               attcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCca                                                                                                >  1:91004/1‑272 (MQ=255)
                                                                                                                                                                                                    tAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTGGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGata                                                                          >  2:41942/1‑289 (MQ=255)
                                                                                                                                                                                                           gtggtgGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGATATCCGCCGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTTTAACCCAGACCGCTGGCTTGCCCCAGCTTCGGATACATACTTATCTTAGtca                                                                   >  2:177637/1‑287 (MQ=255)
                                                                                                                                                                                                               tggGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGc                                                                <  1:177637/288‑1 (MQ=255)
                                                                                                                                                                                                                     tACTGGGATTTGAAGCAAAAGGTTGTCAGATGTATGCCGTCGGTCTGTTGACTTGGGGGATTAGTTGCTGGCGGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGTGCGTTAAAATATTACATGTTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCcaggca                                                          <  1:391142/288‑1 (MQ=255)
                                                                                                                                                                                                                     tACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCcaggca                                                          <  2:180565/288‑1 (MQ=255)
                                                                                                                                                                                                                        tGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGgt                                                      >  1:157489/1‑289 (MQ=255)
                                                                                                                                                                                                                                      aaaGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGa                                        >  1:169419/1‑289 (MQ=255)
                                                                                                                                                                                                                                                         ccGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTg                      <  2:169419/288‑1 (MQ=255)
                                                                                                                                                                                                                                                            tCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                                                                                                        >  2:97205/1‑187 (MQ=255)
                                                                                                                                                                                                                                                            tCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                                                                                                        <  1:97205/187‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gTGACGTGGGTGATTAGTCTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACgt                                                                                                                                      <  2:606/165‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gTGACGTGGGTGATTAGTCTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACgt                                                                                                                                      >  1:606/1‑165 (MQ=255)
                                                                                                                                                                                                                                                                       aCGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACg        <  1:407365/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTaca  <  1:127235/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTaca  <  2:157489/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                      
TTCGTTGTGAGTCACAAAATGGCAACCTTAAGTTAAATTTATTAACAATTACGAAAATGTTGTCTACACCACAGACGGTTGTTGCAGAATATGCAAGGATGTTGTTTTTCGTTAACGGAGCTGCCATGAATCTGCCTGTAAAAATCCGCCGTGACTGGCACTACTATGCGTTCGCCATTGGCCTTATATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACA  >  NC_000913/399468‑400024

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