Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,451 A→G 5.0% intergenic (+75/‑403) fimB → / → fimE tyrosine recombinase/inversion of on/off regulator of fimA/tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,531,4510AG 5.0% 1145.6 / 34.2 338intergenic (+75/‑403)fimB/fimEtyrosine recombinase/inversion of on/off regulator of fimA/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (152/169);  new base G (8/9);  total (160/178)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 346 total aligned reads displayed.

ACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTC  >  CP009273/4531302‑4531600
                                                                                                                                                     |                                                                                                                                                     
 ccGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttat                                                                                                                                                      <  1:696345/150‑1 (MQ=255)
     cAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                                                                                                                  >  1:64426/1‑150 (MQ=255)
       gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                                                                                                                <  1:856435/150‑1 (MQ=255)
          aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                                                                                                             <  2:4239774/150‑1 (MQ=255)
          aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                                                                                                             <  2:4239601/150‑1 (MQ=255)
          aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                                                                                                             >  2:1341566/1‑150 (MQ=255)
          aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                                                                                                             <  1:871711/150‑1 (MQ=255)
          aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                                                                                                             <  1:4128568/150‑1 (MQ=255)
            gCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGtt                                                                                                                                           <  1:767799/150‑1 (MQ=255)
              aGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAt                                                                                                                                         <  1:2389212/150‑1 (MQ=255)
              aGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAt                                                                                                                                         >  2:3545564/1‑150 (MQ=255)
              aGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAt                                                                                                                                         <  1:1473360/150‑1 (MQ=255)
               gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACCGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAtt                                                                                                                                        <  1:3229404/150‑1 (MQ=255)
               gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAtt                                                                                                                                        <  1:952755/150‑1 (MQ=255)
                ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTg                                                                                                                                       <  1:29489/150‑1 (MQ=255)
                  cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTa                                                                                                                                     <  1:37583/150‑1 (MQ=255)
                  cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTa                                                                                                                                     <  1:36611/150‑1 (MQ=255)
                    tttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAg                                                                                                                                   >  2:2159004/1‑150 (MQ=255)
                    tttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGACTAACCATTTATTAATTGATGTTATTGTAAg                                                                                                                                   >  1:2464151/1‑150 (MQ=255)
                        tACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGtgtg                                                                                                                               <  2:4523505/150‑1 (MQ=255)
                         aCGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGa                                                                                                                              <  1:1105550/150‑1 (MQ=255)
                         aCGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGa                                                                                                                              <  1:29042/150‑1 (MQ=255)
                          cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAt                                                                                                                             >  1:2422687/1‑150 (MQ=255)
                           ggCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAtt                                                                                                                            >  1:2644915/1‑150 (MQ=255)
                             cATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCt                                                                                                                          >  2:707596/1‑150 (MQ=255)
                              aTCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCtt                                                                                                                         <  2:3059327/150‑1 (MQ=255)
                               tCTGGGATAGAGCCAGAGGACGACAGCGTCACTCTGTTTTATAGAGATAGTTGAGATACCAGGGATGTTGTTAATAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTTATTGATGTTATTGTAAGTGTGTTTCTTa                                                                                                                        <  1:4050700/150‑1 (MQ=255)
                                cTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTaa                                                                                                                       <  1:18155/150‑1 (MQ=255)
                                   ggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTa                                                                                                                    <  1:4447241/150‑1 (MQ=255)
                                   ggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTa                                                                                                                    >  1:2737027/1‑150 (MQ=255)
                                       agagCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCat                                                                                                                <  2:2592014/150‑1 (MQ=255)
                                            cAGAGGACGACAGCGTCACGCTGTTTTATATAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCAtaataa                                                                                                           <  2:2708629/150‑1 (MQ=255)
                                            cAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCAtaataa                                                                                                           >  1:1955046/1‑150 (MQ=255)
                                               agGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATaaaaa                                                                                                        >  2:3999006/1‑148 (MQ=255)
                                                  acgacAGCGTCACGCTGTTTTATAGAGATAGTTTAGATACCTTGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTa                                                                                                     <  1:394339/150‑1 (MQ=255)
                                                  acgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAATCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTa                                                                                                     <  1:4304863/150‑1 (MQ=255)
                                                  acgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTa                                                                                                     <  1:3255607/150‑1 (MQ=255)
                                                   cgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTaa                                                                                                    <  1:3248196/150‑1 (MQ=255)
                                                        gCGTCACGCTGTATTATAGAGATAGTTGAGATACCAGGGATTGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTaa                                                                                               <  2:905343/150‑1 (MQ=255)
                                                           tCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCa                                                                                            >  2:1540757/1‑150 (MQ=255)
                                                            cACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCat                                                                                           >  2:2300902/1‑150 (MQ=255)
                                                              cGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCatat                                                                                         >  1:3501834/1‑150 (MQ=255)
                                                               tcTGTTTTATATATATAGTTGAGATTCCAGGGATGTTGTTAAAAACTACATTATTCTTAAGATATATAAAATATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACTTTAAGTTAACCATATc                                                                                        <  2:2460209/149‑1 (MQ=255)
                                                               gCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATc                                                                                        >  1:3390880/1‑150 (MQ=255)
                                                               gCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATc                                                                                        <  1:1853627/150‑1 (MQ=255)
                                                                  gTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCAt                                                                                     >  1:2588056/1‑150 (MQ=255)
                                                                    tttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATAc                                                                                   >  1:1155280/1‑150 (MQ=255)
                                                                     ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTGACCATATCCATACa                                                                                  >  2:3392524/1‑150 (MQ=255)
                                                                      tataGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACaa                                                                                 <  2:1023582/150‑1 (MQ=255)
                                                                       ataGAGATAGTTGAGATACCAGGGATGGTTTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACaaa                                                                                <  2:2607042/150‑1 (MQ=255)
                                                                        taGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACaaaa                                                                               >  2:338179/1‑150 (MQ=255)
                                                                             aTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAATTAACCATATCCATACAAAATACaa                                                                          <  2:2296845/150‑1 (MQ=255)
                                                                                     gaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTAt                                                                  >  1:316159/1‑150 (MQ=255)
                                                                                      aTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATg                                                                 <  1:4212681/150‑1 (MQ=255)
                                                                                      aTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATg                                                                 <  1:2881262/150‑1 (MQ=255)
                                                                                        aCCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGtt                                                               >  2:371051/1‑150 (MQ=255)
                                                                                          cAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttct                                                             <  2:2618044/150‑1 (MQ=255)
                                                                                           agtgATTGTGTTAAAAACTACATTATTATTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                            <  2:4350722/147‑1 (MQ=255)
                                                                                           aGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                            <  1:4521191/150‑1 (MQ=255)
                                                                                           aGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                            <  1:2282310/150‑1 (MQ=255)
                                                                                             ggATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTTATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCa                                                          <  2:2892616/150‑1 (MQ=255)
                                                                                             ggATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAATTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCa                                                          >  1:3018927/1‑150 (MQ=255)
                                                                                              gATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCaa                                                         >  1:1281942/1‑150 (MQ=255)
                                                                                               aTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCaaa                                                        >  1:4374583/1‑150 (MQ=255)
                                                                                                 ggTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaat                                                      <  2:3622164/150‑1 (MQ=255)
                                                                                                  gtgtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaata                                                     >  1:1801208/1‑150 (MQ=255)
                                                                                                  gtgtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaata                                                     <  2:2070830/150‑1 (MQ=255)
                                                                                                   tgtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaataa                                                    <  1:4271506/150‑1 (MQ=255)
                                                                                                   tgtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAATTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaataa                                                    >  1:2361864/1‑150 (MQ=255)
                                                                                                       aaaaaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAataaataa                                                >  1:2133570/1‑150 (MQ=255)
                                                                                                        aaaaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATATTCTTCAAAAtaaataaa                                               >  2:1209371/1‑150 (MQ=255)
                                                                                                         aaaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAAc                                              >  2:2375296/1‑150 (MQ=255)
                                                                                                          aaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACa                                             >  1:354312/1‑150 (MQ=255)
                                                                                                               cATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAATGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtc                                        >  2:4107267/1‑150 (MQ=255)
                                                                                                                attattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtca                                       >  2:3819474/1‑150 (MQ=255)
                                                                                                                attattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtca                                       >  1:1617736/1‑150 (MQ=255)
                                                                                                                 ttattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcat                                      >  1:1847546/1‑150 (MQ=255)
                                                                                                                  tattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcatt                                     >  2:833796/1‑150 (MQ=255)
                                                                                                                  tattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcatt                                     >  2:4396854/1‑150 (MQ=255)
                                                                                                                  tattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcatt                                     <  1:167502/150‑1 (MQ=255)
                                                                                                                        tACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATaaa                               >  2:2988412/1‑150 (MQ=255)
                                                                                                                          cGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAtt                             <  1:2333616/150‑1 (MQ=255)
                                                                                                                          cGAGATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATaaagt                             <  2:3839793/150‑3 (MQ=255)
                                                                                                                           gATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAttt                            >  2:2532025/1‑150 (MQ=255)
                                                                                                                            atatatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTa                           >  1:1018858/1‑150 (MQ=255)
                                                                                                                              atatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATTAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTaca                         <  2:1783074/150‑1 (MQ=255)
                                                                                                                                 tCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACAt                      <  2:1113987/150‑1 (MQ=255)
                                                                                                                                 tCAAAAATGATTAACAATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACAt                      <  2:3130015/150‑1 (MQ=255)
                                                                                                                                    aaaaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCAc                   <  1:1306969/150‑1 (MQ=255)
                                                                                                                                     aaaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACt                  >  1:3299350/1‑150 (MQ=255)
                                                                                                                                           ttAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAt            >  1:597968/1‑150 (MQ=255)
                                                                                                                                           ttAAACATTTATTAATTGATGTTATTGTAAGTTTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAt            <  2:2773502/150‑1 (MQ=255)
                                                                                                                                            tAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtc           <  1:750453/150‑1 (MQ=255)
                                                                                                                                            tAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtc           >  2:2325730/1‑150 (MQ=255)
                                                                                                                                             aaCCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtct          >  2:1858453/1‑150 (MQ=255)
                                                                                                                                             aaCCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtct          >  2:3649074/1‑150 (MQ=255)
                                                                                                                                               ccATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTcc        <  1:3203663/150‑1 (MQ=255)
                                                                                                                                               ccATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAAATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTcc        <  2:3908786/150‑1 (MQ=255)
                                                                                                                                                 aTTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAATTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTg      >  1:1544010/1‑150 (MQ=255)
                                                                                                                                                  tttATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGt     <  2:1130947/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTC  >  CP009273/4531302‑4531600

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: