Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,537,714 A→C 15.1% V351G (GTG→GGG)  nanM ← N‑acetylneuraminic acid mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,537,7140AC15.1% 182.1 / 2.5 66V351G (GTG→GGG) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (45/9);  new base C (2/10);  total (47/19)
Fisher's exact test for biased strand distribution p-value = 2.10e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAG  >  NC_000913/4537486‑4537946
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
gcgcATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTgc                                                                                                                                                                                                                                  <  2:1066249/238‑1 (MQ=255)
   ccttttCCCTCGCCGGTAAGGGCAGTACTTCCCCAAACCCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAAAGATGGGGCGATATTAGTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCACTGAAGCGGTAGTTTTGTACTGTGACTTTATTATCCTGCACAGTGATCAAAACTGAATCCGTCACCGCTTTgccggc                                                                                                                                                                                                                              <  2:688215/235‑3 (MQ=255)
                        gCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGATTTGTACTGTGACTTTATTATCCTTAACAGTGATCAAACCTGAATCCGTCACCGCTTTGCCGCCTTCAGTTTCACCGCCaataat                                                                                                                                                                                                         >  1:540755/1‑239 (MQ=255)
                             ccTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCt                                                                                                                                                                                                                                       <  2:125984/205‑1 (MQ=255)
                             ccTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCt                                                                                                                                                                                                                                       >  1:125984/1‑205 (MQ=255)
                             ccTTCGCCAGAGGCCATGGCATTAGACTGACCAGCAACGGTAAAAACATAGTAATAATCGGGCAATATTATTGGTTTCATATAAATTCCGCGTAATGGAAATGACAAATTGATAGAAGGGCCCCTTTACCGGGGCCTTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTAAACCCGTCCCCGTTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATa                                                                                                                                                                                                    <  1:1099312/239‑1 (MQ=255)
                              cTTCGCAAAACGCCATGGCATTAGACTGACCAGAAACGGTAAAAACATAGTAATAATCGGGCGATTTTATTGCGTTCATATAAATTTCGGGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTCCGGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTAATCAAAACTGAATCCGTCCCCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAg                                                                                                                                                                                                   <  1:352974/239‑1 (MQ=255)
                              cTTCCCCATACCCCATGGCATTAGACTGACCAGCAACGGTAAGAACATGGTAAAAACGGGGGATTATTATTGGGTTCATAAAAATTCCGCGTAATTGAAATGAAAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCTTTAGTTTTGTACTGTGATTTTATTATCCTTCACAGTGATAAAAACTGAATCCGTCCCCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAg                                                                                                                                                                                                   <  1:93527/239‑1 (MQ=255)
                                  gCCATACGCCATGGCTTTAGACTGACCAGCACCGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTATCAGTGATCAAAACTGAATCCGTCACGGCTTTCACGCCGGCAGTTTCACCCCAAAAAATCAAAAGACta                                                                                                                                                                                               >  2:1066687/1‑239 (MQ=255)
                                              ggCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGc                                                                                                                                                                                   >  1:543746/1‑239 (MQ=255)
                                               gCATGAGACTGACAAGCAACGGAAAAAAAATAGTAAAAATCGGGCAATATTTTGCGGTTCATAAAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGCCCCTGAAGCGTTAGTTTTGTCCTGTGACTTTATTATCCTTCACAGTGATCAAAACTAAATCCGTCCCCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAggg                                                                                                                                                                                    <  1:709018/237‑1 (MQ=255)
                                                          aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCCCCGCTTTGCCCCCTGCAGTTTCACCCCCAATAATCAATAGACTATTATTCCCGGGCAATGTTaac                                                                                                                                                                          >  2:709018/1‑234 (MQ=255)
                                                          aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTgccgcc                                                                                                                                                                                                                              <  2:768535/184‑1 (MQ=255)
                                                          aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTgccgcc                                                                                                                                                                                                                              >  1:768535/1‑184 (MQ=255)
                                                          aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGCCTATTATTCCAGGGCACTGATACTCCg                                                                                                                                                                       >  1:204705/1‑239 (MQ=255)
                                                          aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGCCTATTATTCCAGGGCAATGATACCCCg                                                                                                                                                                       >  1:706139/1‑239 (MQ=255)
                                                          aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGGGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAACCAATAAACAATTATTCCAGGGCAATGAAACCCCg                                                                                                                                                                       >  2:13898/1‑239 (MQ=255)
                                                                          aCATAATAAAAAGCGGGCAATTTTATTGCTTTCATATAAATTTCGCGTAATTAAAATGACAAATTGAAAGCAGGCCCCTTTTCCCGGGGCCCTGAAGCGTTAGTTTTGTCCTGTGACTTTATTATCCTTCACAGTGATAAAAACTGAATCCGTCACCGTTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCg                                                                                                                                                       <  1:696325/239‑1 (MQ=255)
                                                                              aGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATCGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATaa                                                                                                                                                   >  1:467579/1‑239 (MQ=255)
                                                                                   taatCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATACTTCACAGTGATCAAAACTGAATCCGCCACCGCTTTGCCGCCTGCAGTTTTACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACATTCCGATAAttaa                                                                                                                                               >  1:862287/1‑236 (MQ=255)
                                                                                        gggCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTCTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAAATGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGTGATAATTCACCCg                                                                                                                                           >  1:231804/1‑237 (MQ=255)
                                                                                          gCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATCCTCCGTAGGCCCGACCTTGCGATAATTAACCCGAttt                                                                                                                                       >  1:685037/1‑239 (MQ=255)
                                                                                          gCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGa                                                                                                                                                      >  1:362652/1‑224 (MQ=255)
                                                                                                          ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                                                                                                                                                                                            >  1:10409/1‑172 (MQ=255)
                                                                                                          ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGCTACTCCGTAGGCCCGACCTTCCGATAATTCACCCGATTTATCCCATTTCCCGt                                                                                                                         >  1:25779/1‑239 (MQ=255)
                                                                                                          ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGCTACCCCCTAGGCCCGACCTTGCGATAATTAACCCGATCTATACCATTTCCCGt                                                                                                                         >  2:79317/1‑239 (MQ=255)
                                                                                                          ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGCTTTATCCCATTTCCCGt                                                                                                                         >  1:313084/1‑239 (MQ=255)
                                                                                                          ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGt                                                                                                                         >  1:325168/1‑239 (MQ=255)
                                                                                                          ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCCCAGGCCCAACCTTGCGATAATTCACCCGATTTATCCCATTTCCGGt                                                                                                                         >  2:646809/1‑239 (MQ=255)
                                                                                                          ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATAAATAGACTATTATTCCAGGGCAATGATACCCCTTAGCCCGGACCTTGCGTAAATTACCCCGCTTTATCCCATTTccctt                                                                                                                         >  1:148773/1‑237 (MQ=255)
                                                                                                               ataAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACATTTCGGAAATTCACCCGCTTTATACCATTTCCCGTTAtac                                                                                                                    >  2:627456/1‑237 (MQ=255)
                                                                                                                   attgtgcgTAATTGGAATGTCAAATTGTATACATGTGCCCCTTTACGGGGCCCTGGACGTTTGATTTTTTGTCTGGGCTTTTATTACCTTACTGATTGAAAAAACCTTAACCCCACCCCCTTTTCCCGCTTGACGTTTCCCCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCaa                                                                                                                 <  2:25779/234‑1 (MQ=255)
                                                                                                                       cgcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATa                                                                                                                                                                            >  1:148626/1‑175 (MQ=255)
                                                                                                                       cgcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTAATCCCATTTCCCGTTATCCCAAAAAta                                                                                                            >  2:293286/1‑238 (MQ=255)
                                                                                                                       cgcgTAATCGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCACTGATCCTCCGTAGCCCCGACCTTGCGTTAATTACCCCGGTTTATCCCATTTCCCTTTAGTCCaaaccta                                                                                                            >  2:681226/1‑235 (MQ=255)
                                                                                                                       cgcgTAATGGAATGAAAAATTTGATAGCAGGCCCCCTTTCCCGGGCCCCGGAAGCGTTAGTTTTGTACGGTGATTTTATTACCCTTCACAGTGATCAAAACTGAATCCGTCCCCGCTTGGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATg                                                                                                            <  1:611345/239‑1 (MQ=255)
                                                                                                                           aaaTTGAAATGACAAAGTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGGTTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTTATAAAAACTGAATCCGTCACCGTTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAata                                                                                                        <  1:293286/238‑1 (MQ=255)
                                                                                                                               tttaaTGACTAATTTATAGGAGGGCCCCTTTCCCGGGCGCACTAAAGGTTTTTTTTTTCTCTGTCTCTTTATTTCCCTTCACAGTAACAAAAATTAACTCCGCCCCGATGTGCCCGCCTGCATTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCa                                                                                                     <  2:313084/236‑1 (MQ=255)
                                                                                                                                   aTGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTACTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTg                                                                                                  >  1:219086/1‑237 (MQ=255)
                                                                                                                                      aCAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTACCCCGATTTATCCCATTTCCCGTTATGCCACAGATGAATATCAGTCTtata                                                                                            >  1:603214/1‑239 (MQ=255)
                                                                                                                                       cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCg                                                                                                                                                                       >  1:624125/1‑164 (MQ=255)
                                                                                                                                                    cAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGa                                                                                                                                                                                                  >  1:267826/1‑126 (MQ=255)
                                                                                                                                                       ggCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCATTAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCTCCCGATTTATCCCATTTCCCGTTATGCCACAGCTTAATATCAGCGGTAAAAGGTTTTTACAGGcc                                                                             >  2:723067/1‑239 (MQ=255)
                                                                                                                                                       ggCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAgg                                                                               >  1:311387/1‑237 (MQ=255)
                                                                                                                                                       ggCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAAAATCAGTGCTAAATGATTTTTTAAgg                                                                               >  2:408768/1‑237 (MQ=255)
                                                                                                                                                              ttACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTTCCGCCTGCCGTTTCCCCGCCAATAATCAATAGACTCTTATTCCCGGGCAATGATACTCCGTAGGCCCGACCTTGCGTTCATTCACCCGGTTTATCCCATTTCCCCTTTTGCCAAAAGTGAAATTCAATGCCTATTTATTTTTTCCCGCCTTCAtcc                                                                      >  2:257013/1‑237 (MQ=255)
                                                                                                                                                                 cgggggCCTGGAAGCGTTAGTTTGTCACTGTGTCTTAATTACCCTCCCCGGTAACCAAACCTAAATCCGGCCCCGTTTTGCCGCCGGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCa                                                                   <  1:646809/237‑1 (MQ=255)
                                                                                                                                                                 ccGGGGCCCTGAAGCGTTTGTTTTGTACTTTGATTTTATTATCCTTCACAGTGATAAAAACTGAATCCGTCACGGCTTTGCCGCCGGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCa                                                                   <  1:257013/239‑1 (MQ=255)
                                                                                                                                                                 ccGGGCCCCTGAAGCGTTAGTTTTGTACTGTGAATTTATTATCCTCCACATTGATAAAAACTGAATCCGTCCCCGCTTTGCCGCCTGCAGTTCCCCCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCa                                                                   <  1:408768/239‑1 (MQ=255)
                                                                                                                                                                     ggCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGt                                                               >  1:344568/1‑239 (MQ=255)
                                                                                                                                                                       gCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAatat                                                                                                       >  1:783583/1‑199 (MQ=255)
                                                                                                                                                                        cccTGAAGCGTTAGTTTTGTACTGTGACTTTATTATACTTCAAAGTGACAAAAACTGAATCAATACCCGATTAGCCGACTGCACATACACCGACACAACTCCACAGACAAATATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTCTCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGt                                                               >  1:594693/1‑237 (MQ=255)
                                                                                                                                                                                  ttAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGCTAATTCACCCGATTTATCCCATTTCCCGTTAGGCCAAAGATGAATATCAGTCCTAATTGATTTTTTCAGTCCTTCATGCGCCTATTTCTTACCGt                                                      >  2:1060333/1‑236 (MQ=255)
                                                                                                                                                                                            tACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGtt                                                              >  1:1085912/1‑219 (MQ=255)
                                                                                                                                                                                                  tgACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTccc                                                                                                                           >  1:120077/1‑154 (MQ=255)
                                                                                                                                                                                                  tgACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGa                                                                                         >  1:1039625/1‑187 (MQ=255)
                                                                                                                                                                                                        tattatCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCt                               >  1:556505/1‑239 (MQ=255)
                                                                                                                                                                                                            atCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGa                           >  1:650548/1‑239 (MQ=255)
                                                                                                                                                                                                                ttCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGCTTTATCCCATTTCCCGTTATGCCAAAGACGAATATCAGTGCTATATGATTTTTTCCGGCCTTCATACGCATAGTTCTTACCGTTTTGGAAATTTTTTCGTGAACCTTGGAAAcc                       >  2:541591/1‑239 (MQ=255)
                                                                                                                                                                                                                  cacaGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATCCGCa                                                                   >  2:907484/1‑193 (MQ=255)
                                                                                                                                                                                                                  cacaGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTAATCCCATTTCCCGTTATGCCAAAGATGAATATCAGCGCTATATGATTTTTTCAGGCCTTCATGCGCATATTTCTTACCGTTCTGGTAATTTTCTCTGGAAACtttaacacga                      >  2:607220/1‑231 (MQ=255)
                                                                                                                                                                                                                  cacaGTGATCAAAACTAAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCa                                                                   <  1:907484/193‑1 (MQ=255)
                                                                                                                                                                                                                       ttacaaaaaCTGAATCCGTCACCGTTTTCGCCTCTGCATTTTCCCCGCCAATAATCAATAGACTATTATACCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCt                <  1:175038/234‑1 (MQ=255)
                                                                                                                                                                                                                        aaacAAAACGGAATCCGGCCCCGTTTTCCCCCCTGGATTTCAACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCt                <  1:1060333/235‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAAtata    >  1:139830/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTTCGATACTTCAACCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCCATATGATTTTTTCAGGCCTTCATGCCAATATTCATTGCCGTTCTGGAAATTTTCCTGTTGACCTTTGATTCCGGCCCCTCCCCCaaaacaaaa  >  2:1036328/1‑233 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
GCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAG  >  NC_000913/4537486‑4537946

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