Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 372,934 C→T 5.4% intergenic (+167/+57) yaiL → / ← frmB DUF2058 family protein/S‑formylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273372,9340CT 5.4% 230.2 / 3.5 75intergenic (+167/+57)yaiL/frmBDUF2058 family protein/S‑formylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base C (39/31);  new base T (2/2);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

TGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACTCAACCCCTGATCAACCATGATTTCCGCAACGCGTTGA  >  CP009273/372682‑373204
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                              
tgtTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAAGCCGGATaa                                                                                                                                                                                                                                                          >  1:406362/1‑275 (MQ=255)
          cgcgcTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGAATTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGGAGTGTAGACCGGCTCAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCCTCTTCTCAGGCCTACGGCTTAAGGGCGCCGTGTAGGCCGGCTAAGGGGtgcccg                                                                                                                                                                                                                                                >  1:207973/1‑270 (MQ=38)
                                          acCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCa                                                                                                                                                                                                                >  1:353257/1‑275 (MQ=255)
                                                         acaAAGTGCCTGATGATTGTATGGTGTAAACATAGATAAGTCTTGCGCGGATGCCTCAGGAGACGCTATGCGCGTCTTATCAGGGCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGgccggataaggcgttcacgccgcatccggcagtcgtgcgcagatgcctgatgcgacgctaacgcgt                                                                                                                                                                                                 >  2:650044‑M2/2‑209 (MQ=255)
                                                                         tttGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTCTTATCAACGCATATTCAGTGTATTAGCGTGGTAGGCAata                                                                                                                                                                                   >  1:226382/1‑273 (MQ=255)
                                                                                          aGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACg                                                                                                                                                                                                                  >  1:15020/1‑225 (MQ=255)
                                                                                            aTAAGTCTTGCGCAGTGGCCTGATGGTACGCTATGCGCGTCTTATCAGGCCTACGGTTGATGGGCGAAGTGGAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATGGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGga                                                                                                                                                               <  1:154444/274‑1 (MQ=255)
                                                                                            aTAAGACGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTTCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCCATAAACCTg                                                                                                                                                                 >  1:255338/1‑272 (MQ=255)
                                                                                              aaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCata                                                                                                                                                                                                              <  2:15020/225‑1 (MQ=255)
                                                                                                     tcgcTGATGCCTGATGCGACGCTATGCGAGTCTTATCAGGCCTACGGTTTATGGGTGAAGTGTAGACCGGATAATGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaataa                                                                                                                                                     <  1:404979/274‑1 (MQ=255)
                                                                                                     gcgcAGAGGCATGATGCGACGCTATGCGCGTCTTATCAGACCTACGGTTTATGGGCGAAGTGTAGACCGGAGAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTGATCAGGCCTACGGCGTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATc                                                                                                                                                                                                                      <  1:467860/210‑1 (MQ=255)
                                                                                                       gcAGATGCCGGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATAtt                                                                                                                                                                                                            <  1:6588/218‑1 (MQ=255)
                                                                                                           ttGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCCACGCATATTCAATTTATTGGCTTGGTCAGCAATCTGCTCGCCCATAAAACTGGAGACAAAATAATAGCtg                                                                                                                                               <  1:429802/274‑1 (MQ=255)
                                                                                                                 gtttcGAGGCTATGCGCTTATTATAAGGCCTACGGTTTATGGGGGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCTATGACTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCa                                                                                                                                         <  1:409860/271‑1 (MQ=255)
                                                                                                                       aCGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                                                                                                     >  2:362068/1‑177 (MQ=255)
                                                                                                                           taagcgcGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACAAATTCTCGCGCCACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGt                                                                                                                                                                                                                               <  1:381356/177‑1 (MQ=17)
                                                                                                                           tATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTc                                                                                                                                                                                                                                 <  1:362075/177‑1 (MQ=255)
                                                                                                                             tGCGCGTCTTATCAGGCCTACTGTTTATGGTCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAAAAAAAGCTGTGATAATAAAACTCgcg                                                                                                                              >  1:236502/1‑272 (MQ=255)
                                                                                                                             tGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCANGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCAAATAAAACTGGAGACAAAATAATAGCTGTGAAAAAAAACCTCGTg                                                                                                                              >  1:676171/1‑274 (MQ=255)
                                                                                                                               cgcgTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTCTTATCAACGCATATTCAGTTTCTTGGCGTGGTTGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGCTCATAACCCTCATGCTa                                                                                                                           >  1:22264/1‑275 (MQ=255)
                                                                                                                                      tATCGGGCCTAAGGTGTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCGTTGCGCGTCTTATCAGGCCTACAAATTCTCGCGCCACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGt                                                                                                                                                                                                                               <  1:551611/169‑1 (MQ=11)
                                                                                                                                         cAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTAc                                                                                                                                                                                                                                                                                <  2:407133/116‑1 (MQ=255)
                                                                                                                                          agggatACGGTTTATGGGCGAAGTGTAGGCGGGATAAGGCGTTCACGCCGCCTGCGGCAGTCGTGCGCCATGACTGATGCGACGCTGTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACg                                                                                                                <  1:82572/270‑1 (MQ=255)
                                                                                                                                          aGGCCTACGGTTTATGGGCGAAGTGTAGCCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTAc                                                                                                                                                                                                                                                                                <  1:628551/115‑1 (MQ=255)
                                                                                                                                            gCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCAGCATCCGGc                                                                                                                                                                                                                                     >  1:638556/1‑156 (MQ=255)
                                                                                                                                              cTACGGTTTATGGGCGAAGTTCAGACCGGGTAAGGCTTTCACGCCGCATCCGGCAGTCGTGCTCCATGCCTGATGCTACGCTTTGCGCGTCTTATCAGTCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTc                                                                                                              <  2:676157/273‑1 (MQ=255)
                                                                                                                                                aCGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTc                                                                                                                                                                                                                                 <  2:638542/156‑1 (MQ=255)
                                                                                                                                                              atccgtAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATc                                                                                                                                                                                                                      <  2:651932/149‑1 (MQ=255)
                                                                                                                                                               agggtAGCCCGAATCATGCGTTGCCGCCGCATCCGTCAGTCGCTCGCCATGCCTTATTGGACGTTTTGCGCTTTTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTAGTGCTATTATAAAAGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGAGCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc                                                                                           <  2:632500/272‑1 (MQ=255)
                                                                                                                                                                           aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTAc                                                                                                                                                                                                                                                                                >  1:74528/1‑82 (MQ=255)
                                                                                                                                                                           aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCNATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCa                                                                                >  1:540962/1‑274 (MQ=255)
                                                                                                                                                                           aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCaa                                                                               >  1:639836/1‑275 (MQ=255)
                                                                                                                                                                               ggCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGc                                                                                                                                                                                                                                                                            <  2:74527/82‑1 (MQ=255)
                                                                                                                                                                                    tCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                                                                                                     >  2:327944/1‑116 (MQ=255)
                                                                                                                                                                                    tCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAataata                                                                                                                                                   >  2:86275/1‑198 (MQ=255)
                                                                                                                                                                                       ggccgccgccGGCNGTCGTGCGCTCTGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTc                                                                    <  2:540952/273‑1 (MQ=255)
                                                                                                                                                                                        gccgcATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCtg                                                                                                                                               <  1:86276/198‑1 (MQ=255)
                                                                                                                                                                                        gccgcACCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTc                                                                                                                                                                                                                                 <  1:327950/116‑1 (MQ=255)
                                                                                                                                                                                         ccgcACCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTNGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGACGTGTTCATAACCCTCGTGATAACGGAtt                                                                                                                    <  1:514493/224‑1 (MQ=255)
                                                                                                                                                                                            cATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAg                                                                                                                                                                                                                                   >  2:318092/1‑110 (MQ=255)
                                                                                                                                                                                                cGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGt                                                                                                                                                                                                                               <  1:318098/110‑1 (MQ=255)
                                                                                                                                                                                                        tgcgcagaTGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAATTCTCACGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTa                                                   >  1:358577/8‑274 (MQ=25)
                                                                                                                                                                                                         tGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGTGCTAGTATCAACGCATATTCAGTCTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAAGAATAGCTGTGATCATAACCCTCGTGCTAACGGATGAACGGATTGAAATTAATCTCCTGGCAGCCCTTTTCACGACTTGGAGTCAGCAGCTggcccgatgtc                                                  >  2:110040/1‑264 (MQ=37)
                                                                                                                                                                                                               aTGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGNCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAatcat                                           <  1:562010/275‑1 (MQ=255)
                                                                                                                                                                                                                tGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCCATATGCTCGCCAATAAAACTGGGGACa                                                                                                                                                           >  2:320222/1‑162 (MQ=255)
                                                                                                                                                                                                                     gATGCGACGCTTTGAGCGTCTTATCAGGCCTACGGCTTTCGGGCGAGgtgt                                                                                                                                                                                                                                                                     >  1:288449/1‑51 (MQ=255)
                                                                                                                                                                                                                      aTGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGgccggataaggcgttcacgccgcatccggcagtcgtgcgcagatgcctgatgcgacgctaac                                                                                                                                                                                                     <  1:547884‑M2/114‑63 (MQ=255)
                                                                                                                                                                                                                      aTGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACg                                                                                                                                                                                                                  >  2:446372/1‑101 (MQ=255)
                                                                                                                                                                                                                      aTGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAATCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGAAGCTGTTCTGCGTAAAAATCATCACTCaa                                    >  1:415347/1‑275 (MQ=255)
                                                                                                                                                                                                                      aTGCGACGCTAACGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCtat                                                                                                                                                                                                                          >  1:350341/1‑92 (MQ=14)
                                                                                                                                                                                                                      aTGCGACGCTAACGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGAGTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTACGTAAAAATCATCACTCCAc                                   >  1:188501/1‑275 (MQ=255)
                                                                                                                                                                                                                       tGCGACGCTTTGCGCGTCTCATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTTTTCTGCGTaaaaacggtctcttat                                    >  1:581431/1‑263 (MQ=255)
                                                                                                                                                                                                                          gACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGgccggataaggcgttcacgccgcatccggcagtcgtgcgcagatgcctgatgcgacgctaacgcgt                                                                                                                                                                                                 >  2:547873‑M2/1‑48 (MQ=255)
                                                                                                                                                                                                                          gACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCata                                                                                                                                                                                                              <  1:446379/101‑1 (MQ=255)
                                                                                                                                                                                                                            cGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAataata                                                                                                                                                   >  2:588898/1‑158 (MQ=255)
                                                                                                                                                                                                                            cGCTTTGCGCGTCTTATCAGGCCTACGGATTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAc                                                                                                                                                                   >  1:531414/1‑142 (MQ=255)
                                                                                                                                                                                                                                ttGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGga                                                                                                                                                               <  2:531405/142‑1 (MQ=255)
                                                                                                                                                                                                                                ttGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCtg                                                                                                                                               <  1:588911/158‑1 (MQ=255)
                                                                                                                                                                                                                                 tGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCTGGATAAGGCGtt                                                                                                                                                                                                                                                    >  1:244177/1‑56 (MQ=255)
                                                                                                                                                                                                                                      gTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGgccggataaggcgttcacgccgcatccggcagtcgtgcgcagatgcctgatgcgacgctaacgcgtcttatcagtcctacaaattcccgcgcca                                                                                                                                                                     <  2:523749‑M2/130‑95 (MQ=255)
                                                                                                                                                                                                                                       tCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCtat                                                                                                                                                                                                                          >  1:736424/1‑76 (MQ=255)
                                                                                                                                                                                                                                        cTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGtt                                                                                                                                                                                                                                                    >  2:154354/1‑49 (MQ=255)
                                                                                                                                                                                                                                        cTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                                                                                                                      >  1:414629/1‑63 (MQ=255)
                                                                                                                                                                                                                                        cTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGCCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGAGTGGAAGGCAATATTCTCGCCAATACAACTGGAGACAAAATAATAGCTGTCATCATTACCCTCGTGATACCGGAt                                                                                                                     >  1:479402/1‑176 (MQ=37)
                                                                                                                                                                                                                                          tATCAGCCCTACGGTTTATGGGCGAAGTGTAGgccggataaggcgttcacgccgcatccggcagtcgtgcgcagatgcctgatgcgacgctaacgcgtctaatcacgcctacaacttcccgcgccatccg                                                                                                                                                                 >  1:523758‑M2/1‑32 (MQ=255)
                                                                                                                                                                                                                                           aTCAGGCCTACGGTTTATGGGCGAAGTGTAGaccggataaggcgttcacgccgcatccggcagtcgtgcgcagatgcctgatgcgacgctaac                                                                                                                                                                                                     <  2:294428‑M2/93‑63 (MQ=255)
                                                                                                                                                                                                                                           aTCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATc                                                                                                                                                                                                                      <  2:736410/76‑1 (MQ=255)
                                                                                                                                                                                                                                            tCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                                                                                                                                <  1:154355/49‑1 (MQ=255)
                                                                                                                                                                                                                                            tCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                                                                                                                                                                  <  2:414622/63‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCCTACGGTTTATGGGCGAAGTGTAGaccggataaggcgttcacgccgcatccggcagtcgtgcgcagatgcctgatgcgacgctaacgcgt                                                                                                                                                                                                 >  1:294433‑M2/1‑27 (MQ=255)
                                                                                                                                                                                                                                                gCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCaa                                                                                                                                                                                     >  1:85526/1‑104 (MQ=255)
                                                                                                                                                                                                                                                   tACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCt                                                                                     >  2:299082/1‑197 (MQ=255)
                                                                                                                                                                                                                                                    aCGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATg                                                                                                                                                                                 <  2:85525/104‑1 (MQ=255)
                                                                                                                                                                                                                                                       gCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc                                                                                 <  1:299087/197‑1 (MQ=255)
                                                                                                                                                                                                                                                          tACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCa                                                                                                                                                                                                                >  2:59238/1‑67 (MQ=255)
                                                                                                                                                                                                                                                          tACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGa                                                                                                          >  1:431128/1‑169 (MQ=255)
                                                                                                                                                                                                                                                          tACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACTCAACCCCTGATCAACCATGATTTCCGCAACGCGTTGa  >  2:95343/1‑273 (MQ=255)
                                                                                                                                                                                                                                                            aCGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCtat                                                                                                                                                                                                                          >  2:207890/1‑56 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                              
TGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACTCAACCCCTGATCAACCATGATTTCCGCAACGCGTTGA  >  CP009273/372682‑373204

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