Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,609,235 G→A 7.9% R397H (CGT→CAT)  deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,609,2350GA 7.9% 237.4 / 7.9 83R397H (CGT→CAT) deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base A (4/3);  total (42/41)
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.80e-01

TGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTG  >  CP009273/4608992‑4609492
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                 
tGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGTAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATTGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACCCCATCAATTACAGCAGCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCATCCGCGGGCGGTTATCCCCTCGAAAAACGaa                                                                                                                                                                                                                                    >  1:696665/1‑275 (MQ=255)
      cGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAaaccac                                                                                                                                                                                                                              >  1:185057/1‑275 (MQ=255)
                        cctattgctttGAGAATTAAGCTAAGTATCGGCCGACATCGATGCTGACGAAATCAGTCTATTCGTATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTTCAATGTGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGTTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGaaa                                                                                                                                                                                                            <  1:645000/266‑1 (MQ=255)
                            cTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAACGACGAAAACAACTGGCAGGAAGCGGCGAAAGCgg                                                                                                                                                                                                        >  1:415381/1‑275 (MQ=255)
                                      aaCTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGc                                                                                                                                                                                              >  1:192169/1‑275 (MQ=255)
                                           cgcgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTa                                                                                                                                                                                         >  1:392814/1‑275 (MQ=255)
                                                       gaggagagggATGGTGAAGAGAGANGTCTATGCGGGTAACGAGGGTTGTGGGAGTCAAGTGGAGCCCCGCGCGTGGGGGATGGGAGGGGTTGCGATGGGCGGCGGACGGCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATa                                                                                                                                                                             <  2:415374/266‑1 (MQ=255)
                                                                                   tATGCTGATACCGAAGTTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTAGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCATc                                                                                                                                                                                                                                                                 >  2:658591/1‑163 (MQ=255)
                                                                                   tATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGTGTCTGGGCGGCCATGTAGATGGTCAGCGTCCGCTGGCGGTTATCCACGCGACAGACGCAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATCAAACTTGCCGATAAAGCACCGCAAAGCACACCAACAGACTa                                                                                                                                                 >  1:634871/1‑275 (MQ=255)
                                                                                   tATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCATc                                                                                                                                                                                                                                                                 >  2:390168/1‑163 (MQ=255)
                                                                                   tATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCATCGCTGTCggca                                                                                                                                                                                                                                                      >  2:473489/1‑171 (MQ=255)
                                                                                   tATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCCCGCGCGGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTTGACTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGCCGGTCAGCGTc                                                                                                                                                                                                                                                                 >  1:146805/1‑163 (MQ=255)
                                                                                       cTGATACCGGAGGTTTTGGGAGTGAAATGGCTACCAGCGCGATGGGGATGGCAGTGGGTGCAATGAGAGGCGTATGCCGTACGGCATCTGACACCTGCGATTAAAGCGTCGGCTTTACTGATATGGCGCGGCGGGGAGACCAGGGTGACGGTCAGCGTcgctg                                                                                                                                                                                                                                                             <  2:146804/163‑5 (MQ=255)
                                                                                       cTGATACCGAAGTTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTAGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCATcgctg                                                                                                                                                                                                                                                             <  1:658600/163‑5 (MQ=255)
                                                                                       cTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGCCCAGGTATACGGTCAGCATCGCTGTCggcagata                                                                                                                                                                                                                                                  <  1:473496/173‑7 (MQ=255)
                                                                                       cTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCATcgctg                                                                                                                                                                                                                                                             <  1:390174/163‑5 (MQ=255)
                                                                                           tACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGAAATTAAACTTGCCGATAAAGCACCGGAAAGCACAACAACTGTCTCTAGCCGTa                                                                                                                                         >  1:259617/1‑275 (MQ=255)
                                                                                                   gTTTTGTCAGTGATATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                                                                                 <  1:448579/275‑1 (MQ=255)
                                                                                                             tGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGACGTATCAGCAAATAACGGgata                                                                                                                       >  1:280366/1‑275 (MQ=255)
                                                                                                              ggactGGGTACCCGCGCTGGGGGGGTGGCAGTGGTTGCTCGTGGGTGCGCACGCCGTCAGGCTTCTAACGGCGTTGATTACAGCGTCGGCTTGACTGAGATGGCGCGTGTGGGCGACCAGTTAGACGGTCAGCGTCCGCTGTCGGTTATCCACGCGAAAGACGAAAACAACTGTCAGGAAGCGGCGAAATCGGTGAAAGCGGCAATTAAACTTGCCTATAAAGCACCGGAAAGCACACCAACTTTCTATCGCCGTATCAGCTAATAACGGtata                                                                                                                       <  2:696656/270‑1 (MQ=255)
                                                                                                                       aCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGCCGGTCAGCGTCCGCTGTCGGTTATCCACGCGAAAGACGAAAACAACTTGCAGGAAGCGGCGCACGCGGTGAAAGCGGCAATTACACTTGCCGATAAAGCACCGGCAAGACCAACAACCGGCCAGCGCCGCAGNGGCCAATACCggggtgcggcggggg                                                                                                              >  2:481465/1‑261 (MQ=255)
                                                                                                                           gcgcggTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGNTTACAGCGNCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                                                                            <  1:481472/256‑1 (MQ=255)
                                                                                                                            cgcgcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGaa                                                                                                                                                                  >  1:217854/1‑217 (MQ=255)
                                                                                                                                cTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGca                                                                                                                                                              <  2:217849/217‑1 (MQ=255)
                                                                                                                                                    acgcggcggnGGACGCCGGCAGGCATCTGACACCATCTATTACAGCGTCGGCTTTCCTGATATTGCGCGTGTGGGCGACCATGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                                                                                <  2:280361/271‑1 (MQ=255)
                                                                                                                                                     ctacgcagcggACGCCTACAGGCACCTGACAGCATGCCTTACAGCGTCGGCTTTTCTGATATGGCGCTGCTGGGCGACCAGGTAGACGGTCAGCGGCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                                                                                <  2:185054/267‑1 (MQ=255)
                                                                                                                                                     aTGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTACGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                                                                                <  1:382853/274‑1 (MQ=255)
                                                                                                                                                      tGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                                                                        <  1:176044/201‑1 (MQ=255)
                                                                                                                                                        ggggtccannaccGCCAAGCATCTGCCTCCATCGATTCGCGCCGCTGCTTTGTTGAGATCGTGGGTCTGGGCGACCAGGTAGACTGTCAGCGTCCTCTGGCGGTGAGCCACGCGAAAGACGAAAACAACGGGCAGGAAGCGGCGAAATCGGTGAAAGCGGCAATTAAACGTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAGTTAAATTTGAAGCACTCAGTACGGTGaa                                                                              <  2:75183/262‑1 (MQ=255)
                                                                                                                                                        ggcggcGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAAt                                                                                                                                >  2:659799/1‑223 (MQ=255)
                                                                                                                                                         gcggcgGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGt                                                                                                                                                                                                                                                                  >  2:172783/1‑92 (MQ=255)
                                                                                                                                                            gcggACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                                                                            <  1:659808/223‑1 (MQ=255)
                                                                                                                                                            gcggACGCCGTCAGGCACCTGACACCATCGATTCCAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCATGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGGATCa                                                                                                                                      <  1:699083/213‑1 (MQ=255)
                                                                                                                                                             cggACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGc                                                                                                                                                                                                                                                              <  1:172786/92‑1 (MQ=255)
                                                                                                                                                             cggACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAGGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCACTTAAACTTGACGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGAGCATAACGGTATACTGCTCGGATCATTGAAATTTGAAGAACTGAGTACGGAGAACatat                                                                         >  1:75184/1‑273 (MQ=255)
                                                                                                                                                                  gCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACGGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCANACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTaa                                                                                    >  1:556519/1‑256 (MQ=255)
                                                                                                                                                                    cGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctgat                                                                                                             >  2:586215/1‑230 (MQ=255)
                                                                                                                                                                        aGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAtt                                                                                                         <  1:586224/230‑1 (MQ=255)
                                                                                                                                                                            aTCTGACACCATCGATTTCAGCTTTGTCTTTATTGTTATTGTTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCTAAAGACGAAAACAACTTGCAGGAAGCGGNGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCATCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCCCTGAGTACGGAGAACATATTAAACGTGCATttatta                                                        <  1:610983/275‑1 (MQ=255)
                                                                                                                                                                              cTGCCACCATCGCTTACAGCGCCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCCGCGACAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGCGACAGCGGCAATTAAACTAGCCGCTAAAGCAcc                                                                                                                                                                      >  2:310709/1‑163 (MQ=255)
                                                                                                                                                                                  cacCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCGGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGaa                                                                                                                                                                  <  1:310715/163‑1 (MQ=255)
                                                                                                                                                                                       tCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAg                                                                                                                                     >  1:458385/1‑187 (MQ=255)
                                                                                                                                                                                           ttACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                                                                                 <  2:458378/187‑1 (MQ=255)
                                                                                                                                                                                               aGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGaaa                                                                                                                                                                 >  1:679731/1‑151 (MQ=255)
                                                                                                                                                                                               aGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGa                                                                                                                                  >  2:521107/1‑182 (MQ=255)
                                                                                                                                                                                                 cGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCgg                                                                                                                                                                                               >  1:643105/1‑119 (MQ=255)
                                                                                                                                                                                                   tCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGcaca                                                                                                                                                            >  1:723377/1‑152 (MQ=255)
                                                                                                                                                                                                   tCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGcac                                                                                                                                                             <  2:679722/151‑1 (MQ=255)
                                                                                                                                                                                                   tCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAg                                                                                                                                     >  1:472203/1‑175 (MQ=255)
                                                                                                                                                                                                   tCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAg                                                                                                                                     >  2:576382/1‑175 (MQ=255)
                                                                                                                                                                                                   tCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGa                                                                                                                                  >  1:556783/1‑178 (MQ=255)
                                                                                                                                                                                                   tCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                                                                              <  1:521114/182‑1 (MQ=255)
                                                                                                                                                                                                   tCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgat                                                                                                                  >  1:646264/1‑194 (MQ=255)
                                                                                                                                                                                                    cGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACgg                                                                                                                           >  2:121765/1‑184 (MQ=255)
                                                                                                                                                                                                     ggCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAAt                                                                                                                                                                                           <  2:643096/119‑1 (MQ=255)
                                                                                                                                                                                                       cTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                                                                        <  2:723367/152‑1 (MQ=255)
                                                                                                                                                                                                       cTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                                                                                 <  1:576391/175‑1 (MQ=255)
                                                                                                                                                                                                       cTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                                                                                 <  2:472196/175‑1 (MQ=255)
                                                                                                                                                                                                       cTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                                                                              <  2:556774/178‑1 (MQ=255)
                                                                                                                                                                                                       cTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctga                                                                                                              <  2:646255/194‑1 (MQ=255)
                                                                                                                                                                                                        tttACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGtata                                                                                                                       <  1:121766/184‑1 (MQ=255)
                                                                                                                                                                                                         ttACAGATATGGCGCGTCAGGGCGACCAGGTCGACGGTCAGCGTCCGCTGGCGGTTATCAACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGgcg                             <  2:192166/273‑1 (MQ=255)
                                                                                                                                                                                                                 atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAaacaa                                                                                                                                                                                                                               >  1:102409/1‑71 (MQ=255)
                                                                                                                                                                                                                 atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGcc                                                                                                                                                                                 >  1:475117/1‑117 (MQ=255)
                                                                                                                                                                                                                 atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACagaaga                   >  1:183072/1‑275 (MQ=255)
                                                                                                                                                                                                                 atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGAACTCATTCGGCATCGGCGCTACagaaga                   >  1:317488/1‑275 (MQ=255)
                                                                                                                                                                                                                     cgcgTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTgg                                                                                                                                                                                                                           <  2:102408/71‑1 (MQ=255)
                                                                                                                                                                                                                     cgcgTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATa                                                                                                                                                                             <  2:475110/117‑1 (MQ=255)
                                                                                                                                                                                                                     cgcgTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTAc                                                                                    >  2:398748/1‑206 (MQ=255)
                                                                                                                                                                                                                     cagggcggggCGAAAAGGTAGACGGTCACCGTCCGCGGGCGGTTATCCACGCGAAAGAAGAAAAAAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCa               <  2:614016/268‑1 (MQ=255)
                                                                                                                                                                                                                       cgTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCgt                                                                                                                                             >  1:634500/1‑146 (MQ=255)
                                                                                                                                                                                                                         tCTGGGCGACCAGGTAGACGGTCAGCGTTCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                                                                                <  1:398755/206‑1 (MQ=255)
                                                                                                                                                                                                                         tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATGCACGCGAAAGCCGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTTCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCCGCGAATAACGGAATACGGATCTGATCATCTAACTTTGAAGCACTAAGTAAGGAGAACAGATGAAACGTGCCTTTATTATGGTGCTGAACTCATTCGTCATAGGCGCTACatacaattctaaa            >  2:588469/1‑262 (MQ=255)
                                                                                                                                                                                                                         tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGcaca                                                                                                                                                            >  1:228784/1‑130 (MQ=255)
                                                                                                                                                                                                                         tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTATATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTg                                                                                                 >  2:663463/1‑189 (MQ=255)
                                                                                                                                                                                                                           tGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGaa                                                                      <  1:333344/214‑1 (MQ=255)
                                                                                                                                                                                                                             ggCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                                                                        <  2:228779/130‑1 (MQ=255)
                                                                                                                                                                                                                             ggCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTATATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGc                                                                                             <  1:663472/189‑1 (MQ=255)
                                                                                                                                                                                                                             aacGACCAGGTAGACGGTCAGCGGCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCt        <  2:183069/272‑1 (MQ=255)
                                                                                                                                                                                                                              gCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTg                                                                                                                                                                                   >  2:450546/1‑102 (MQ=255)
                                                                                                                                                                                                                                  ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGa                                                                                                                                                                               <  1:450553/102‑1 (MQ=255)
                                                                                                                                                                                                                                  ccAGGTAGACGGTCAGCGTCCGCTGGCGGCTCTCCTCGCGCAAGACGAAAACAAATGGCAGTAAGCAGCGAAAGCGGTGAACGCGGCACTTAAAATTGTCGATCAAGCACAGGAAAGCACACCAACTGTCTATCGCCTTATCAGCGAATAACGGTATACTGATCTTATCATTGAAATTTGAAGCACTGAGTAAGGAGAACATATGAAACGTGCATTTAAGCTGGTGGTGGACTCATCAGGCAACGACGCTCACGACGCTGCACAAAGCTTTGGTg  >  2:625842/1‑275 (MQ=255)
                                                                                                                                                                                                                                   cAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGNTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACagaaga                   >  1:614025/1‑257 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                 
TGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTG  >  CP009273/4608992‑4609492

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: