Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,587,922 G→A G64G (GGC→GGT ugpQ ← glycerophosphodiester phosphodiesterase UgpQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,587,9220GA98.3% 209.1 / ‑4.1 60G64G (GGC→GGTugpQglycerophosphodiester phosphodiesterase UgpQ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/29);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

CAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCA  >  NC_000913/3587682‑3588168
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
cagcGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTgg                                                                                                                                                                                                                                              >  1:211030/1‑251 (MQ=255)
    ggcggcgTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACg                                                                                                                                                                                                                                          <  2:322556/251‑1 (MQ=255)
    ggcggcgTCATACCGGCCCAAAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGCGCCGGTGGCGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGCCTCACCTTTAAACATTTTGATGTACCAACTGCCCGAATCCACGCGCAGTAAATCCTGCCCGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTCGTACg                                                                                                                                                                                                                                          >  2:347199/1‑251 (MQ=255)
       ggcgTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTc                                                                                                                                                                                                                                       >  1:142972/1‑251 (MQ=255)
                       acaGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATgg                                                                                                                                                                                                                       <  1:311166/251‑1 (MQ=255)
                          gcccgcgTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGaac                                                                                                                                                                                                                    <  2:11115/248‑3 (MQ=255)
                          gCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGaac                                                                                                                                                                                                                    <  2:191767/251‑3 (MQ=255)
                              gcgTGCCGCCAGCGCCACCATTTTGCCCGTTAATTGCCCGGTGCCGGGGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCATCTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCACAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGa                                                                                                                                                                                                                <  2:182419/251‑1 (MQ=255)
                                gTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAg                                                                                                                                                                                                              <  2:335582/251‑1 (MQ=255)
                                    cgccAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGCTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGAtct                                                                                                                                                                                                          >  1:406943/1‑251 (MQ=255)
                                       cAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGc                                                                                                                                                                                                       >  1:97877/1‑251 (MQ=255)
                                       cAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGc                                                                                                                                                                                                       >  1:21249/1‑251 (MQ=255)
                                           gcCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGATGGGTTTGCTTTCGATATTCGCCCTCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCACGCATCCACGCGCAGTCAATCCTGCCAGTTAAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATc                                                                                                                                                                                                   >  2:394356/1‑251 (MQ=255)
                                                        ccGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTc                                                                                                                                                                                      <  1:3579/251‑1 (MQ=255)
                                                             aaTGGCCCGGTGCCGGTGGTTGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTc                                                                                                                                                                                 <  2:395671/251‑1 (MQ=255)
                                                              aTGGCCCGGCGCCGGTGGTGGGTTTGATTTCGATATTCGCCATAATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGTAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCCACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGATCAGGAAGATCTCGCCATCTTTCGCTAACTTCGCGTc                                                                                                                                                                                 >  2:58988/1‑250 (MQ=255)
                                                                              ggtggGTTTGATTTCGAGATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTCCCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTAt                                                                                                                                                                <  1:130845/251‑1 (MQ=255)
                                                                                  ggTTTGATTTCCATATTCTCCATCCTCCCGTGTTCGCCGCAACGTTCCGCCACCGGCGAAAGCAACGGTCACGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCATTAAATCCTGCCAGTTCAGTTCTCCCGCGCCGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAATATCTCGCCATCGTTCGATAACTTCGCGTCAAATTCGATCATCTTATGAcc                                                                                                                                                            <  1:400782/251‑1 (MQ=255)
                                                                                   gTTTGATTTCTATATTCGCCATCAGCCCGTGTTCGCGGCAACGTTCCGCCACCTTCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTTCTGTACCAACTGCCCGCATGCACGCGCAGTAAATCCTGCCATGTCAGTTCACCCGCGTCTCCCCAACCGTTTCTGGTACTTTCGTGATTTTCGTCATGGAACAGAAAGTTCTCTCCATCTTTCTATAACTTCGCGTCAAATTCGATCATCTTATGACCg                                                                                                                                                           <  1:279606/251‑1 (MQ=255)
                                                                                    tttGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGt                                                                                                                                                          >  2:269994/1‑251 (MQ=255)
                                                                                      tGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGCGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTAt                                                                                                                                                        >  1:215831/1‑251 (MQ=255)
                                                                                        aTTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTAttt                                                                                                                                                      <  1:252943/251‑1 (MQ=255)
                                                                                                    gCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTc                                                                                                                                          >  2:26017/1‑251 (MQ=255)
                                                                                                    gCCATCACCCCCTGTGCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGTTAGCGGCTCACCTTTAAACATGTTGCTGTCCCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCCCCCGCGACGCACCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAGCTTCGAGTAAAATTCGATCATCTTATGACCGTATTTTGCCACGACGTc                                                                                                                                          <  2:355098/251‑1 (MQ=255)
                                                                                                    cccATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTc                                                                                                                                          <  2:216240/250‑1 (MQ=255)
                                                                                                            cccGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACg                                                                                                                                  >  1:81566/1‑251 (MQ=255)
                                                                                                                             gTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTATACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCggg                                                                                                                 >  1:298371/1‑251 (MQ=255)
                                                                                                                                  gCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAg                                                                                                            <  2:97877/251‑1 (MQ=255)
                                                                                                                                       cggCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGGAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCCAATTCGATCATCTTATGACCGAATTTTGCACCGACGCCGATTGACGCAAGGGAGTTTTCAGTGGCCAGCTTAc                                                                                                       >  2:237338/3‑251 (MQ=255)
                                                                                                                                       cTGCGAAAGCAACGGTAGCGGCTGACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCCGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTAc                                                                                                       >  2:36931/1‑251 (MQ=255)
                                                                                                                                         gCGAAAGCACCGGTAGCTGCTCCCCTTTAAACATTTTGCTGTTCCAACTGCCCGCATCCCCGCGCATTAAATCCTGCCAGTTCAGTCCCCCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGTGGCCAGCTTAcc                                                                                                      <  2:8530/250‑1 (MQ=255)
                                                                                                                                           gAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTAccgcc                                                                                                   <  2:211030/251‑1 (MQ=255)
                                                                                                                                               gCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGTGCCAGCTTACCGCCGCCa                                                                                               <  1:330230/251‑1 (MQ=255)
                                                                                                                                                   cGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGAt                                                                                           >  1:129487/1‑251 (MQ=255)
                                                                                                                                                     gTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGa                                                                                         <  1:135486/251‑1 (MQ=255)
                                                                                                                                                       aGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGc                                                                                       >  1:151320/1‑251 (MQ=255)
                                                                                                                                                                ccTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTGTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCTATGGACGCGAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATg                                                                                <  1:394356/249‑1 (MQ=255)
                                                                                                                                                                    tAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTCGCTGGGACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTAGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAACGAAGGTGCgggg                                                                           >  2:334345/1‑250 (MQ=255)
                                                                                                                                                                        cATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATa                                                                        >  2:47934/1‑249 (MQ=255)
                                                                                                                                                                         aTTTTGCTGTACCAACTGCCCGCATCCCCGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTGCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAgg                                                                     <  2:250133/251‑1 (MQ=255)
                                                                                                                                                                         aTTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAgg                                                                     <  1:49342/251‑1 (MQ=255)
                                                                                                                                                                            ttGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGcc                                                                   >  2:200465/1‑250 (MQ=255)
                                                                                                                                                                                tGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGCGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGCGATAAGGCCAGtt                                                               >  2:334244/1‑250 (MQ=255)
                                                                                                                                                                                  tACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGCCGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGTGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACt                                                            >  2:290664/1‑251 (MQ=255)
                                                                                                                                                                                   aCCAACTGCCCGCATCCCCGCGCCGTAAATCCTGCCAGATCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCAATCTTTCGATAACCTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCCACGTCGATTGACGCCAGGTTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACt                                                            <  2:142972/250‑1 (MQ=255)
                                                                                                                                                                                     cAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCAt                                                         >  1:73853/1‑251 (MQ=255)
                                                                                                                                                                                        cTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACt                                                      >  1:287222/1‑251 (MQ=255)
                                                                                                                                                                                          gCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCg                                                    <  2:123309/251‑1 (MQ=255)
                                                                                                                                                                                             cGCATTCACGCGCAGTAAATCCTGCCAGTTCAGTTCAACCGCGACCCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTg                                                 >  2:258108/1‑251 (MQ=255)
                                                                                                                                                                                                 tCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGAACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCt                                             >  1:210398/1‑251 (MQ=255)
                                                                                                                                                                                                       cgcAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTc                                        >  1:94834/1‑250 (MQ=255)
                                                                                                                                                                                                            tAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATc                                  >  1:157728/1‑251 (MQ=255)
                                                                                                                                                                                                             aaaTCCTGCCAGTTCAGTTCACCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATc                                  <  2:94834/250‑1 (MQ=255)
                                                                                                                                                                                                               aTCCTGCCAGTTCAGTTCACCCGCGACGACCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCaa                                >  1:121787/1‑250 (MQ=255)
                                                                                                                                                                                                                        gttcagttcaCACGCGACGCCCCAACCGTTGCTGGCACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGc                        >  2:252220/1‑249 (MQ=255)
                                                                                                                                                                                                                           cagttcaCCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTg                   <  1:214868/251‑1 (MQ=255)
                                                                                                                                                                                                                           cagttcaCCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTg                   <  2:303208/251‑1 (MQ=255)
                                                                                                                                                                                                                                caCCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGAtttt              <  1:200465/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATg    >  1:375682/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cccAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCa  <  1:244217/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
CAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCA  >  NC_000913/3587682‑3588168

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