Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,150,316 G→A E219K (GAG→AAG)  yghA → NADP(+)‑dependent aldehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,150,3160GA94.6% 221.2 / NA 74E219K (GAG→AAG) yghANADP(+)‑dependent aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (41/29);  minor base C (4/0);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCC  >  NC_000913/3150083‑3150555
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
ccATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTg                                                                                                                                                                                                                                >  1:302181/1‑251 (MQ=255)
 cATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTg                                                                                                                                                                                                                                >  1:177776/1‑250 (MQ=255)
   ttCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAat                                                                                                                                                                                                                             >  1:363731/1‑251 (MQ=255)
     ccGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGTCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAata                                                                                                                                                                                                                            <  1:258460/250‑1 (MQ=255)
      cGGATATTGCAGACCTCCCCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGAATTCGGGTGAATAtt                                                                                                                                                                                                                          >  2:295127/1‑251 (MQ=255)
      cGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATAtt                                                                                                                                                                                                                          <  1:565550/251‑1 (MQ=255)
           atTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTcgc                                                                                                                                                                                                                     <  2:177776/251‑1 (MQ=255)
              gccgACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTTGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGACGCAAGGAATTCGGGCGAATATTGTCGCGcc                                                                                                                                                                                                                  >  2:449375/4‑251 (MQ=255)
                      cACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGccc                                                                                                                                                                                                            <  1:33743/249‑1 (MQ=255)
                          aGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTg                                                                                                                                                                                                      >  2:454184/1‑251 (MQ=255)
                           gCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGATCTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTgg                                                                                                                                                                                                     >  1:479587/1‑251 (MQ=255)
                            cGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGa                                                                                                                                                                                                    <  1:428396/251‑1 (MQ=255)
                            cGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGa                                                                                                                                                                                                    <  1:485302/251‑1 (MQ=255)
                             gAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGAc                                                                                                                                                                                                   >  1:299971/1‑251 (MQ=255)
                             gAACAGTTTCAAAAGACCTTTGCCATTAACGTTTGCGCGCTGTGCTTGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTTCAAGTCTCATCTCCACTCCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGGCTATGCGGCTACGAAGGCGGCGAGTCTGAACTACAGCCGTGGCTTGGCAGAACAGGGCGCGAAGAAAGGGATTCGGGTGAATATGTTCGCGCCAGGCCCGATCTGGAc                                                                                                                                                                                                   <  2:515716/251‑1 (MQ=255)
                              aaCAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGAc                                                                                                                                                                                                   >  2:578231/1‑250 (MQ=255)
                              aaCAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGATCTGTCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGAGCAAGTATCATCACCAATTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACGATGCGGCTACGCAGACGGCGATTCTGAACTACAGCCGTGTCTTGGCAAAACAGGTCGCGAGGAAAGATATTCGGGTGCATATTGTCGCGCCCGTCCCGACCTggcac                                                                                                                                                                                                  >  2:308861/1‑248 (MQ=255)
                                cAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGc                                                                                                                                                                                                >  1:479694/1‑251 (MQ=255)
                                 aGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCa                                                                                                                                                                                               <  1:369579/251‑1 (MQ=255)
                                    ttCAAAAGACCTTTGCCATTACCGTTTTCGCGCTGTTCTGGCTAACCCAGGAAACGATCCCCCTGCTACCGAATGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGAATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACtc                                                                                                                                                                                            >  2:526780/1‑250 (MQ=255)
                                                ttGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAAGTCGCGCAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAACTTTCcggg                                                                                                                                                                                >  2:131840/1‑250 (MQ=255)
                                                  gcaattcaagTTTTCGCGCTGTTCCGGCTAACCCAGGAAGCAATCCCCCTGCTACCGAATGGAGCAGATATCTTCACCACGTCGTCACTCCAGGCATACCAGCCAAGTCCGCATTTACTGTACTATCCGGCTACGAAGGCGGCGATTCTGCACTACAGCCGTGGCTTGGCAAAACATGACGCGCAGACAGGAATGCGGGCGAATCTTGTCGCGGCATGCCCGCTCTGCACAGCCCGGCAAATTTCcggcgg                                                                                                                                                                              >  2:436787/10‑251 (MQ=255)
                                                   ccATTCACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGCTCCCCCTGCTACCGAAAGGCGCAAGTATCATCCCCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCTTTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAACGGTATTCGGGTGAATATTGTCGAGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcggcgg                                                                                                                                                                              >  2:241179/1‑250 (MQ=255)
                                                   ccATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGAATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcggcggg                                                                                                                                                                             >  2:559946/1‑250 (MQ=255)
                                                        aaCGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGAGTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAAc                                                                                                                                                                        <  1:296058/251‑1 (MQ=255)
                                                          cGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGACGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGc                                                                                                                                                                      >  1:499397/1‑251 (MQ=255)
                                                             tttCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAACGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAgg                                                                                                                                                                   >  1:481772/1‑251 (MQ=255)
                                                                 gcgcTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCACGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATa                                                                                                                                                                >  1:25740/1‑250 (MQ=255)
                                                                           ggCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCg                                                                                                                                                        <  1:400415/248‑1 (MQ=255)
                                                                           ggCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGACAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCAGATCTGGACCGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATccagc                                                                                                                                                       >  2:136718/1‑246 (MQ=255)
                                                                                   cAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCCCTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGAGGCTTGGCAAAACAGGTCGCGAAGAAAAGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGGTTGGTc                                                                                                                                              >  1:563512/1‑250 (MQ=255)
                                                                                    aGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTc                                                                                                                                              >  2:58402/1‑249 (MQ=255)
                                                                                                       aCCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCCGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCAGATGAAACGTg                                                                                                                         >  2:329274/1‑251 (MQ=255)
                                                                                                         cGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCg                                                                                                                       <  1:407280/251‑1 (MQ=255)
                                                                                                          gAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCgg                                                                                                                      >  1:316148/1‑251 (MQ=255)
                                                                                                             aGGTGCAAGTCTCAGCACCACTTCCTCAATCCAGGCATACCAGCCAAGTCCGCATTTAATGGACTATGCGGCTACGAAGGCGGCGATTCTAAACTACAGCCGTGGCTTGGCAAAACAGGCCGCGAAGAAAGGTATTCGGGAGAATATTGTCGCGCCAGGCCCGATCTGGACACCACTGCAAATTTCCGGCGGACACACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGCACCGTGCGGGGc                                                                                                                   >  2:224873/1‑251 (MQ=255)
                                                                                                                 gCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAAcc                                                                                                               <  1:116572/251‑1 (MQ=255)
                                                                                                                     gTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACcggcg                                                                                                           >  1:133232/1‑251 (MQ=255)
                                                                                                                     gTATCATCACCACTTCGTCAATCCAGGCATAACAGACAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGCGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCACGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACcggcg                                                                                                           >  2:201410/1‑251 (MQ=255)
                                                                                                                           tcaCCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGCCCCGCTCTGGACAGCACTGCAAATTTCCGGCGGACAACCGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTg                                                                                                     >  1:268401/1‑251 (MQ=255)
                                                                                                                             accacTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGc                                                                                                   <  1:429083/251‑1 (MQ=255)
                                                                                                                                   ccGTCAATCCAGGGAGACCAGCCAAGTCCGCATTTACTGGGCTATGCGTCTACGAAGGCGGCGATTCTGAACTACAGCCGGGGCGTGGCAAATCAGGTCGCGAAGAAAGGTAGTCGGGTGAATATTGTCGCGCCAGGCCCGTTCGGGTCAGCGCTGCGAATTTGCGGCGGACAGATGCAGGATGAGAGCCCGCGGGTGGGTCAGCAAACGCGGATGAAACGTGCggggttatcggcatagcggggccctg                                                                                              <  2:268401/249‑23 (MQ=255)
                                                                                                                                       caatcagggcatacacgacAAGTACGCCTTTACTGGACTATACGGCTACGAAGGCGGCCATTCTGAACTCCAGCCGTGGCTTGGCAAAACAGGTCGCGCAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCAGCAGATTGGTCAGCAAACGCCGATGAATCGTGCGGGGCAACCGGCGGAACTGGCCCCTGTttat                                                                                         >  2:250381/19‑251 (MQ=255)
                                                                                                                                        aacccAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATg                                                                                        <  2:302181/248‑1 (MQ=255)
                                                                                                                                         aTCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGt                                                                                       <  1:153040/251‑1 (MQ=255)
                                                                                                                                               gCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGCAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGCATATGTTTATCt                                                                                 >  2:204618/1‑251 (MQ=255)
                                                                                                                                               gCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCt                                                                                 <  1:718/251‑1 (MQ=255)
                                                                                                                                                cATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGTTCAATATTGTCGCGCCAGGCCCGATCt                                                                                                                                                                                                       <  2:493226/132‑1 (MQ=255)
                                                                                                                                                cATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCt                                                                                                                                                                                                       >  1:493226/1‑132 (MQ=255)
                                                                                                                                                 atgccAGCCAAGTCCGCATTTACTGGCCTCTGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCttg                                                                               <  2:446570/248‑3 (MQ=255)
                                                                                                                                                      aGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGCATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGa                                                                                                                                  >  2:197483/1‑195 (MQ=255)
                                                                                                                                                      aGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGa                                                                                                                                  <  1:197483/195‑1 (MQ=255)
                                                                                                                                                       gCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTc                                                                         >  1:227236/1‑251 (MQ=255)
                                                                                                                                                         cAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAg                                                                       <  2:481772/251‑1 (MQ=255)
                                                                                                                                                         cAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCACGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTc                                                                         >  1:74464/1‑249 (MQ=255)
                                                                                                                                                                 cATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCACACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCg                                                                 >  1:111996/1‑249 (MQ=255)
                                                                                                                                                                       cTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGt                                                         >  1:281441/1‑251 (MQ=255)
                                                                                                                                                                         ggACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATAATTTTTTCTGGCAAGTCAGGCGTCGAGGTACgccc                                                       >  2:309115/1‑248 (MQ=255)
                                                                                                                                                                                  ggCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACGGGCCCCTGGATATGTTTATCTGGCAAGGCCGGAGGCGAGCTACGGCACCGcgga                                                >  2:489674/1‑246 (MQ=255)
                                                                                                                                                                                   gCTACGAAGACGGCGCTTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGt                                             >  1:451192/1‑251 (MQ=255)
                                                                                                                                                                                         aaGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGATTCGAGCTACGTCACCGCAGAAGTGCACgg                                       <  1:203174/251‑1 (MQ=255)
                                                                                                                                                                                           ggcggcGATTCTGAACTACAGCCGTGGCTTGTCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCg                                     <  2:438805/251‑1 (MQ=255)
                                                                                                                                                                                            gcggcgATACTGGCCTACAGCCGTGGCTTGGAAAAACAGGTCGCGCAGAAAGGTATTCGGGTGCATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGGTTGGTCAGCAAACGCCGATGAAACGTGCGGAGCAACCGGCGGAACTGGCCCCTGTATATGTTTAGCCGGCAAGGCAGGATACGAGCTACGCCACCGCCGAAGAGCACGGccg                                    >  2:333823/1‑249 (MQ=255)
                                                                                                                                                                                               gcgATTCTGAACTACAGACGTGGCTTGGCAAACCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCCCGGCGtgtg                                 <  2:479694/251‑1 (MQ=255)
                                                                                                                                                                                                  aTTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTAGCGCCCGGCCCGATCTGGACAGAACTGCAAATTTCCGCCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTAAATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTgcgg                              >  1:570155/1‑251 (MQ=255)
                                                                                                                                                                                                           tACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGACGTGCACGGCGTGTGCGGCGGCGAGc                      >  2:183504/1‑250 (MQ=255)
                                                                                                                                                                                                             cAGCCGTGGCTTGGCAAACCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCTGACAAACGCAGTATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCATGAGTCGAGCTACGTCACCACAGAAGTGCACGGCGTGTGCGGCGGCGAGCAtt                   <  1:171250/251‑1 (MQ=255)
                                                                                                                                                                                                              aGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCAttt                  <  2:33888/251‑1 (MQ=255)
                                                                                                                                                                                                                   tGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGtt             >  1:114612/1‑251 (MQ=255)
                                                                                                                                                                                                                    ggCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGtt             <  1:526780/250‑1 (MQ=255)
                                                                                                                                                                                                                       ttGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGc                                           >  2:39858/1‑217 (MQ=255)
                                                                                                                                                                                                                       ttGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGc                                           <  1:39858/217‑1 (MQ=255)
                                                                                                                                                                                                                         ggCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTaaaaa        >  1:474491/1‑250 (MQ=255)
                                                                                                                                                                                                                              aaCAGGTCGCGAAGAAAGGTCATCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACCGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTGGGTCAGCAAACGCCGATGAAACGTGCGGGGCACCCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCAGTTAGGTTAAAAAAATGcc  <  2:65830/251‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
CCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCC  >  NC_000913/3150083‑3150555

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