Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 157,585 C→T D50N (GAT→AAT)  folK ← 2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913157,5850CT98.4% 215.5 / NA 65D50N (GAT→AAT) folK2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (42/21);  total (43/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGT  >  NC_000913/157354‑157829
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
acaGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTggg                                                                                                                                                                                                                                   >  1:220900/1‑251 (MQ=255)
   gCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTg                                                                                                                                                                                                                                >  1:238693/1‑251 (MQ=255)
       aaaTCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCggg                                                                                                                                                                                                                            >  1:158698/1‑250 (MQ=255)
       aaaTCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGt                                                                                                                                                                                                                            >  1:328242/1‑251 (MQ=255)
       aaaTCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGt                                                                                                                                                                                                                            >  1:136709/1‑251 (MQ=255)
       aaaTCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGt                                                                                                                                                                                                                            >  2:435216/1‑251 (MQ=255)
       aaaTCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGt                                                                                                                                                                                                                            <  1:274940/251‑1 (MQ=255)
           ccACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCACACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGAGGTaaa                                                                                                                                                                                                                         >  2:67691/1‑250 (MQ=255)
            cACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAAc                                                                                                                                                                                                                       >  1:391336/1‑251 (MQ=255)
                       atatCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACgaagaa                                                                                                                                                                                                                 >  2:118548/1‑246 (MQ=255)
                       atatCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACgaagaa                                                                                                                                                                                                                 <  1:118548/246‑1 (MQ=255)
                               gTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGCCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAAtgt                                                                                                                                                                                                    <  1:586176/251‑1 (MQ=255)
                                 gCGGAATGGTCAGGAGCTCAGTATTTGTCACTTCATTACCAAACAGCATGCTGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTAAGCTTTGCGGCCGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGGGCAGGCGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGGGGGGTGCGGTAAACCGGAGAAACGGTAAGAAtgcgg                                                                                                                                                                                                  >  2:56836/1‑248 (MQ=255)
                                  cGGAACGGTCAGGCGTTCAGTCTTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGc                                                                                                                                                                                                 >  2:413413/1‑251 (MQ=255)
                                              gagTTCAGTATTTATCACTTCATTACCACACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCCGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTACGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGAGCGGGAAAACGACGACACGGTAAGAATGCGGCTTTCAGGGat                                                                                                                                                                                       >  2:529274/3‑249 (MQ=255)
                                               cGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATc                                                                                                                                                                                    >  1:170009/1‑251 (MQ=255)
                                                      tATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCATGCGTGGTCCCCAGCGTTCAGCTTTGAGGACGCGAACTTGCTGAAATTCAATACGATGTGTGTGATTGAGTCGCTCTTCAGGTGCAAGAGCGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGACCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTAAGGGAGATCGCCTAAt                                                                                                                                                                             >  1:351574/1‑251 (MQ=255)
                                                       aTTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAAt                                                                                                                                                                             <  1:292224/250‑1 (MQ=255)
                                                         ttATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCTTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGc                                                                                                                                                                           <  1:197479/250‑1 (MQ=255)
                                                                  cATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGgc                                                                                                                                                                  >  1:180959/1‑250 (MQ=255)
                                                                   aTTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGCAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTACTGCTTTCAGGgc                                                                                                                                                                  >  1:39258/1‑249 (MQ=255)
                                                                             aGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGCCCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGGAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTGCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTg                                                                                                                                                      <  2:153922/251‑1 (MQ=255)
                                                                                atgatgTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTaa                                                                                                                                                                              <  1:291507/224‑1 (MQ=255)
                                                                                atgatgTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTaa                                                                                                                                                                              >  2:291507/1‑224 (MQ=255)
                                                                                         agaTCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGATTTTTTCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGc                                                                                                                                                     <  1:464241/240‑1 (MQ=255)
                                                                                         agaTCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGGGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGc                                                                                                                                                     >  2:464241/1‑240 (MQ=255)
                                                                                            tCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGCAGCGACCTTGCTGCAATTCAATACGCTGTGTGGGATTGAGTAGCTCGTCAGGTGCAAGAGCGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCTGCAGCGGTGGGGTGCGGTAAAACGAAGACAAGGTAAGCATGTGGCTTTCCGGGTTATCGCCTAATGTTTTCAGGGCAGCTGTGGCCTGACCCAGCGGAGGGGc                                                                                                                                       >  2:588016/1‑251 (MQ=255)
                                                                                             ccAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCTAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCGCCAGCGGAGAGGcc                                                                                                                                      <  1:104139/251‑1 (MQ=255)
                                                                                             ccAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTAAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCTTTGACCTGCTCCCGAGGAGAGGcc                                                                                                                                      >  2:322534/1‑251 (MQ=255)
                                                                                                  gtgcgtgGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGAtt                                                                                                                                 >  1:283739/1‑251 (MQ=255)
                                                                                                                    tCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAatat                                                                                                                <  2:158698/250‑1 (MQ=255)
                                                                                                                      aGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACg                                                                                                             <  1:549038/251‑1 (MQ=255)
                                                                                                                      aGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGAAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGCGGCCAGATTGCTGCCTATGGCAATATACg                                                                                                             >  1:365860/1‑251 (MQ=255)
                                                                                                                       gCTTTGCGGACGCGACCTTGCTGCAATTCACTACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGc                                                                                                            >  2:142881/1‑251 (MQ=255)
                                                                                                                                      ccTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCAACGGCTGCGTTTAAGTACTCGGGTTGATTTTGCGGCCCAAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGCGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGACCA‑ATTGTTGCCTATGGCCATATAAGCCACTggaactccgtac                                                                                            >  1:67099/1‑240 (MQ=255)
                                                                                                                                         tgctgcAATTAAATACGCGGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGGTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTAc                                                                                            <  2:328242/249‑1 (MQ=255)
                                                                                                                                           ctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTc                                                                                        >  1:156082/1‑251 (MQ=255)
                                                                                                                                            tgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGAAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTcc                                                                                       >  1:318449/1‑250 (MQ=255)
                                                                                                                                              cAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATg                                                                                                  <  2:229788/238‑1 (MQ=255)
                                                                                                                                              cAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATg                                                                                                  >  1:229788/1‑238 (MQ=255)
                                                                                                                                                aTTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCacgac                                                                                   >  1:60106/1‑251 (MQ=255)
                                                                                                                                                   cAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACgtg                                                                                <  2:382676/251‑1 (MQ=255)
                                                                                                                                                          ctgtgtggcatcgagtcgccctttcggTGAAAGCGAGGTTCCCACAGCCACGGCTGCGTTTGAGACATCGGATTGATCTTGCGGCCCCAGCGGCGGGGTGCGGCGAAAAGACGAAAAGGTCAGAATGTGGCTTTCAGGGATATCGACTAATGCTTTCAGGGCAGCATTAACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCACTATACGCCACTGTCATGCCGTACCCTCACGCCGTGGcgcgcg                                                                         >  2:301086/26‑251 (MQ=255)
                                                                                                                                                           tgtgtgtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGATTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGt                                                                        >  1:21161/1‑251 (MQ=255)
                                                                                                                                                              gtgtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAAGCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGttt                                                                      <  1:19384/250‑1 (MQ=255)
                                                                                                                                                                gtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGCCGTGGTGCGCGTTTGCg                                                                   >  2:478316/1‑251 (MQ=255)
                                                                                                                                                                  gATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCAGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGt                                                                  <  1:155227/250‑1 (MQ=255)
                                                                                                                                                                    ttGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGGGGTGCGCGTTTGCGTggc                                                               >  2:538644/1‑250 (MQ=255)
                                                                                                                                                                       aGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCCGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAACCGACGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGACTATGGCAATATACGCCACTGTCCTGCGGTACCCTCCCGACGTGGTGCGCGTTTGCGTGgacg                                                             >  1:526355/1‑250 (MQ=255)
                                                                                                                                                                            ctctTCAGGTGCAAGAGAGGATTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGCGGGGTGCGGTACAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGCTATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCCATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGt                                                       >  2:292586/1‑251 (MQ=255)
                                                                                                                                                                                  aGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTacgacg                                                 >  1:205779/1‑251 (MQ=255)
                                                                                                                                                                                   ggTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCTGGTTGATGTTGCGGCTCCAGCGGTAGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCGTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTTCTGCCTATGGCAATATACGCCACTGTCATGCGGTAGCCTCACGGCGTGGTGCGCGTTTGCGTGGACGACGAGTacgacga                                                <  1:56836/251‑1 (MQ=255)
                                                                                                                                                                                   ggTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTacgacgc                                                >  2:125020/1‑250 (MQ=255)
                                                                                                                                                                                   ggTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTacgacg                                                 >  1:219237/1‑250 (MQ=255)
                                                                                                                                                                                    gTGCAAGAGAGGGTTCCAGCGCCACGGCTGCGTGTAAGGAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGTTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGTGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTacgacgac                                               <  1:588016/251‑1 (MQ=255)
                                                                                                                                                                                              ggTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACgg                                     >  1:254537/1‑251 (MQ=255)
                                                                                                                                                                                                         ccACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCa                          <  1:538644/251‑1 (MQ=255)
                                                                                                                                                                                                          cACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAg                         <  1:441393/251‑1 (MQ=255)
                                                                                                                                                                                                           aCGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGc                        >  1:113021/1‑251 (MQ=255)
                                                                                                                                                                                                               cTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGt                    >  2:144466/1‑251 (MQ=255)
                                                                                                                                                                                                                 gCGTTTAAGTAATCGGGTTGATTTTGCGGCCAAAGCGGTGGGGTGCGGTAAAACGCAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGGACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGCCGCGGGTCCGGATCTTCATCCAGCTCgatg                  >  2:60649/1‑248 (MQ=255)
                                                                                                                                                                                                                     ttAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGc               >  1:284178/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACTACGAGTACGACGACGCgg                                           <  1:155254/210‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCgg                                           >  2:155254/1‑210 (MQ=255)
                                                                                                                                                                                                                                 gTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCGCGTTGAGCATCCCTTTTTGGt  <  1:478316/251‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGT  >  NC_000913/157354‑157829

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