Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 910,811 T→G 23.6% H80P (CAC→CCC)  poxB ← pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913910,8110TG23.6% 94.0 / 6.3 66H80P (CAC→CCC) poxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (22/27);  new base G (11/6);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTT  >  NC_000913/910581‑911046
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
cAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGTAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACTCGATTCAGGTGGCAATCGAACAGGCCGTCGATTAAGTGCAGGtg                                                                                                                                                                                                                                         >  1:828214/1‑237 (MQ=255)
cAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAAAAGGCCGTTGATTAAGTGAAGGt                                                                                                                                                                                                                                          >  1:1201703/1‑237 (MQ=255)
           ccgcGGTTAAGCCCCGCTTTGCCATTGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAACCCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCa                                                                                                                                                                                                                            <  1:305634/236‑1 (MQ=255)
           ccccGGTTAAGCACCGCTTTGGGAATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGCCTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCa                                                                                                                                                                                                                            <  1:570158/236‑1 (MQ=255)
           ccACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCa                                                                                                                                                                                                                            <  1:1122712/239‑1 (MQ=255)
                       acccttttGCGCATGGCAATCGCCAGTCTTTGTGGGATTTCCCCTGGGCTGGAAACCAGTTCGTAATAGGGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGGGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGgcg                                                                                                                                                                                                                 <  2:1178228/232‑1 (MQ=255)
                          gtttggccattGGGAACCGCCATTACTTGGGGGATCTGCCCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGAc                                                                                                                                                                                                            <  1:916631/229‑1 (MQ=255)
                           cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTTCCGGGGCCGCACGATCCGGCGCAGAc                                                                                                                                                                                                            >  1:195886/1‑238 (MQ=255)
                           cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCCCGATCCGGCCCAGAc                                                                                                                                                                                                            >  1:1025636/1‑238 (MQ=255)
                           cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGGGCAGAc                                                                                                                                                                                                            >  1:68069/1‑238 (MQ=255)
                           cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGAc                                                                                                                                                                                                            >  1:489825/1‑238 (MQ=255)
                           cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGAc                                                                                                                                                                                                            >  1:607300/1‑238 (MQ=255)
                           cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGGGCAGGTTGCCGGGGCCGCACGATCCGGGCCAGAc                                                                                                                                                                                                            >  1:578938/1‑238 (MQ=255)
                           cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATACAGGGCAGGTGGCCGGGGCCGCACGATCAGGCGCAAAc                                                                                                                                                                                                            >  1:626204/1‑238 (MQ=255)
                              tGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATCTGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGGGCAGGTTGCCGGGGCCGCACGTACCGGCGCAGACcccc                                                                                                                                                                                                        >  2:707770/1‑239 (MQ=255)
                                cgcATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGCTTGCGGCGGCAATCGCACAGGCCGTTGATTAATGGCAGGTTGCCGGGGCAGCACGATCCGGCGCAGACCGCCAg                                                                                                                                                                                                      >  1:83642/1‑239 (MQ=255)
                                  gttGGCAGTCGCGTGTATTGGTGGGTTCTGCCCGGGGCGGGAAACCAGCTAGAAATAGTGACTACATTCGCGGAATAGCTTTTGTGGGTGGGTTTCCTGGAAATAGCGGAGGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGtt                                                                                                                                                                                                    <  2:437758/237‑1 (MQ=255)
                                  cATGGCAATCGCTAGTCCTTGTGGGATTTGCTCCGGGCTGGAACCCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTGCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:372520/239‑1 (MQ=255)
                                   ctGGCAATCGCCAGTACTTGGGGGATCTGCCCCGGGCTGGAACCCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGGGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:175489/237‑1 (MQ=255)
                                   atggacaTCGCCAGTACTTGGGGGATCTGCTCGGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGGGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:1371652/232‑1 (MQ=255)
                                   agggCAATCGCCGGGTCTTGTGGGAACTGCCCCGGGCTGGAAACCAGCTCGCAATAGTGACTACAGTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:837906/236‑1 (MQ=255)
                                   aTGGCAATCGCCAGTACTTGTGGGATCTGCCCCGGGCTGGAACCCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:707770/238‑1 (MQ=255)
                                     ggCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCCCTGGAGGGAATATGAGCGGAAATCCCCCGTCCCGAACCTGGGTTGGGGGGACAACGGAACGGCGCGTGATGGAAATGCAGGTTGCCCGGGCCCGCCGACCCGGCGCAGAACGCCaatcggc                                                                                                                                                                                                 >  2:837906/1‑233 (MQ=255)
                                     ggCAATCGCCAGTACTTGTGGGATCGGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGACCGTGCTTGCGATGGCAATCGAAAAGGCCGTTGATTAATTCAAGGCTGCCGGGGCCGCACgcaccgccaccagaccccagtgac                                                                                                                                                                                                   >  2:207562/1‑215 (MQ=255)
                                       caaccccttgttcTTTTGGGTTTTGCTCGGGGTTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGTTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGTTTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCg                                                                                                                                                                                               <  2:489825/227‑1 (MQ=255)
                                          tCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAACCCAGCTGGCAATAGTGACTACATTCGCGGAATAGCTCTTGGGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTa                                                                                                                                                                                            <  1:541941/239‑1 (MQ=255)
                                          tCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAATGGCAGGTTGCCGGGGCCGCACCATCCgggccaaccccccaggtccccggta                                                                                                                                                                                            >  2:1031685/1‑216 (MQ=255)
                                                             ttGCCCCGGGGTGGAACCACGCTCGAAATGGGGATTGCATTCGCGAAATAGCTCTTGTGGGTTGGTTTCGGGGAAATGGCCGCTGCCAATGTCGCGGAAGGGAATAGGAGCGGCAATCGCCAGTCCCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCa                                                                                                                                                                         <  2:644671/237‑1 (MQ=255)
                                                              tGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCCGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCacc                                                                                                                                                                       >  1:550534/1‑237 (MQ=255)
                                                                cTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGCGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTCATTTTGGCTTCGGCGccgacg                                                                                                                                                                      >  1:925698/1‑234 (MQ=255)
                                                                cTCAGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACTGGCCGTTGATTAAGGGCAGGTTGCCGGGGCCGCACGATCCGGCGCCGCCAGCACGTATTCAGCTAGGTTGGGCTTCAGGGccgacg                                                                                                                                                                      >  2:41758/1‑234 (MQ=255)
                                                                   cGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGGTTAAGGGCAGGTTGCCGGGGCCGCACGACCCGGGGCAGCCCGCCCGTTCTCCCCGAAGtgggggttcgtggcccagggcaa                                                                                                                                                                   >  2:1092345/1‑217 (MQ=255)
                                                                   cGGGCTGGAAACCAGCTCGCAATAGCGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGGGCAGGTTGCCGGGGCCGCACGATCCGGGGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCaa                                                                                                                                                                  >  1:701241/1‑239 (MQ=255)
                                                                        gggaacccggctcgatatgtttaCTACATTGGGGGAATAGCTGTTGGGGGGGGTTGTCCTGAAAAGACCCGCTGCCAATTCCGGCGGGGGAAATTGGGGGGGCAACCGCCAGTCCCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGCCCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGgc                                                                                                                                                              <  2:1025636/216‑1 (MQ=255)
                                                                          gAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTTCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGAGCACAACCGCCAATTCTCCCCGAACTTTGGCTTCAACGCACAGGGCAAAGggcgc                                                                                                                                                           >  2:140565/1‑239 (MQ=255)
                                                                               cAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGGTTAAGGGCCGGTTGCCGGGGCCGCCCGATCCGGCGCCGACCTCCGGTTCTCCCGTAAGTTGTGTTCCAGCCCCACCGGCAAAAGCCCCCAcc                                                                                                                                                       >  2:222733/1‑237 (MQ=255)
                                                                                  cTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAggaccccccgtct                                                                                                                                                    >  1:240051/1‑227 (MQ=255)
                                                                                              aCTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCa                                                                                                                                                                                                                            <  1:948694/156‑1 (MQ=255)
                                                                                              aCTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCa                                                                                                                                                                                                                            >  2:948694/1‑156 (MQ=255)
                                                                                              aCTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGGGGCgggggg                                                                                                                                       >  1:73694/1‑239 (MQ=255)
                                                                                               cTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAGGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCCGTTCTCCGCTAAGTTGGGCTTCACCGCAAGGGGCAAAGGCCGCCACTCCTCCGggggggggggg                                                                                                                                      >  2:523538/1‑238 (MQ=255)
                                                                                                ttcctttgccggatcgctgtgggggggggggTTCCTGGAAAGTACCGCCTGCAAATTTCGCGGGGGGAAAAATGGACGGAAATGCCAAGTACGGGACGGTGATGGGGGGGGCAAGCAAACAGGCCGTTGATTAAGGGCAGGTTGCGGGGGCCGCACGATCCGGCGCAGACCGCCGGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGa                                                                                                                                      <  2:1112667/216‑1 (MQ=255)
                                                                                                tacatcgggggattggCTCTTGGGGGGGGTTTTCCGGGAAATGCCCGCTGCCAATTTCGGTGGAGGGAATTTGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGAc                                                                                                                                     <  1:867004/224‑1 (MQ=255)
                                                                                                   aTTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCACCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGGGGTgtgggct                                                                                                                                  >  1:1112517/1‑232 (MQ=255)
                                                                                                   aTTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGGGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCGGCGAAGTTGGGCTTCAGCGCCGGCGGCAAAAGCCCCCACTTCTTCGgggggggggggagc                                                                                                                                   >  2:242008/1‑234 (MQ=255)
                                                                                                     tCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGGGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCTGGGCggggggggacgc                                                                                                                                 >  2:663232/1‑232 (MQ=255)
                                                                                                                ctctTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACCGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGATCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGACGGGGTGACATCCCCTCCgtggtg                                                                                                                      >  1:486058/1‑238 (MQ=255)
                                                                                                                    tgggggtgggtTCCCGGGAAATTGCCGCTGCCATTTCCGTTGGAGGGAATATGGGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCAt                                                                                                                 <  1:66576/236‑1 (MQ=255)
                                                                                                                     gtgggtggGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCggggggggcttccccccatgggggccat                                                                                                                 >  1:1001820/1‑216 (MQ=255)
                                                                                                                     gtgggtggGCTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCCAAGGTGTGCTTCAGCGCAAACAGCAAAGACCTCCACTTCTTCGGGGCGGGTggggccgcaccactcagtgcccccg                                                                                                                >  2:455742/1‑216 (MQ=255)
                                                                                                                     gtggggggGTTTTCTGGAAATAGCCCGTGCCAATTGTGGGGGGGGGAAATTGGGGGGGAATTCGCACTGACGGGACGGTGGTTGGGGGGGTAATCGAACAGGCCGTTGATTAAGGCCCGGTTGCCGGGGCCGCAGGATCCGGCGAAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCAAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCAt                                                                                                                 <  2:1014587/239‑1 (MQ=255)
                                                                                                                          tgggtTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCCCGATCCGGCGCAGACCGCCAGTTATCCGCTAAGTTGGGCTTCAGCGCCAGCGGCAAAGGCCCCCACTTCTTCGTGGGGGGGGGGCATcacccccgaggggcccccgcggtt                                                                                                           >  2:414982/1‑216 (MQ=255)
                                                                                                                                tCCTGGAAATGGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACt                                                                                                     <  1:242008/239‑1 (MQ=255)
                                                                                                                                cccggaaaaTGGCGGCTGCCAATTTCGCTGGAGGGAATATGGGCGGCAATCGCCAGTACCGGACGGTGATTGCGGGGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCAt                                                                                                                 <  1:1368767/221‑1 (MQ=255)
                                                                                                                                 gcggtgaaattcccgtggCCATTTTCGGTGGGGGTAATAGGGGGGGCAACCGCCGGGACGGGAGGGTGTTGGCGGTGCCAATGGATCAGGCCGTTGATGAAGTCGGGTTGGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACt                                                                                                     <  2:607300/222‑1 (MQ=255)
                                                                                                                                  tccgtgaaagaccGGTGCCAATTTCCCTGGGGGGAAATTGGGGGGGAATCCGCACGGCCCGGAACTGGATGGGGTGGGCAATGGAACAGGCCGTTGATTAAGTGCAGGTGGCCGGGGCCGCCCGTTCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGTCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACt                                                                                                     <  2:550534/228‑1 (MQ=255)
                                                                                                                                              cggtGCAAATTTCGTTGGAGGGAATATGAGCGGCAATCGCCAGGACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGAc                                                                                          <  1:1007576/234‑1 (MQ=255)
                                                                                                                                                  cgcattttCGATGGGGGGAATAGGGGGGGAAATCGACGGGCCGGGAAGGGGATGGGGGTGCCAAGCGAAAGGGCCGTTGATTAAGGGCAGGTTGCGGGGCCCGCACGACCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTc                                                                                     <  2:291038/233‑1 (MQ=255)
                                                                                                                                                               gAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGACATCCACTCGATGGGGCCCCTGCGCTTAGAAATGTCCCTAAAACCGCTctcgaagaggcccc                                                                        >  1:291038/1‑225 (MQ=255)
                                                                                                                                                                  ggaaatggAGCGGCAATCGCAAGTACCGGAACGTGATGGGGGTGGCAATCGAAAAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGa                                                                     <  1:918274/231‑1 (MQ=255)
                                                                                                                                                                  gaaaTAGGGGCGCCACTCGCAAGTACCGGAACGGTATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGa                                                                     <  1:455742/236‑1 (MQ=255)
                                                                                                                                                                    aaTAGGAGCGGCAATCGCCATTACCGGACCGTGATGGCGGTGCAAATCGAGCAGGCCGTTGATTAAGTGCAGGTGTCCGGGCCCGCACGAGCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGGGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGa                                                                                                                           <  1:668625/182‑1 (MQ=255)
                                                                                                                                                                             ctgtaaatcgcaatgtcgggAACGTGATTGCGGTGGAAACCGACCAGGCAGTTGCTTAAGGGCAGGTTGCCGGGGCCGCACGATCCGGCCAAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGa                                                            <  1:1031685/220‑1 (MQ=255)
                                                                                                                                                                                             ccGGAATGTTGTTGGGGGGGCAATCGACCCGGCCGTTGATTAGGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGcc                                           <  1:523538/239‑1 (MQ=255)
                                                                                                                                                                                                    gttgttGCGGTGCCAATGAACCAGGCGTTGGATTGAGGGGCGGTCGCGGGGCCCGCACTGCCCGGCGCAGCCCGCCAGTCCTCCGCTAGGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGa                                    <  2:195886/235‑1 (MQ=255)
                                                                                                                                                                                                                                       aaGGGCAGGTTGCGGGGGCCGCACGACCCGGCCCAGCCCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCgtttgttt  <  1:663232/239‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
CAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTT  >  NC_000913/910581‑911046

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