Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,513,424 A→T 46.8% V213V (GTT→GTA hofQ ← DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,513,4240AT46.8% 84.9 / inf 363V213V (GTT→GTAhofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base A (89/104);  new base T (79/91);  total (168/195)
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.35e-01

Only 100 of 369 total aligned reads displayed.

GGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTGCGCTACCCACTGTTCAAGCGCGCTTAA  >  CP009273/3513275‑3513573
                                                                                                                                                     |                                                                                                                                                     
  aCTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTAc                                                                                                                                                     <  1:4648741/150‑1 (MQ=255)
    tGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGc                                                                                                                                                   >  2:4141667/1‑150 (MQ=255)
    tGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGc                                                                                                                                                   >  1:3595100/1‑150 (MQ=255)
     ggCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGcc                                                                                                                                                  <  2:3473634/150‑1 (MQ=255)
      gCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGnca                                                                                                                                                 <  1:2963398/150‑1 (MQ=255)
      gCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGcca                                                                                                                                                 <  1:573529/150‑1 (MQ=255)
      gCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGcca                                                                                                                                                 <  1:1603723/150‑1 (MQ=255)
       cGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGccac                                                                                                                                                >  2:4180367/1‑150 (MQ=255)
         ataataTCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCATTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTTCCAAGCGTGGTGACTTGCCCTACGccacca                                                                                                                                              <  2:2768145/150‑1 (MQ=255)
          taataTCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAg                                                                                                                                             >  2:2813072/1‑150 (MQ=255)
          taataTCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCATCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAg                                                                                                                                             <  2:1475789/150‑1 (MQ=255)
           aataTCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGc                                                                                                                                            >  2:3742422/1‑150 (MQ=255)
              atCCAGCTGCTGTTTTTGTTCGAGCGCGGACAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCgtg                                                                                                                                         >  2:3583862/1‑150 (MQ=255)
              atCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCgtg                                                                                                                                         >  1:1822139/1‑150 (MQ=255)
                   gctgctGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTgag                                                                                                                                    >  1:2152099/1‑148 (MQ=255)
                       ctgTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATc                                                                                                                                <  2:3322252/150‑1 (MQ=255)
                        tgTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCg                                                                                                                               <  1:1154302/150‑1 (MQ=255)
                        tgTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCg                                                                                                                               <  1:3152168/150‑1 (MQ=255)
                           ttttGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTTCGCCACCAGCGTGTTGCGCATCGGcc                                                                                                                            <  1:1974383/150‑1 (MQ=255)
                           ttttGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGcc                                                                                                                            <  2:217113/150‑1 (MQ=255)
                           ttttGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGcc                                                                                                                            >  1:3520350/1‑150 (MQ=255)
                            tttGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCa                                                                                                                           >  2:482650/1‑150 (MQ=255)
                                  cGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTcc                                                                                                                     <  1:1564697/150‑1 (MQ=255)
                                  cGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTcc                                                                                                                     <  1:3013374/150‑1 (MQ=255)
                                   gAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCa                                                                                                                    <  1:3445863/150‑1 (MQ=255)
                                   gAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCa                                                                                                                    <  1:1042730/150‑1 (MQ=255)
                                      cncgGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCAttt                                                                                                                 >  1:3807371/1‑150 (MQ=255)
                                      cgcgGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCAttt                                                                                                                 >  1:2683246/1‑150 (MQ=255)
                                       gcgGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTc                                                                                                                <  2:4269956/150‑1 (MQ=255)
                                        cgGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCa                                                                                                               <  1:1330328/150‑1 (MQ=255)
                                            aaGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGcc                                                                                                           >  2:2070759/1‑150 (MQ=255)
                                            aaGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGcc                                                                                                           >  1:4041080/1‑150 (MQ=255)
                                             aGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCt                                                                                                          <  2:3542694/150‑1 (MQ=255)
                                              gCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTa                                                                                                         <  2:1866793/150‑1 (MQ=255)
                                              gCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTa                                                                                                         <  1:3934489/150‑1 (MQ=255)
                                              gCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTa                                                                                                         <  1:2801520/150‑1 (MQ=255)
                                              gCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTa                                                                                                         >  2:4530407/1‑150 (MQ=255)
                                               cTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTaa                                                                                                        >  2:2312290/1‑150 (MQ=255)
                                               cTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTaa                                                                                                        <  1:4180162/150‑1 (MQ=255)
                                               cTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTaa                                                                                                        <  2:3114889/150‑1 (MQ=255)
                                                 cAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACt                                                                                                      >  2:736390/1‑150 (MQ=255)
                                                 cAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACt                                                                                                      <  1:853225/150‑1 (MQ=255)
                                                  aaGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTc                                                                                                     <  2:1675224/150‑1 (MQ=255)
                                                   aGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCa                                                                                                    <  2:3704773/150‑1 (MQ=255)
                                                   aGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCa                                                                                                    <  2:4095834/150‑1 (MQ=255)
                                                    gATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCAc                                                                                                   <  2:858617/150‑1 (MQ=255)
                                                     aTCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACg                                                                                                  >  1:4423812/1‑150 (MQ=255)
                                                         aGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGATAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCAATTCACGCCTAACTCACGCaaa                                                                                              >  1:752721/1‑150 (MQ=255)
                                                           cAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACt                                                                                            <  1:2532169/150‑1 (MQ=255)
                                                           cAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACt                                                                                            <  1:140872/150‑1 (MQ=255)
                                                            aaGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACtt                                                                                           <  1:1263085/150‑1 (MQ=255)
                                                               cGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACttttt                                                                                        <  2:18318/150‑1 (MQ=255)
                                                                gTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACtttttt                                                                                       <  1:2321702/150‑1 (MQ=255)
                                                                  ccGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCa                                                                                     <  2:1203618/150‑1 (MQ=255)
                                                                   cGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCAt                                                                                    <  2:1367746/150‑1 (MQ=255)
                                                                   cGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCAt                                                                                    >  1:4168223/1‑150 (MQ=255)
                                                                   cGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCAt                                                                                    >  2:1272025/1‑150 (MQ=255)
                                                                    gTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCAtt                                                                                   >  1:4271826/1‑150 (MQ=255)
                                                                    gTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCAtt                                                                                   >  1:1692730/1‑150 (MQ=255)
                                                                     ttGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTa                                                                                  >  2:1287674/1‑150 (MQ=255)
                                                                     ttGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTa                                                                                  >  2:4405457/1‑150 (MQ=255)
                                                                     ttGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTa                                                                                  <  1:2225292/150‑1 (MQ=255)
                                                                        aTGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATg                                                                               <  1:4008026/150‑1 (MQ=255)
                                                                             cccAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGAc                                                                          <  1:3689528/150‑1 (MQ=255)
                                                                             cccAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGAc                                                                          <  2:1033850/150‑1 (MQ=255)
                                                                              ccAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACa                                                                         >  2:3960358/1‑150 (MQ=255)
                                                                              ccAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACa                                                                         <  1:4714396/150‑1 (MQ=255)
                                                                               cAATGTTAAATCCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACCCACGCAAACTTTTTTCATTAATGGTGACaa                                                                        <  2:2602100/150‑1 (MQ=255)
                                                                               cAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACCCACGCAAACTTTTTTCATTAATGGTGACaa                                                                        <  2:3035994/150‑1 (MQ=255)
                                                                                    ttAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATgc                                                                   <  1:2875069/150‑1 (MQ=255)
                                                                                     tAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATgcg                                                                  >  1:4099518/1‑150 (MQ=255)
                                                                                     tAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATgcg                                                                  >  1:1455666/1‑150 (MQ=255)
                                                                                     tAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATgcg                                                                  <  1:3320825/150‑1 (MQ=255)
                                                                                         aCCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGa                                                              >  2:676012/1‑150 (MQ=255)
                                                                                         aCCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGa                                                              >  2:643060/1‑150 (MQ=255)
                                                                                          ccGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGAc                                                             >  2:4752590/1‑150 (MQ=255)
                                                                                           cGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTCACTCACGCCAACTTTTTGCATTAATGGTGACAATATGCGCCGACa                                                            >  2:2873495/1‑150 (MQ=255)
                                                                                            gACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAg                                                           <  2:1543653/150‑1 (MQ=255)
                                                                                            gACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAg                                                           <  1:1903854/150‑1 (MQ=255)
                                                                                                  cGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAAc                                                     <  2:527914/150‑1 (MQ=255)
                                                                                                   gTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAAcc                                                    >  1:1070656/1‑150 (MQ=255)
                                                                                                         gCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCg                                              <  2:1998798/150‑1 (MQ=255)
                                                                                                         gCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCg                                              <  2:3284069/150‑1 (MQ=255)
                                                                                                            gTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGAcc                                           <  2:2922122/150‑1 (MQ=255)
                                                                                                            gTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGAcc                                           <  2:133240/150‑1 (MQ=255)
                                                                                                              cGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCgg                                         >  2:4583263/1‑150 (MQ=255)
                                                                                                              cGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCgg                                         <  1:2437530/150‑1 (MQ=255)
                                                                                                                 tACGTAGAGGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGAATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAg                                      >  1:1375218/1‑150 (MQ=255)
                                                                                                                   cGGAGAGGTCGCGACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGAt                                    <  2:4691053/150‑1 (MQ=255)
                                                                                                                      agagGTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCa                                 <  1:3381200/150‑1 (MQ=255)
                                                                                                                        agGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCAtt                               >  1:2458032/1‑150 (MQ=255)
                                                                                                                         ggTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCAttt                              <  1:201622/150‑1 (MQ=255)
                                                                                                                         ggTCGCTACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCAttt                              <  1:203561/150‑1 (MQ=255)
                                                                                                                           tCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTgc                            >  1:1745411/1‑150 (MQ=255)
                                                                                                                               tACCAAGCGTGGTGACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTGCGCTa                        >  1:2561907/1‑150 (MQ=255)
                                                                                                                               tACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTGCGCTa                        <  2:2561567/150‑1 (MQ=255)
                                                                                                                                   acgCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTGCGCTACCCa                    <  2:2488119/148‑1 (MQ=255)
                                                                                                                                         ggtgACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTGCGCTACCCACTGTTc              >  2:4338957/1‑150 (MQ=255)
                                                                                                                                          gtgACTTGCCCTACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTGCGCTACCCACTGTTCa             <  1:218093/150‑1 (MQ=255)
                                                                                                                                                  cccAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTGCGCTACCCACTGTTCAAGCGCGCt     <  1:1213989/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTGGTGACTTGCCCAACGCCACCAGCGTGTTGCGCATCGGCCAGCGTCCATTTCACGCCTAACTCACGCAAACTTTTTTCATTAATGGTGACAATATGCGCCGACAGCTCAACCTGCCCGACCGGCAGATCCATTTGCGCTACCCACTGTTCAAGCGCGCTTAA  >  CP009273/3513275‑3513573

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