Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,461 G→A 24.5% intergenic (+85/‑393) 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,4610GA24.5% 552.3 / 183.6 261intergenic (+85/‑393)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 G (103/93);  new base A (39/25);  total (143/118)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 263 total aligned reads displayed.

TGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCGCACTTTTT  >  CP009273/4531313‑4531610
                                                                                                                                                    |                                                                                                                                                     
tGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGt                                                                                                                                                      <  1:2689477/150‑1 (MQ=255)
  cAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTa                                                                                                                                                    <  2:2394630/150‑1 (MQ=255)
  cAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTa                                                                                                                                                    <  1:1985800/150‑1 (MQ=255)
   aGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGAtattat                                                                                                                                                   >  1:3365308/1‑150 (MQ=255)
    gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAtt                                                                                                                                                  >  1:139332/1‑150 (MQ=255)
     ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTg                                                                                                                                                 >  2:1714747/1‑150 (MQ=255)
     ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAttt                                                                                                                                                 >  2:876626/1‑149 (MQ=255)
     ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTg                                                                                                                                                 >  2:167718/1‑150 (MQ=255)
        gTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTaa                                                                                                                                              >  1:803689/1‑150 (MQ=255)
        gTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTaa                                                                                                                                              <  2:2541431/150‑1 (MQ=255)
          ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgt                                                                                                                                            <  1:1438627/150‑1 (MQ=255)
              aCGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGa                                                                                                                                        <  2:1807404/150‑1 (MQ=255)
               cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAt                                                                                                                                       <  1:1348824/150‑1 (MQ=255)
                ggCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGAtt                                                                                                                                      <  2:158191/150‑1 (MQ=255)
                     cTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTaa                                                                                                                                 >  1:3652148/1‑150 (MQ=255)
                               gCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCataat                                                                                                                       <  2:1766092/150‑1 (MQ=255)
                                  agagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATaaa                                                                                                                    >  2:3206643/1‑150 (MQ=255)
                                      gacgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAtt                                                                                                                <  1:3091311/150‑1 (MQ=255)
                                      gacgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAtt                                                                                                                <  2:361635/150‑1 (MQ=255)
                                        cgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTaa                                                                                                              <  1:2133026/150‑1 (MQ=255)
                                         gacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg                                                                                                             <  2:419479/150‑1 (MQ=255)
                                         gacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg                                                                                                             >  1:933080/1‑150 (MQ=255)
                                            aGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTa                                                                                                          >  1:3985330/1‑150 (MQ=255)
                                              cGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAATTAAc                                                                                                        >  1:1746628/1‑150 (MQ=255)
                                               gTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAcc                                                                                                       >  1:818052/1‑150 (MQ=255)
                                                 cACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCat                                                                                                     >  1:610385/1‑150 (MQ=255)
                                                       gTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCAt                                                                                               >  1:3302980/1‑150 (MQ=255)
                                                        ttttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATa                                                                                              <  1:246549/150‑1 (MQ=255)
                                                          ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACa                                                                                            <  1:1785367/150‑1 (MQ=255)
                                                            ataGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACaaa                                                                                          >  2:3562577/1‑150 (MQ=255)
                                                               gagaTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATa                                                                                       >  2:4000669/1‑150 (MQ=255)
                                                                 gaTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACa                                                                                     >  1:2995375/1‑150 (MQ=255)
                                                                  aTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACaa                                                                                    <  1:971965/150‑1 (MQ=255)
                                                                  aTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACaa                                                                                    >  1:3920853/1‑150 (MQ=255)
                                                                   tAGTTGAGATACCAGGTATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAAt                                                                                   <  1:2117945/150‑1 (MQ=255)
                                                                      ttGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGt                                                                                <  1:3687807/150‑1 (MQ=255)
                                                                        gagaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGttt                                                                              >  2:2102039/1‑150 (MQ=255)
                                                                         agaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTa                                                                             <  2:2027601/150‑1 (MQ=255)
                                                                          gaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTAt                                                                            <  2:2720092/150‑1 (MQ=255)
                                                                          gaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTAt                                                                            >  2:623336/1‑150 (MQ=255)
                                                                           aTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATg                                                                           >  2:1467529/1‑150 (MQ=255)
                                                                            tACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGt                                                                          <  2:396995/150‑1 (MQ=255)
                                                                             aCCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGtt                                                                         >  1:2949805/1‑150 (MQ=255)
                                                                              ccAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGAGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttc                                                                        >  1:1076954/1‑150 (MQ=255)
                                                                              ccAGGGAAGGTGTTAAAAACTACATTATTCTTACTATATATCAAAAATTAATAACCATTTAATAATTGATGTTATTGTAAGTGTGATTCTAAATTATCATAATAAACATAAAGTTAACCATATCCATACAAAATACAATGGTTTATGttc                                                                        >  1:1887011/1‑150 (MQ=255)
                                                                               cAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttct                                                                       >  2:1044795/1‑150 (MQ=255)
                                                                                aGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                                      >  2:1937170/1‑150 (MQ=255)
                                                                                aGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                                      <  1:2485079/150‑1 (MQ=255)
                                                                                aGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                                      >  1:103170/1‑150 (MQ=255)
                                                                                  ggATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATAAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGATAATGttcttaa                                                                    >  1:3048014/1‑148 (MQ=255)
                                                                                     tGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACAATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCaaaa                                                                 <  1:3104480/150‑1 (MQ=255)
                                                                                      ggTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaat                                                                >  2:2804065/1‑150 (MQ=255)
                                                                                      ggTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaat                                                                >  2:3047771/1‑150 (MQ=255)
                                                                                      ggTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaat                                                                >  2:549709/1‑150 (MQ=255)
                                                                                         gtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaataaa                                                             <  1:536106/150‑1 (MQ=255)
                                                                                          ttAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaataaat                                                            >  1:3582319/1‑150 (MQ=255)
                                                                                          ttAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaataaat                                                            >  1:261605/1‑150 (MQ=255)
                                                                                          ttAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaataaat                                                            <  2:1838619/150‑1 (MQ=255)
                                                                                           tAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAaataaata                                                           >  1:1167384/1‑150 (MQ=255)
                                                                                             aaaaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAAtaaataaa                                                         >  2:3399660/1‑150 (MQ=255)
                                                                                                aCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaa                                                      >  2:3351184/1‑150 (MQ=255)
                                                                                                aCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaa                                                      >  2:2768525/1‑150 (MQ=255)
                                                                                                   aCATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAt                                                   >  1:1786861/1‑150 (MQ=255)
                                                                                                    cATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtc                                                  <  1:2769991/150‑1 (MQ=255)
                                                                                                     attattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtca                                                 <  1:2320949/150‑1 (MQ=255)
                                                                                                      ttattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcat                                                <  2:1588414/150‑1 (MQ=255)
                                                                                                         ttCTTACGATATATCAAAAATTATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcattca                                             <  2:1621033/150‑1 (MQ=255)
                                                                                                          tCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcattcat                                            <  2:3843804/150‑1 (MQ=255)
                                                                                                           cTTAAGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATa                                           <  2:124207/150‑1 (MQ=255)
                                                                                                              aCGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAt                                        >  2:709079/1‑150 (MQ=255)
                                                                                                              aCGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAt                                        >  2:3189665/1‑150 (MQ=255)
                                                                                                               cGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAtt                                       >  1:1017101/1‑150 (MQ=255)
                                                                                                                gATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAttt                                      >  2:1053163/1‑150 (MQ=255)
                                                                                                                gATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAttt                                      <  1:73239/150‑1 (MQ=255)
                                                                                                                gATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAttt                                      >  2:1314351/1‑150 (MQ=255)
                                                                                                                 atatatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTa                                     <  1:1946555/150‑1 (MQ=255)
                                                                                                                 atatatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTa                                     >  1:2361377/1‑150 (MQ=255)
                                                                                                                 atatatCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTa                                     >  2:1156320/1‑150 (MQ=255)
                                                                                                                  tatatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTac                                    <  1:472076/150‑1 (MQ=255)
                                                                                                                   atatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTaca                                   >  2:3704036/1‑150 (MQ=255)
                                                                                                                   atatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTaca                                   <  2:222491/150‑1 (MQ=255)
                                                                                                                       cAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATc                               <  1:2233816/150‑1 (MQ=255)
                                                                                                                         aaaaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCAc                             >  2:2715227/1‑150 (MQ=255)
                                                                                                                           aaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACtt                           >  2:326983/1‑150 (MQ=255)
                                                                                                                               aTTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTaaaa                       <  1:4067098/150‑1 (MQ=255)
                                                                                                                               aTTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTaaaa                       <  1:1654922/150‑1 (MQ=255)
                                                                                                                                 tAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtc                     >  2:494148/1‑150 (MQ=255)
                                                                                                                                   aCCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtctc                   <  2:1382262/150‑1 (MQ=255)
                                                                                                                                   aCCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtctc                   >  2:1199359/1‑150 (MQ=255)
                                                                                                                                     cATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATAAAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCt                 >  1:3255339/1‑150 (MQ=255)
                                                                                                                                      aTTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTg                >  2:2904962/1‑150 (MQ=255)
                                                                                                                                      aTTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTg                <  1:1114289/150‑1 (MQ=255)
                                                                                                                                      aTTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTg                <  1:3914413/150‑1 (MQ=255)
                                                                                                                                        ttattaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGtt              >  1:2121267/1‑150 (MQ=255)
                                                                                                                                         tattaATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGttt             >  2:2461735/1‑150 (MQ=255)
                                                                                                                                          attaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTc            >  2:2259314/1‑150 (MQ=255)
                                                                                                                                            taATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCg          >  2:664653/1‑150 (MQ=255)
                                                                                                                                             aaTTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCGc         >  1:55731/1‑150 (MQ=255)
                                                                                                                                             aaTTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCGc         >  1:3779464/1‑150 (MQ=255)
                                                                                                                                                 gATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCGCACtt     <  2:325628/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCGCACTTTTT  >  CP009273/4531313‑4531610

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: