Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 147,946 G→T 31.5% I47I (ATC→ATA yadK ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273147,9460GT31.5% 507.0 / 275.1 304I47I (ATC→ATAyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (103/105);  new base T (42/54);  total (145/159)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 307 total aligned reads displayed.

TGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTC  >  CP009273/147797‑148095
                                                                                                                                                     |                                                                                                                                                     
 gATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGa                                                                                                                                                      >  2:110760/1‑150 (MQ=255)
     tAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTcc                                                                                                                                                  <  1:2051398/150‑1 (MQ=255)
        cacaTCGCCAGTTACCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCt                                                                                                                                               <  2:1212616/150‑1 (MQ=255)
          caTCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCttt                                                                                                                                             <  2:2809265/150‑1 (MQ=255)
           aTCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCtttt                                                                                                                                            >  1:3437834/1‑150 (MQ=255)
            tCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTg                                                                                                                                           <  1:787351/150‑1 (MQ=255)
                   ttCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAATATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGc                                                                                                                                    >  2:2293974/1‑150 (MQ=255)
                   ttCCCTGTAATGCACAGTCGGATACAGTGATAGTGAAATCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCAACCCAAAGACGCCAAGATAAAAATTACACAGGTTACTGATTCCGCTTTTGCAGGTGc                                                                                                                                    >  1:628683/1‑150 (MQ=255)
                    tCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCAGCTTTTGCAGGTGCt                                                                                                                                   >  2:3152149/1‑150 (MQ=255)
                       cTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGtt                                                                                                                                >  2:1552413/1‑150 (MQ=255)
                       cTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGtt                                                                                                                                >  1:293027/1‑150 (MQ=255)
                         gTAATGCACAGACGGATACAGTGAGAGTGAACTCGAGTCCACCTGGTTGAGAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACTCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGtttt                                                                                                                              <  2:2766736/150‑1 (MQ=255)
                             tGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTAc                                                                                                                          >  2:2294328/1‑150 (MQ=255)
                             tGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTAc                                                                                                                          >  2:246524/1‑150 (MQ=255)
                                  aGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTgg                                                                                                                     <  2:1543908/150‑1 (MQ=255)
                                    tCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAAGATCCCACCCCACCGACGCCAAGATCAATATGACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGca                                                                                                                   <  1:2065962/150‑1 (MQ=255)
                                     cGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGcaa                                                                                                                  >  1:3645769/1‑150 (MQ=255)
                                      ggATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGcaac                                                                                                                 >  1:4884033/1‑150 (MQ=255)
                                      ggATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGcaac                                                                                                                 <  1:1925466/150‑1 (MQ=255)
                                       gATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGcaaca                                                                                                                <  1:2903171/150‑1 (MQ=255)
                                        aTACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGcaacaa                                                                                                               >  1:2805180/1‑150 (MQ=255)
                                          aCAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAgg                                                                                                             <  1:1226094/150‑1 (MQ=255)
                                            aGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTc                                                                                                           <  2:3022828/150‑1 (MQ=255)
                                            aGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTGTTGCAGGTGCTGTTGTTTGAATTGGCACCAAGGTc                                                                                                           >  2:4497126/1‑150 (MQ=255)
                                             gTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCa                                                                                                          <  1:4650221/150‑1 (MQ=255)
                                               gATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAat                                                                                                        <  2:1644212/150‑1 (MQ=255)
                                                aTAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAata                                                                                                       <  1:4683639/150‑1 (MQ=255)
                                                   gTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATcc                                                                                                    <  1:238605/150‑1 (MQ=255)
                                                   gTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATcc                                                                                                    >  1:1796739/1‑150 (MQ=255)
                                                    tGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCt                                                                                                   >  1:4494443/1‑150 (MQ=255)
                                                         tctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCa                                                                                              <  1:5102920/150‑1 (MQ=255)
                                                         tctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCa                                                                                              <  1:2609523/150‑1 (MQ=255)
                                                         tctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCa                                                                                              <  1:77938/150‑1 (MQ=255)
                                                          ctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAg                                                                                             <  2:1467161/150‑1 (MQ=255)
                                                            tttCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCa                                                                                           >  2:627327/1‑150 (MQ=255)
                                                             ttCCACCTGGTATATAACTTTCAGGAGTAAAATTACCTTAAAAATATCCCACACCAACTACGCCAAGATCAATATTACCCTGTTTACTGATTCCGCTTTTTCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTTTCCAGCaa                                                                                          >  2:1538885/1‑150 (MQ=255)
                                                                cACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATgg                                                                                       >  2:2236856/1‑150 (MQ=255)
                                                                cACCTGGTTGATAACAATCAGGACTAACATTACCAGAAAAATATCCCACCACAACGACGCCAAGATCAATATTAACCTGGTTAATGATTCCGCTTTAGAAGGTGCTGTTATTTACATTGGAAACAATGTCAATATCCTGTCCAGCAATgg                                                                                       >  2:4277930/1‑150 (MQ=255)
                                                                  ccTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCt                                                                                     >  1:4810407/1‑150 (MQ=255)
                                                                   cTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                                    >  2:4157575/1‑150 (MQ=255)
                                                                    tGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTgg                                                                                   >  2:763769/1‑150 (MQ=255)
                                                                     ggTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                  >  2:1081252/1‑150 (MQ=255)
                                                                          aTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                             >  2:3408752/1‑150 (MQ=255)
                                                                          aTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                             <  2:178376/150‑1 (MQ=255)
                                                                          aTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                             <  2:1409282/150‑1 (MQ=255)
                                                                            aaCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGa                                                                           <  1:1473812/150‑1 (MQ=255)
                                                                             aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                          <  1:3577943/150‑1 (MQ=255)
                                                                              cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                         <  1:4602731/150‑1 (MQ=255)
                                                                              cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                         <  1:5160075/150‑1 (MQ=255)
                                                                               tttCAGGAGTAACATTACGTGACAAATATCCCACCCCAACGACGCCACGATCACTATTACCCTGGTGACTTAGTCCGCTGTTGCGGGTGCGGGGGTTTACCTTGGCCAGCAGGTCCAGCTCCGGTCCAGCTCTGGCGGGCCAgcgggtgc                                                                        >  2:4840681/1‑143 (MQ=255)
                                                                                ttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGTACAGGCGATGCa                                                                       >  1:3489526/1‑150 (MQ=255)
                                                                                     gAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACa                                                                  <  2:3939097/150‑1 (MQ=255)
                                                                                       gTAACATTACCTGAAAAATCTCCCACCCCAACGACGCCAAGATCAATCTTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAACCAGt                                                                <  2:4778902/150‑1 (MQ=255)
                                                                                       gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGt                                                                <  2:2537560/150‑1 (MQ=255)
                                                                                       gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTAACCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTATTACATTGGCAACAAGGTCAAAATCCTGTCCAGCAATGGCTGGACAGGCCATGCAAAACAGt                                                                >  2:689981/1‑150 (MQ=255)
                                                                                         aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                              >  2:2673394/1‑150 (MQ=255)
                                                                                         aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                              <  1:4270191/150‑1 (MQ=255)
                                                                                         aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                              >  1:2776337/1‑150 (MQ=255)
                                                                                          aCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAAt                                                             <  2:4993165/150‑1 (MQ=255)
                                                                                           cATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATg                                                            >  2:488495/1‑150 (MQ=255)
                                                                                           cATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATg                                                            <  1:586653/150‑1 (MQ=255)
                                                                                            aTTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGa                                                           >  1:3815431/1‑150 (MQ=255)
                                                                                            aTTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGa                                                           <  1:33723/150‑1 (MQ=255)
                                                                                             ttACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAACATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                          <  2:781939/150‑1 (MQ=255)
                                                                                              tACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                         >  1:3043544/1‑150 (MQ=255)
                                                                                                 cTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaa                                                      <  2:3697530/150‑1 (MQ=255)
                                                                                                  tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                     <  1:151711/150‑1 (MQ=255)
                                                                                                   gAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaac                                                    >  1:3317563/1‑150 (MQ=255)
                                                                                                   gAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaac                                                    <  1:3185994/150‑1 (MQ=255)
                                                                                                    aaaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                                   <  2:4601029/150‑1 (MQ=255)
                                                                                                         tatCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatt                                              >  2:1093463/1‑150 (MQ=255)
                                                                                                          atCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                             <  1:4355786/150‑1 (MQ=255)
                                                                                                          atCCCACCCCAACGACGCCAAGATCAATAGTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                             <  2:2985804/150‑1 (MQ=255)
                                                                                                           tCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattat                                            <  2:1895601/150‑1 (MQ=255)
                                                                                                                ccccAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTATTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                       <  2:2079490/150‑1 (MQ=255)
                                                                                                                ccccAACGACGCCAAGATCAATATTACCCTGGTTACTTAGTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                       <  2:4682150/150‑1 (MQ=255)
                                                                                                                 cccAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                                      <  1:842594/150‑1 (MQ=255)
                                                                                                                 cccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                                      <  2:5154297/150‑1 (MQ=255)
                                                                                                                 cccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                                      <  1:985545/150‑1 (MQ=255)
                                                                                                                   cAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                    <  2:157767/150‑1 (MQ=255)
                                                                                                                    aaCGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcat                                   >  1:4306754/1‑150 (MQ=255)
                                                                                                                       gacgCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAg                                >  1:5190272/1‑150 (MQ=255)
                                                                                                                           ccAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTAATCAGCTTa                            >  2:1180621/1‑150 (MQ=255)
                                                                                                                            cAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAc                           >  1:1138144/1‑150 (MQ=255)
                                                                                                                            cAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAc                           >  1:3755245/1‑150 (MQ=255)
                                                                                                                            cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAc                           >  1:49622/1‑150 (MQ=255)
                                                                                                                             aaGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACg                          <  1:2198724/150‑1 (MQ=255)
                                                                                                                                   aaTATTACCCTGGTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                    <  1:3985533/150‑1 (MQ=255)
                                                                                                                                   aaTATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                    >  1:4110584/1‑150 (MQ=255)
                                                                                                                                    atatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                   >  1:3677961/1‑150 (MQ=255)
                                                                                                                                    agatTACCCTGGTTACTGAGTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                   <  2:4818305/148‑1 (MQ=255)
                                                                                                                                         aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATg              <  2:4409862/150‑1 (MQ=255)
                                                                                                                                             tGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCata          >  2:3832242/1‑150 (MQ=255)
                                                                                                                                              ggTTACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCATCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatat         <  2:3532319/150‑1 (MQ=255)
                                                                                                                                              ggTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACTAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTATCTTCATCAGCTTACGTTGAGTGGGATGCatat         >  2:1983167/1‑150 (MQ=255)
                                                                                                                                              ggTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCTTCAGCTTACGTTGAGTGGGATGCatat         >  2:5243874/1‑150 (MQ=255)
                                                                                                                                              ggTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatat         >  1:5089829/1‑150 (MQ=255)
                                                                                                                                               tttACTGATTCAGCTTTTTCAGGTGCTGTTTTTTACCTTGGAAACAAGGTCAATATCCTGTCCTGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatata        <  2:4561247/149‑1 (MQ=255)
                                                                                                                                                ttACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatatat       <  1:3168540/150‑1 (MQ=255)
                                                                                                                                                 tACTTATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCACAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCAtatata      >  2:1344608/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTC  >  CP009273/147797‑148095

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