Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,719,676 C→T 9.4% E16K (GAA→AAA)  pdxH ← pyridoxamine 5'‑phosphate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077791,719,6760CT 9.4% 242.2 / 3.6 81E16K (GAA→AAA) pdxHpyridoxamine 5'‑phosphate oxidase
Reads supporting (aligned to +/- strand):  ref base C (11/61);  new base T (6/2);  total (18/63)
Fisher's exact test for biased strand distribution p-value = 8.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TTTCGATTTGATGTGCTTTACGGCTGCCGAGGTTGGTGTAAAACACCATGCCTTTTTCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCC  >  NC_007779/1719439‑1719911
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
tttCGATTTGATGTGCTTTACGGCTGCCGAGGTTGGTGTAAAACACCATGCCTTTTTCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGGTCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCa                                                                                                                                                                                                                                               <  1:521019/239‑1 (MQ=255)
tttCGATTTGATGTGCTTTACGGCTGCCGAGGTTGGTGTAAAACACCATGCCTTTTTCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCa                                                                                                                                                                                                                                               <  1:570747/239‑1 (MQ=255)
                       cTGCCGAGGTTGGTGTAAAAAACCAGTCTTTTTTCGTCGAAATGTTGGAGAAAAACGATGCGCTGAAAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGGGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTtgct                                                                                                                                                                                                                        <  2:858760/238‑1 (MQ=255)
                         gCCGAGGTTGGTGTAAAACACCATGCCTTTTTCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAAAAGGGGGGGGGGGGAACCCGCCCTTGCGGGGTTACGGCGCAGAAGCGCCGTTttctgcc                                                                                                                                                                                                                     >  1:550649/1‑238 (MQ=255)
                                    tgtAAAACCCCATGCCTTTTTCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCCCCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATtcgt                                                                                                                                                                                                              <  1:1029874/234‑1 (MQ=255)
                                      tAAAACACCATGCCTTTTTCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATtcgt                                                                                                                                                                                                              <  1:489226/232‑1 (MQ=255)
                                           ccccttgtgtttttCGTGGTAATGTTTGAGAAAAACGTGGCGCTGTTAAGGATAACTATGTTCATCCACGGTAGCACCCCCCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACa                                                                                                                                                                                                  <  2:1637105/230‑1 (MQ=255)
                                              cATGCCTTTTTCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGAACGCGGCGGCGGAAACCCCCCTTTGTGTATTTACGGCGCAGAAGGGCGCGTTGTTGCCAATCTTCGGTTTTAGACACgg                                                                                                                                                                                              >  1:1637105/1‑239 (MQ=255)
                                                      tttCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGAACGCGGGGGCGGAACCCCCCCTTTGTGGATTTACGGCGCAGAAGCGCGGTTTGCTTCAAATCGTCGTTTTCATACATGGTTTTCttt                                                                                                                                                                                      >  1:96024/1‑239 (MQ=255)
                                                       ttCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTa                                                                                                                                                                                     <  1:1054866/239‑1 (MQ=255)
                                                       ttCGTCGTAATGTGTGATTAAAACGATGCGCGTATAAGTCTGTCTATTTTCCGCTCCGGTATCTACGACCATCGCGTTTGGGTCCGCCAGTTTGGGTTCACAAGCCTGAGAGAGCCAACTTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTa                                                                                                                                                                                     <  2:1774543/239‑1 (MQ=255)
                                                            cgtAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGa                                                                                                                                                                                 <  1:1541836/238‑1 (MQ=255)
                                                                      gAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCCCCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTggg                                                                                                                                                                      <  1:459549/239‑1 (MQ=255)
                                                                       ggTAAAGCGGTGCGTGGAAAGGGTTGACCAGTTTCATCCACGGAGGCGACCCCATTGGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTggg                                                                                                                                                                      <  2:1407759/237‑1 (MQ=255)
                                                                                                  cgaTGTTCTGCCACGGTAGCGCCCACAGTCGGGGTGGGGCCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGTCGCTATTGTGCGCCGCCCCTGGaaa                                                                                                                                          <  2:1294917/237‑1 (MQ=255)
                                                                                                  ccATGTTCATCCCCGTTAGCGACCTCCATCGCGGTAGGGTCTGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTTACGCTATTGTGCGCCGCCCCTGGaaa                                                                                                                                          <  1:1871303/239‑1 (MQ=255)
                                                                                                  ccATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGATCGCGGCGGGGGAACCCCCCCTTTTGGGATTTCCGGCGCCGATGGGCGGGTTTCTTCAATTTCTCGTTTTCAGACATTGGTTTTTTTTCGGATTGTTCGTGGGGGACGCTATTTTGCGGCGCCCCttggaa                                                                                                                                          >  1:1556998/1‑234 (MQ=255)
                                                                                                      gTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGATTGCGGGGGGGGAAACCCCCCTTTGTGGTTTTACGGGGCCGAAGGGCGGGTTGTTGGAATTTGTCGTTTTTAGa                                                                                                                                                                                                    >  1:499044/1‑178 (MQ=255)
                                                                                                      gTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATtcgt                                                                                                                                                                                                              <  1:701761/168‑1 (MQ=255)
                                                                                                         cATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCa                                                                                                                                   <  1:1770250/239‑1 (MQ=255)
                                                                                                           tCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGAACGCGGGGGCGGGACCCGCCCTTTGTGGGTTTTCGTCGCAGGAGCGCGGATTGCTTGAAATCCTTCGTTTTAGACATGGGTTTTTTTTAGTGTTTGTTAGGGGGGGGTGTTGTTTGGTGCCGCCCCGCGAAAAATCTCa                                                                                                                                   >  2:2119593/1‑238 (MQ=255)
                                                                                                           tCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAAc                                                                                                                                 <  1:1308709/239‑1 (MQ=255)
                                                                                                           tCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGAGCGGGGGCCGGGAAACCCCCTTTGGGGTGTTTACGGCGCAGAATGGCGCGTGTCCGTAATTCCGTCCTCTACAACACGGGTTTCTATTCGGGTTGTTCGTGGGGGAGGCCCTTGTGCGCCGCCCCCGGAGAAATCTCacc                                                                                                                                 >  2:1415120/1‑237 (MQ=255)
                                                                                                             cACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGc                                                                                                                               <  1:489009/239‑1 (MQ=255)
                                                                                                                                  ggTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAg                                                                                                          <  1:1757278/239‑1 (MQ=255)
                                                                                                                                   gTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAg                                                                                                          <  1:2119593/238‑1 (MQ=255)
                                                                                                                                   gTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTAGTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAg                                                                                                          <  1:739961/238‑1 (MQ=255)
                                                                                                                                    tAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTggg                                                                                                                                                                      <  1:292316/177‑1 (MQ=255)
                                                                                                                                         tCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTtgc                                                                                                                                                                                                                          <  1:467232/123‑1 (MQ=255)
                                                                                                                                           cGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTGGGTGTAGTCACGGCGCAGATGCGCGATTTGCGGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGaaa                                                                                                                                          <  1:1610183/198‑1 (MQ=255)
                                                                                                                                            gCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGa                                                                                                <  1:1864834/239‑1 (MQ=255)
                                                                                                                                             ccAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAACGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACt                                                                                                               <  1:1162092/223‑1 (MQ=255)
                                                                                                                                                     gCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCa                                                                                                                                                                           <  1:1110022/155‑1 (MQ=255)
                                                                                                                                                      cTTCACGAGCCTGAGAGAGCCAATGTTCAGAAAGGGTTAATGGATGGGGGGGAAGGTCGGGGCGGCGAAACCCGCATTTGGTGGATTCACGGCGCAGATGCTCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTGTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCACCGCTGTGGATTTTGTAACTGACAGTTAGTCAAGACGATGCggg                                                                                       <  2:550649/238‑2 (MQ=255)
                                                                                                                                                      cTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTTACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGt                                                                                       <  1:1926410/238‑1 (MQ=255)
                                                                                                                                                      cTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGt                                                                                       <  1:208211/238‑1 (MQ=255)
                                                                                                                                                      cTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGt                                                                                       <  1:662413/238‑1 (MQ=255)
                                                                                                                                                      cTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGt                                                                                       <  1:1366635/238‑1 (MQ=255)
                                                                                                                                                      cTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGt                                                                                       <  1:1415120/238‑1 (MQ=255)
                                                                                                                                                      cTTCACAAACCTGAGAGAGGCTACGTTCAAAAAGTGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAGCCCGCCTTTGGTGTATTCACGGCGCAGATGCGAGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGt                                                                                       <  2:11707/238‑1 (MQ=255)
                                                                                                                                                       ttCACAAGCCTGAGGGAGCCAACGTTAAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTcg                                                                                     <  2:357876/239‑1 (MQ=255)
                                                                                                                                                       ttCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTcg                                                                                     <  1:988132/239‑1 (MQ=255)
                                                                                                                                                       ttCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTcg                                                                                     <  1:1367123/239‑1 (MQ=255)
                                                                                                                                                       ttCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTcg                                                                                     <  1:342056/239‑1 (MQ=255)
                                                                                                                                                       ttCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTcg                                                                                     <  1:107739/239‑1 (MQ=255)
                                                                                                                                                       ttCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTcg                                                                                     <  1:486903/239‑1 (MQ=255)
                                                                                                                                                              gCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGAc                                                                                                                                                                                                   <  1:664997/123‑1 (MQ=255)
                                                                                                                                                              gCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGa                                                                                                                                                                                 <  1:691473/140‑1 (MQ=255)
                                                                                                                                                              gCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAg                                                                                                          <  1:1893066/211‑1 (MQ=255)
                                                                                                                                                               ccTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTa                                                                                                                                                                                     <  1:837447/135‑1 (MQ=255)
                                                                                                                                                                   agagagCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTggg                                                                                                                                                                      <  1:416024/146‑1 (MQ=255)
                                                                                                                                                                    gagagCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACa                                                                         <  1:1126941/238‑1 (MQ=255)
                                                                                                                                                                        gCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTc                                                                      <  1:948406/237‑1 (MQ=255)
                                                                                                                                                                                  aaaaaGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTATCCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTa                                                                                                                  <  1:1688917/183‑1 (MQ=255)
                                                                                                                                                                                         gTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGaaa                                                                                                                                          <  1:423337/152‑1 (MQ=255)
                                                                                                                                                                                          ttAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGAcc                                                   <  1:956822/238‑1 (MQ=255)
                                                                                                                                                                                          ttAATGGACCGGGGGGAAGATGGCGGGGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGAcc                                                   <  2:1333030/238‑1 (MQ=255)
                                                                                                                                                                                                        gggAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGAt                                      <  1:1838569/238‑1 (MQ=255)
                                                                                                                                                                                                        gggAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGAt                                      <  1:838394/238‑1 (MQ=255)
                                                                                                                                                                                                          gAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTc                                    <  1:1360631/238‑1 (MQ=255)
                                                                                                                                                                                                            aGATCGGGTCGGTTTAAACCACCTTTGGTGTATTCACGTGGTGGGTGGGGGATTTGGTTGTATTGGTCGTTATCAGACATGGTGTTCGTTACTGTTTGTCAGTGGGTGACGCTGTTGTGCTCCGCCCCTGGAAAAATCTCCACGCTGTGGATTTTGTATCTGACAGTTATTCAAGAAGATGCGGTCGCGTTTATAGACAGTCGCATCATCGCCTTTCGACCAGAAATCATAGAGTCCg                                  <  2:2024356/238‑1 (MQ=255)
                                                                                                                                                                                                             agtgcggggcggcgTAGCCGGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGTTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGt                                 <  2:1234954/232‑1 (MQ=255)
                                                                                                                                                                                                              aTCGCGGCGGCGTAACCCGCCTTTTGTGTACTCCCGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTACCTGACAGTTATTCAAGACGATGCGTTCGCGTTTATAGACAGTCTCCTCATCGCCTTTCGACCAGAAAAAATAAATTCCGGCa                               >  2:754018/1‑239 (MQ=255)
                                                                                                                                                                                                                     gcggcgTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGCTTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTaa                         >  2:1578949/1‑239 (MQ=255)
                                                                                                                                                                                                                     gcggcgTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACt                                                                                                               >  1:1416140/1‑153 (MQ=255)
                                                                                                                                                                                                                     gcggcgTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACt                                                                                                               <  2:1416140/153‑1 (MQ=255)
                                                                                                                                                                                                                         cgTAACCCGTCGTTGGTGTATTCACGGCGCAGGTGCGCGATTTTCTGCAATTCGTCGTTATCAGACATGGGTTGCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAAcgc                      <  2:663605/238‑1 (MQ=255)
                                                                                                                                                                                                                         cgTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGGCAGTGTAAcgc                      >  2:179035/1‑238 (MQ=255)
                                                                                                                                                                                                                                  ccTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAACGCGCGCCagag             >  1:1999568/1‑238 (MQ=255)
                                                                                                                                                                                                                                   cTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAACGCGCGCCAGAgg            >  2:1234722/1‑238 (MQ=255)
                                                                                                                                                                                                                                   cTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAACGCGCGCCAGAgg            <  1:1957315/238‑1 (MQ=255)
                                                                                                                                                                                                                                      tGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAg                                                                                                          >  2:51018/1‑141 (MQ=255)
                                                                                                                                                                                                                                      tGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAg                                                                                                          <  1:51018/141‑1 (MQ=255)
                                                                                                                                                                                                                                      tGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTcg                                                                                     >  1:449610/1‑162 (MQ=255)
                                                                                                                                                                                                                                      tGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGAATTTGTAACTGACAGTTATTCAAGACGATGCGGTcg                                                                                     <  2:449610/162‑1 (MQ=255)
                                                                                                                                                                                                                                        gtgtATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGCGGATTTTGTAACTGACAGTTATTCACGACGATGCGGTCGCGCTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCAGTCAGTGTAACGCGCTCCAAAGGCTGaaa      >  2:2004813/1‑239 (MQ=255)
                                                                                                                                                                                                                                         tgtATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGt                                                                       <  1:99990/173‑1 (MQ=255)
                                                                                                                                                                                                                                         tgtATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGt                                                                       >  2:99990/1‑173 (MQ=255)
                                                                                                                                                                                                                                            tttCCCGGCGCTGATGCGTCATTTGGTGCAATTCGTCGTTATCAGACATTGTTTGCTTTTCGGATTGTCCGTTGGTGACGCTATTGTGCGCCGCCCCTGGAAAACTCTCAAAGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTGTCGACCAGAAAACATAGATTCCGTCAGTGTATCGCGCGCCAGAGGCTGAAATGcc  <  2:4448/238‑1 (MQ=255)
                                                                                                                                                                                                                                            aTTTACGGCGCAGATGCGGGATTTGCTGCAATTCGTCGTGTTCATGCAGTGTTTTCTTTACGGATTGTCATTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGcc  <  2:1409573/239‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGc   >  1:1369391/1‑237 (MQ=255)
                                                                                                                                                                                                                                               cACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGc   >  1:167937/1‑235 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
TTTCGATTTGATGTGCTTTACGGCTGCCGAGGTTGGTGTAAAACACCATGCCTTTTTCGTCGTAATGTTTGAGTAAAACGATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCGGCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACATAGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCC  >  NC_007779/1719439‑1719911

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