Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,529,741 C→T 6.2% intergenic (‑423/‑1033) nanC ← / → fimB N‑acetylnuraminic acid outer membrane channel protein/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,529,7410CT 6.2% 606.2 / 19.0 176intergenic (‑423/‑1033)nanC/fimBN‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base C (79/86);  new base T (4/7);  total (83/93)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 178 total aligned reads displayed.

ACCCAGAGTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTA  >  CP009273/4529593‑4529890
                                                                                                                                                    |                                                                                                                                                     
aCCCAGAGTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCTATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCAcc                                                                                                                                                      <  2:2344109/150‑1 (MQ=255)
 cccAGAGTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCt                                                                                                                                                     >  1:2009325/1‑150 (MQ=255)
  ccAGAGTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCtt                                                                                                                                                    <  2:1062728/150‑1 (MQ=255)
  ccAGAGTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCtt                                                                                                                                                    <  1:1653232/150‑1 (MQ=255)
    agagTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTACCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTta                                                                                                                                                  >  1:3937823/1‑150 (MQ=255)
      agTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTtata                                                                                                                                                <  2:2111793/150‑1 (MQ=255)
       gTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCATCTTTATAc                                                                                                                                               <  1:2248836/150‑1 (MQ=255)
        tCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCATCTTTATAcc                                                                                                                                              <  2:3588600/150‑1 (MQ=255)
           tAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGt                                                                                                                                           <  1:357368/150‑1 (MQ=255)
           tAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGACATATATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGt                                                                                                                                           >  2:4098748/1‑150 (MQ=255)
           tAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGAGATGTTTGACATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGTTGCACCTTTATACCTGt                                                                                                                                           >  1:3669842/1‑150 (MQ=255)
            aaaaTAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTCTCTTGGTGCACCTTTATACCTGtt                                                                                                                                          >  2:4892143/1‑150 (MQ=255)
            aaaaTAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGtt                                                                                                                                          >  1:1857979/1‑150 (MQ=255)
            aaaaTAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGtt                                                                                                                                          <  2:4586086/150‑1 (MQ=255)
            aaaaTAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGtt                                                                                                                                          >  1:3109385/1‑150 (MQ=255)
            aaaaTAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGtt                                                                                                                                          <  2:4586468/150‑1 (MQ=255)
                tAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATAc                                                                                                                                      <  2:4992970/150‑1 (MQ=255)
                 aGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCATCTTTATACCTGTTATAcc                                                                                                                                     >  1:1040025/1‑150 (MQ=255)
                       tCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCa                                                                                                                               >  2:620565/1‑150 (MQ=255)
                       tCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTACATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCa                                                                                                                               <  2:624997/150‑1 (MQ=255)
                          gTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTGTACCAGATCaaaa                                                                                                                            <  1:3657998/150‑1 (MQ=255)
                          gTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCaaaa                                                                                                                            <  1:3657872/150‑1 (MQ=255)
                             gATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATc                                                                                                                         >  2:1252826/1‑150 (MQ=255)
                             gATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATc                                                                                                                         >  2:4850655/1‑150 (MQ=255)
                             gATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATc                                                                                                                         >  2:3834638/1‑150 (MQ=255)
                               tctcATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCAc                                                                                                                       <  1:3518158/150‑1 (MQ=255)
                               tctcATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCAc                                                                                                                       <  1:2235088/150‑1 (MQ=255)
                               tctcATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCAc                                                                                                                       <  1:2235087/150‑1 (MQ=255)
                                ctcATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACg                                                                                                                      >  1:4589293/1‑150 (MQ=255)
                                   aTTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCaa                                                                                                                   <  1:3869921/150‑1 (MQ=255)
                                   aTTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCaa                                                                                                                   <  1:3869962/150‑1 (MQ=255)
                                   aTTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCaa                                                                                                                   >  1:1605008/1‑150 (MQ=255)
                                       tCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCa                                                                                                               <  2:163733/150‑1 (MQ=255)
                                           atTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATaca                                                                                                           >  2:3785373/1‑150 (MQ=255)
                                             tCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATacaac                                                                                                         >  1:60097/1‑150 (MQ=255)
                                              ccTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATacaaca                                                                                                        >  2:5131223/1‑150 (MQ=255)
                                               cTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATAcaacaa                                                                                                       >  2:1664380/1‑150 (MQ=255)
                                                tttCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACaaa                                                                                                      <  2:3451495/150‑1 (MQ=255)
                                                 ttCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACaaaa                                                                                                     <  1:454691/150‑1 (MQ=255)
                                                     gACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCATCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAg                                                                                                 <  2:809373/150‑1 (MQ=255)
                                                     gACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAg                                                                                                 <  2:2057034/150‑1 (MQ=255)
                                                     gACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAg                                                                                                 >  2:4569404/1‑150 (MQ=255)
                                                        tGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGAtt                                                                                              <  1:2782635/150‑1 (MQ=255)
                                                          aaGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTg                                                                                            >  2:2131625/1‑150 (MQ=255)
                                                          aaGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTg                                                                                            <  2:3864872/150‑1 (MQ=255)
                                                           aGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGc                                                                                           <  2:1197180/150‑1 (MQ=255)
                                                              gCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAAt                                                                                        <  1:2660725/150‑1 (MQ=255)
                                                                    tcatcaCCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCgtgt                                                                                  <  1:2464242/150‑1 (MQ=255)
                                                                     catcaCCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATGGACTATCTTGGTGCACCTGTATACCTGTTATACCAGATCAAACATCACGCAATCCATACAACAACACCAGATTTGCAAGTCGTGTc                                                                                 >  2:1816313/1‑150 (MQ=255)
                                                                       tcaCCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTcac                                                                               >  1:1006775/1‑150 (MQ=255)
                                                                          cccGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACaaa                                                                            >  2:1254792/1‑150 (MQ=255)
                                                                          cccGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACaaa                                                                            <  2:383498/150‑1 (MQ=255)
                                                                            cGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAat                                                                          >  1:3556162/1‑150 (MQ=255)
                                                                            cGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAat                                                                          >  2:3391141/1‑150 (MQ=255)
                                                                               tGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATg                                                                       >  2:1322406/1‑150 (MQ=255)
                                                                                gagaAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGt                                                                      >  2:744438/1‑150 (MQ=255)
                                                                                 agaaagaaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTc                                                                     >  1:850218/1‑150 (MQ=255)
                                                                                 agaaagaaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTc                                                                     <  1:3463401/150‑1 (MQ=255)
                                                                                 agaaagaaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTc                                                                     >  1:2836299/1‑150 (MQ=255)
                                                                                   aaagaaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGa                                                                   <  1:3695772/150‑1 (MQ=255)
                                                                                    aagaaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCACTTCGTGTCACACAATATGTCGAt                                                                  >  1:138103/1‑150 (MQ=255)
                                                                                     agaaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATc                                                                 <  2:3166021/150‑1 (MQ=255)
                                                                                     agaaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATc                                                                 <  2:3165938/150‑1 (MQ=255)
                                                                                      gaaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCt                                                                <  2:2779944/150‑1 (MQ=255)
                                                                                       aaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCtt                                                               >  1:4172937/1‑150 (MQ=255)
                                                                                       aaCTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCtt                                                               >  1:3604175/1‑150 (MQ=255)
                                                                                           gAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTc                                                           <  2:1605008/150‑1 (MQ=255)
                                                                                            aaGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCt                                                          <  2:719528/150‑1 (MQ=255)
                                                                                             aGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTa                                                         <  1:689627/150‑1 (MQ=255)
                                                                                              gTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTaa                                                        >  1:5149028/1‑150 (MQ=255)
                                                                                                    tttGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGaa                                                  >  1:1826298/1‑150 (MQ=255)
                                                                                                     ttGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAg                                                 >  2:3920907/1‑150 (MQ=255)
                                                                                                      tGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGa                                                >  2:2561668/1‑150 (MQ=255)
                                                                                                      tGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGa                                                >  2:2561476/1‑150 (MQ=255)
                                                                                                       gCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGAt                                               <  2:3397574/150‑1 (MQ=255)
                                                                                                        ccATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGGTCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATg                                              <  2:1859581/150‑1 (MQ=255)
                                                                                                           tAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCa                                           <  1:1738449/150‑1 (MQ=255)
                                                                                                            aGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCAt                                          <  2:209472/150‑1 (MQ=255)
                                                                                                             gATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCAtg                                         <  2:3297724/150‑1 (MQ=255)
                                                                                                               ttGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCAtgtg                                       <  1:1016436/150‑1 (MQ=255)
                                                                                                                   gagaTGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGc                                   <  1:1281930/150‑1 (MQ=255)
                                                                                                                   gagaTGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTAACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGc                                   >  2:2259439/1‑150 (MQ=255)
                                                                                                                    agaTGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGcc                                  >  2:4527809/1‑150 (MQ=255)
                                                                                                                          ccTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGa                            <  1:1058317/150‑1 (MQ=255)
                                                                                                                           cTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGaa                           >  2:4697615/1‑150 (MQ=255)
                                                                                                                            tGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAac                          >  2:2737355/1‑150 (MQ=255)
                                                                                                                               ctcATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACt                       <  2:418717/150‑1 (MQ=255)
                                                                                                                               ctcATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACt                       <  1:482952/150‑1 (MQ=255)
                                                                                                                                  aTTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGc                    <  1:2161561/150‑1 (MQ=255)
                                                                                                                                  aTTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGc                    <  2:1445111/150‑1 (MQ=255)
                                                                                                                                   ttGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCg                   <  2:4789595/150‑1 (MQ=255)
                                                                                                                                    tGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCgg                  >  1:913104/1‑150 (MQ=255)
                                                                                                                                      aCTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCggtg                >  1:3154098/1‑150 (MQ=255)
                                                                                                                                      aCTATCTTGGTGCACCTGTATACCTGTTAGACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCggtg                <  2:4777157/150‑1 (MQ=255)
                                                                                                                                        tATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCggtggt              <  2:171637/150‑1 (MQ=255)
                                                                                                                                           ctgggTGCACCTTTCGACCTGTTAGACCAGATCAAAAATCACGCAAGCCATACAACAAAACCAGATTTGCAATTCGTGTCACCAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCtt           <  2:138103/147‑1 (MQ=255)
                                                                                                                                           cTTGGTGCATCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCtt           <  2:365587/150‑1 (MQ=255)
                                                                                                                                            ttGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTc          <  2:4417973/150‑1 (MQ=255)
                                                                                                                                              ggTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCaa        <  2:1298051/150‑1 (MQ=255)
                                                                                                                                                   aCCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTATTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACt   >  2:3070763/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCCAGAGTCATAAAATAGAATGTCTGTTGATCTCATTTTCATATTCCTTTCAGACTGAAGTGCGTTTTCATCACCCGATGAGAAAGAACTGAAGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTA  >  CP009273/4529593‑4529890

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: