Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,587,957 G→A L53F (CTC→TTC)  ugpQ ← glycerophosphodiester phosphodiesterase UgpQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,587,9570GA98.3% 188.6 / ‑4.9 59L53F (CTC→TTC) ugpQglycerophosphodiester phosphodiesterase UgpQ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/27);  total (32/27)
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.52e-01

GCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCAGCCACAGCGTGCTGCCTGCCGTCGGGCGCTC  >  NC_000913/3587712‑3588199
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
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                                                  <  2:303208/251‑1 (MQ=255)
                                                                                                                                                                                             cagttcaCCCGCGACGCCCCAACCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTg                                                  <  1:214868/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)
                                                                                                                                                                                                                  accGTTGCTGGTACGTTCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCAGcc                              >  1:253360/2‑250 (MQ=255)
                                                                                                                                                                                                                                 ttCGAGATTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCAGCCACAGCGTGctgcctgc              >  1:286360/1‑251 (MQ=255)
                                                                                                                                                                                                                                     agaTTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCAGCCACAGCGTGCTGCCTGCCGTc          <  1:384222/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gaTTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGACAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCAGCCACAGCGTGCTGCCTGCCGTCg         >  2:19839/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTTGTCGTCATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCAGCCACAGCGTGCTGCCTGCCGTCgg        <  1:269994/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gtcATGGAACAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCAGCCACAGCGTGCTGCCTGCCGTCGGGCGCTc  >  1:402635/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
GCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACCATCAAAAAGATGCAGCTGATTTTCCGCCAGATGCAGCCACAGCGTGCTGCCTGCCGTCGGGCGCTC  >  NC_000913/3587712‑3588199

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: