Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,394,248 C→T R61R (CGC→CGT)  nnr → NAD(P)HX epimerase/NAD(P)HX dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,394,2480CT98.5% 224.9 / NA 68R61R (CGC→CGT) nnrNAD(P)HX epimerase/NAD(P)HX dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (42/25);  minor base G (0/1);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CCCAGTTCCAGCCCCCACTGTTTAATTTTTTGCGCTAACTGATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGCCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGG  >  NC_000913/4394003‑4394494
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         
cccAGTTCCAGCCCCCACTGTTTAATTTTTTGCGCTAACTGATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCActg                                                                                                                                                                                                                                                      <  2:374442/251‑1 (MQ=255)
   aGTTCCAGCCCCCACTGTTTAATTTTTTGCGCTAACTGATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCActggct                                                                                                                                                                                                                                                   >  1:78964/1‑251 (MQ=255)
   aGTTCCAGCCCCCACTGTTTAATTTTTTGCGCTAACTGATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCActggct                                                                                                                                                                                                                                                   >  1:451835/1‑251 (MQ=255)
                          tttttGCGCTAACGGATTGAGATCGGGCGGCTTTGCCATGACGGCCCATACAATGAAGAAAAACCCCGTAAGTATACAACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGTCGGCAGATGTGCTGGGGTTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCTTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGtaat                                                                                                                                                                                                                             <  1:453646/250‑1 (MQ=255)
                                     aCTGATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGgcggcgg                                                                                                                                                                                                                  >  2:367602/1‑249 (MQ=255)
                                     aCTGATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGGCATGGTAATAACGgcggcg                                                                                                                                                                                                                   >  1:87541/1‑249 (MQ=255)
                                        gATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGTTGGc                                                                                                                                                                                                              >  1:444603/1‑251 (MQ=255)
                                        gATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCAAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGAttgg                                                                                                                                                                                                              >  2:427560/1‑248 (MQ=255)
                                           ggagaTCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTAc                                                                                                                                                                                                           <  2:184612/250‑1 (MQ=255)
                                              gaTCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGGTACggg                                                                                                                                                                                                        >  2:111336/1‑249 (MQ=255)
                                                      ggCTCTGACATGACGGACCATAAAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTcgcgcg                                                                                                                                                                                                <  2:451835/251‑1 (MQ=255)
                                                       gCTCTGACATGCCGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCGGTGCGGTCATGGTAATAACGgcggcgggggcctagggggcggcgcg                                                                                                                                                                                               >  2:175806/1‑231 (MQ=255)
                                                         tctGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGGGCTGTGCGGTCATGGGAATAACGgcggcggaggcgccgtgggcgggcgact                                                                                                                                                                                             >  2:261631/1‑229 (MQ=255)
                                                          ctGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACCCTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGCCATGGTAATAACGGCGGCGATGGCTACGTGGGCGAGCGACTg                                                                                                                                                                                            >  1:23062/1‑251 (MQ=255)
                                                               aTGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCa                                                                                                                                                                                        >  1:73949/1‑250 (MQ=255)
                                                               aTGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACAGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGAGGGCTACGTGGTCGCGCGACTGGCCa                                                                                                                                                                                        >  1:412814/1‑250 (MQ=255)
                                                                tGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGCGTCGTTCGGCGCATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACTGCGGCGAGTGCGACGGGGGCGCGCGACTGGCCac                                                                                                                                                                                       >  2:442677/1‑249 (MQ=255)
                                                                  aCGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGGTGGCTACGGGGTCGCGCGACCGGACAAAGc                                                                                                                                                                                    >  2:4226/1‑251 (MQ=255)
                                                                          tACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTCCGTGGTCGCGCGACCGGCCAAAGAGGTCGGc                                                                                                                                                                             >  2:265148/1‑250 (MQ=255)
                                                                           aCAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGGTGGCTACGTGGGCGCGCGACTGGCCAAAGCGGTCGGCA‑t                                                                                                                                                                          >  2:453372/1‑251 (MQ=255)
                                                                                gaagaaAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAgg                                                                                                                                                                     <  2:13643/251‑1 (MQ=255)
                                                                                   gaaAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGGCGCGCGACTGGCCAAAGCGGTCGGCC‑TTGAGG‑Tcc                                                                                                                                                                 >  1:214097/1‑250 (MQ=255)
                                                                                         ccccGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACggtgt                                                                                                                                                           >  1:447646/1‑251 (MQ=255)
                                                                                          cccGTAAGTATGCGTCACGCCGTCTGGTACGGCGAGGATAACCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCACGCACTCTAGTATCTGATGCTTCGAGCTGGCGAGGCCGCATTTCAGGTGTGTCGGTCGGAGTATCCTGACGTCCGGCACTGGCTGGTTCTGTGCGGTCATGGTAATAACGGCGGCTATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCTGCA‑TTGAGG‑ACACgttgtt                                                                                                                                                          <  2:250631/251‑1 (MQ=255)
                                                                                          cccGTAAGTATACCACACACCGTCGGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACgttgtt                                                                                                                                                          <  1:207291/251‑1 (MQ=255)
                                                                                           ccGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCTCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCTCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACTttgttg                                                                                                                                                         <  1:134415/251‑1 (MQ=255)
                                                                                              gaaGTATACCACACACCGGCTGGTAGGCCGACGAGTTCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCGGGCGAGTCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGcc                                                                                                                                                      <  2:214097/250‑1 (MQ=255)
                                                                                               aaGTATACCACACACCGTCTGGTACGCAGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCGATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACTTTGTTGGcc                                                                                                                                                      <  2:121833/250‑1 (MQ=255)
                                                                                                    taCCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGGCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAgg                                                                                                                                                  >  1:336856/1‑249 (MQ=255)
                                                                                                    taCCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGGAATAACGGCGGCGATGGCCACGTGGTGGCGCGCCTTGTCAAAGCCGTCCGCA‑TTGAGG‑TTACGTTTTTGGCCCAgg                                                                                                                                                  >  2:375340/1‑249 (MQ=255)
                                                                                                       cacacaCCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGGCCGGCCAAAACGGTCGGCA‑TTGAGG‑TCACGGTGGTGGCCCAGGAGAGc                                                                                                                                             >  2:13175/1‑251 (MQ=255)
                                                                                                       cacacaCCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGGCGCGCGACCGGCCAAAGCGGTCGGCA‑TAGAGG‑TCACGTTGTTGGCCCAgggaagg                                                                                                                                             >  1:221777/1‑246 (MQ=255)
                                                                                                            cccGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGaaaa                                                                                                                                        <  2:119651/250‑4 (MQ=255)
                                                                                                                      aCGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGGCGCGCGACCGGGCAAAGCGGGCGGCA‑TTGAGG‑TCACGGTGTTGGCCCAGGAGAGCGACAAACCGTTGCCg                                                                                                                              >  1:159289/1‑251 (MQ=255)
                                                                                                                          cgacgaTATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTCTGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACCGGGTCGCGCGACTGGCCCAAGCGGTCGGCA‑TTTAGG‑TCACGTTGTTGGCCCCAGCGAGCGACAAACCGTTTCCGGaa                                                                                                                           >  1:277357/1‑250 (MQ=255)
                                                                                                                             cgaTATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAgg                                                                                                                       <  2:297766/251‑1 (MQ=255)
                                                                                                                               atatCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGGTGGCTACGGGGGCGCGCGACTGGCCCAAGCGGGCGGGA‑TTTGGG‑TCCCGTTGTTGGCCCAGGAGAGCGCAAAACCGTTGCCGGAAGAgggc                                                                                                                     >  2:78583/1‑249 (MQ=255)
                                                                                                                               atatCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTATATGAGATGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACAGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGAGATGGCTACGTGGTCGCCCGACTGGCCAAAAAGGTCGGCA‑TTTAGG‑TCAAGGTGTTGGGGCAAGAAAGAGACAAAACGTTGCCGGGAAAGGcc                                                                                                                     >  1:404660/1‑251 (MQ=255)
                                                                                                                                     gccgcGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGGCGCATTCCAGGTGTGTCGTTCGGCGTATCCTAACGCCCGTCGCTGGCTGGTGCTGTGCTGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAATAGGCCGCGCTg                                                                                                               <  1:264814/251‑1 (MQ=255)
                                                                                                                                           gagaACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCCCGTTGTTGGCCCAGGAGGGCGACAAACCGGTGCCGGAAGAGGCCGCGCTGGCAcgc                                                                                                         >  1:335792/1‑251 (MQ=255)
                                                                                                                                                   gAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGa                                                                                                                            <  1:348405/224‑1 (MQ=255)
                                                                                                                                                   gAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATg                                                                                                 <  2:73949/251‑1 (MQ=255)
                                                                                                                                                   gAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATg                                                                                                 <  1:151367/251‑1 (MQ=255)
                                                                                                                                                   gAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGGAATAACGGCGGCGATGGCTACGGGGGCGCGCGACTGGCCAAAGCGGTCGGGA‑TTGAGG‑GCACGTGGTTGGGCCCGGGGAGGGACAAAACGGTGCCGGgagatcggaagagggtcgggtagggga                                                                                                  >  2:348405/1‑227 (MQ=255)
                                                                                                                                                         gcAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGGTACGTGGGCGCGCGACTGGCCAAAGCGGGCGGCA‑TTGAGG‑GCACGTTGTTGGCCCAGGAGAGGGACAAACCGTTGCCGGAAGAGGGCGCGGCGGCCCGGGAAGCATGGGTaa                                                                                            >  2:21209/1‑250 (MQ=255)
                                                                                                                                                            gATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGGTGGCTACGTGGGCGCGCGCCTGGCCAAAGCGGGCGGCC‑TTGAGG‑GCACGGTGGTGGCCCAGGAGAGCGCCAAACCGTTGCCGGAAGGGGCCGCGCTTGCCAGCGAAGCATGGTTAAAc                                                                                          >  2:359699/1‑249 (MQ=255)
                                                                                                                                                              tgtgCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCgg                                                                                     >  1:416086/1‑251 (MQ=255)
                                                                                                                                                               gtgCTGGGGCTCCCACTCTAGGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCTTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCgg                                                                                     <  1:316132/250‑1 (MQ=255)
                                                                                                                                                                 gCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTAGCGCGACCGGCCAAAGCGGTCGGCA‑TTGGGGGTC‑CGTTGTTGGCCCAGGAGATCGACCAACCGATGCCGGCAGAGGCCGCGATGGCCCGGTAAGCATTGTTAAAC‑GGGGGTg                                                                                  >  2:137713/1‑251 (MQ=255)
                                                                                                                                                                   tGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑GCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGGAGCATGGTTAAAC‑GCggggggg                                                                                >  1:125009/1‑247 (MQ=255)
                                                                                                                                                                    ggggCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGc                                                                                <  1:427560/250‑1 (MQ=255)
                                                                                                                                                                     gggCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGGGGTCGCGCGACTGGCCAAAGCGGGCGGCC‑TTGAGG‑TCACGTTGTTGGCCCAGGAGGACGGCAAACCGGTGCCGGGAGAGGCCCCGCCGGCACCCGCACCATGGGTTAAC‑GCggggggggg                                                                              >  2:170631/1‑245 (MQ=255)
                                                                                                                                                                           acacTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGGTGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATccc                                                                        >  1:47990/1‑250 (MQ=255)
                                                                                                                                                                                    gAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGTTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCTACCGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTgg                                                                                                                                                        <  2:289234/163‑1 (MQ=255)
                                                                                                                                                                                     aGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCCGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAat                                                              >  1:85424/1‑251 (MQ=255)
                                                                                                                                                                                      gCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACAGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAata                                                             <  2:167660/251‑1 (MQ=255)
                                                                                                                                                                                        tGATGCTTCGCGCTAGCGAGGCCGCATTCCAGGTGGGGCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGGGGTGGCGGGACTGGCCAAAGCGGGCGGGAGTTG‑GG‑TGAAGTTGGTGGCCCAGGAGAGGGAAAAAACGGTTCCGGAAAAAGCCGCGCTGGGACGGGAAGAATGTGTAAACAGAGGG‑GGCGGGGTCCATGCTTCGGataat                                                           >  2:333579/1‑249 (MQ=255)
                                                                                                                                                                                          aTGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGCATATTGt                                                         >  1:33336/1‑251 (MQ=255)
                                                                                                                                                                                           tGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGGCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGGGGCGAGATCCATGATTCGAATATTGTc                                                        >  1:276667/1‑251 (MQ=255)
                                                                                                                                                                                                cgctcTGGCGCGGCCTCATTCCAGGTGTGTCGTCCGGCGTATCCTGAAGCCCGTCCCTGGCTGTTGCTGTGCGGTCATGGTAATATCGGCGGCGATGGCTACGTGGTCGAGCGACTGGCCAAAGCGGTCGGCC‑TTGAGG‑TCCCGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAATAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGcc                                                   <  1:389409/247‑1 (MQ=255)
                                                                                                                                                                                                        cGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑GCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGGAGAGGCCGCGCTGGGACGCGGAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGGATATTGTCTGGCCCGAATCgg                                           >  1:101066/1‑251 (MQ=255)
                                                                                                                                                                                                                          gtgtCGTTCGGCGTCTCCTGACGCCCGTCTCTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATg                         <  1:442677/251‑1 (MQ=255)
                                                                                                                                                                                                                              cggtcgtcggaTCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCAGGGTAATAACGGCGGCGATGGCTACGTGGGTGCGCGACTGGCCAAAGCGGTCGGCA‑TGGAGG‑TCACGTTGTTGGCCCAGGAGAGAGACAAACCGTTGCCGGAAGAGTCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCgct                     <  2:122309/241‑1 (MQ=255)
                                                                                                                                                                                                                              cGTNCGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCgct                     <  2:416086/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGGGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTc                 >  1:112864/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gCGTATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCgg               <  1:144279/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tATCCTGACGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCGGGCACGCGAAGCATGGTTAAAC‑GGGGGGGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGGAGATCTGATTGTTGATGCGCTGCTCgggac            >  1:310163/1‑248 (MQ=255)
                                                                                                                                                                                                                                              aCGCCCGTCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTg     <  2:336856/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cccgtcACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCCCGGGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTTTCTGGCCCGAAACGGGAGATCTGATTTTTGATGTGCTGCTCGGTACCGGTTTGCgg  >  2:136100/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cccgtcACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCgg  <  1:251242/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         
CCCAGTTCCAGCCCCCACTGTTTAATTTTTTGCGCTAACTGATTGAGATCGAGGGGCTCTGACATGACGGACCATACAATGAAGAAAAACCCCGTAAGTATACCACACACCGTCTGGTACGCCGACGATATCCGCCGCGGAGAACGCGAGGCGGCAGATGTGCTGGGGCTCACACTCTATGAGCTGATGCTTCGCGCTGGCGAGGCCGCATTCCAGGTGTGTCGTTCGGCGTATCCTGACGCCCGCCACTGGCTGGTGCTGTGCGGTCATGGTAATAACGGCGGCGATGGCTACGTGGTCGCGCGACTGGCCAAAGCGGTCGGCA‑TTGAGG‑TCACGTTGTTGGCCCAGGAGAGCGACAAACCGTTGCCGGAAGAGGCCGCGCTGGCACGCGAAGCATGGTTAAAC‑GCGGGTGGCGAGATCCATGCTTCGAATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGG  >  NC_000913/4394003‑4394494

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