Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,433,219 G→T 7.9% D341Y (GAT→TAT)  tamA → translocation and assembly module for autotransporter export, outer membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,433,2190GT 7.9% 1061.5 / 39.9 286D341Y (GAT→TAT) tamAtranslocation and assembly module for autotransporter export, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (151/112);  new base T (20/3);  total (171/115)
Fisher's exact test for biased strand distribution p-value = 6.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

GAACTCTTATGGTCACAGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGTTTTATCCTGGGGTG  >  CP009273/4433071‑4433368
                                                                                                                                                    |                                                                                                                                                     
  aCTCTTATGGTCACAGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATa                                                                                                                                                    >  1:3221246/1‑150 (MQ=255)
      ttATGGTCACAGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGa                                                                                                                                                <  2:3353842/150‑1 (MQ=255)
       tATGGTCACAGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGaa                                                                                                                                               >  1:1262602/1‑150 (MQ=255)
         tGGTCACAGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATc                                                                                                                                             >  1:77817/1‑150 (MQ=255)
               cattCTGACCACCAGAACTAGTATTACCGCGCCGGAACACACCCTCGACTTCAGCTATAAAATGACGCTGATGAAGAATCCACTGGAACAATATTATTTCGACCAGGGCGTTTTTAATCGCACTGAACTCAACGAAACCGAATCTGACTc                                                                                                                                       >  1:1930728/4‑150 (MQ=255)
               caGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTc                                                                                                                                       >  1:287715/1‑150 (MQ=255)
               caGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTc                                                                                                                                       >  1:288305/1‑150 (MQ=255)
               caGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTc                                                                                                                                       >  1:3529896/1‑150 (MQ=255)
                 gTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCa                                                                                                                                     >  1:3872652/1‑150 (MQ=255)
                 gTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCa                                                                                                                                     >  1:3872646/1‑150 (MQ=255)
                  tCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCAc                                                                                                                                    >  1:64327/1‑150 (MQ=255)
                   cTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACt                                                                                                                                   <  2:3743530/150‑1 (MQ=255)
                   cTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACt                                                                                                                                   <  2:1479921/150‑1 (MQ=255)
                   cTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACt                                                                                                                                   >  1:3698492/1‑150 (MQ=255)
                     gACCAACAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACttc                                                                                                                                 <  2:1274000/150‑3 (MQ=255)
                        caccaGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGACTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACTATACCGAATCTGACTCCACTACGCt                                                                                                                              >  1:916500/1‑150 (MQ=255)
                        caccaGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACTATACCGAATCTGACTCCACTACGCt                                                                                                                              >  1:1037945/1‑150 (MQ=255)
                         accaGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACTATACCGAATCTGACTCCACTACGCtg                                                                                                                             >  2:1023888/1‑150 (MQ=255)
                         accaGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCtg                                                                                                                             >  2:130453/1‑150 (MQ=255)
                            aGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCtggtg                                                                                                                          >  1:1145340/1‑150 (MQ=255)
                              tACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGc                                                                                                                        >  1:2547430/1‑150 (MQ=255)
                               aCTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCt                                                                                                                       >  1:27284/1‑150 (MQ=255)
                                cTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCtt                                                                                                                      <  1:1896786/150‑1 (MQ=255)
                                cTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCAACTGGAACAAGATTATTTGGTGCAGGGCGGTTTTAAGCGAACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCtt                                                                                                                      >  2:3715589/1‑150 (MQ=255)
                                     aTTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACTATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGc                                                                                                                 >  1:1506059/1‑150 (MQ=255)
                                     aTTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGc                                                                                                                 >  2:565130/1‑150 (MQ=255)
                                     aTTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCAACTACGCTGGTGGCTTCTCGc                                                                                                                 >  2:4211296/1‑150 (MQ=255)
                                       ttCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTa                                                                                                               <  1:985293/150‑1 (MQ=255)
                                         ccGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACt                                                                                                             <  2:1558719/150‑1 (MQ=255)
                                             gcCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGa                                                                                                         >  1:1389914/1‑150 (MQ=255)
                                               cGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGAtc                                                                                                       >  2:2098986/1‑150 (MQ=255)
                                                ggAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGAtct                                                                                                      >  1:3287387/1‑150 (MQ=255)
                                                ggAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGAtct                                                                                                      <  1:3615369/150‑1 (MQ=255)
                                                  aaCAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGAtctct                                                                                                    <  1:1180272/150‑1 (MQ=255)
                                                   aCAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGAtctctc                                                                                                   <  1:796218/150‑1 (MQ=255)
                                                    cAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTcc                                                                                                  <  1:3944589/150‑1 (MQ=255)
                                                    cAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTcc                                                                                                  >  1:1968576/1‑150 (MQ=255)
                                                    cAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTcc                                                                                                  >  1:1968711/1‑150 (MQ=255)
                                                    cAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTcc                                                                                                  >  1:1973309/1‑150 (MQ=255)
                                                       aCCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGc                                                                                               <  1:1781261/150‑1 (MQ=255)
                                                            cGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTg                                                                                          >  2:1277269/1‑150 (MQ=255)
                                                            cGACTTCAGCTATAAAATGCCGCTGCGGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTg                                                                                          >  2:1083744/1‑150 (MQ=255)
                                                                 tCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGc                                                                                     >  2:2279954/1‑150 (MQ=255)
                                                                    gCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTg                                                                                  <  2:1796130/150‑1 (MQ=255)
                                                                      tataAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCAGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGcc                                                                                >  2:4480821/1‑150 (MQ=255)
                                                                        taAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCAt                                                                              >  1:2143049/1‑150 (MQ=255)
                                                                         aaaaTGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCAtt                                                                             >  2:3140616/1‑150 (MQ=255)
                                                                          aaaTGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTa                                                                            >  1:1589043/1‑150 (MQ=255)
                                                                             tGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACTATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAAcc                                                                         <  1:3191820/150‑1 (MQ=255)
                                                                             tGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAAcc                                                                         >  2:460491/1‑150 (MQ=255)
                                                                               ccGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTg                                                                       >  2:2378368/1‑150 (MQ=255)
                                                                               ccGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTg                                                                       >  2:3952257/1‑150 (MQ=255)
                                                                               ccGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTg                                                                       >  2:2512258/1‑150 (MQ=255)
                                                                                cGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTgc                                                                      <  2:2105440/150‑1 (MQ=255)
                                                                                  ctgctgAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTgcgc                                                                    >  1:53490/1‑150 (MQ=255)
                                                                                    gctgAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTg                                                                  >  2:607100/1‑150 (MQ=255)
                                                                                      tgAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGa                                                                >  2:1926327/1‑150 (MQ=255)
                                                                                        aagaaTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACTATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGt                                                              >  1:4385047/1‑150 (MQ=255)
                                                                                        aagaaTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGt                                                              <  2:3819900/150‑1 (MQ=255)
                                                                                         agaaTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGtc                                                             <  2:3596334/150‑1 (MQ=255)
                                                                                           aaTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACTATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGtctc                                                           >  1:1164704/1‑150 (MQ=255)
                                                                                           aaTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGtctc                                                           >  2:1334140/1‑150 (MQ=255)
                                                                                           aaTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGtctc                                                           <  2:304428/150‑1 (MQ=255)
                                                                                           aaTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGtctc                                                           <  1:3356944/150‑1 (MQ=255)
                                                                                            aTCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCg                                                          >  1:2952589/1‑150 (MQ=255)
                                                                                               cACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGAcc                                                       >  1:2585829/1‑150 (MQ=255)
                                                                                               cACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGAcc                                                       <  2:2047363/150‑1 (MQ=255)
                                                                                                aCTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCa                                                      <  1:2197792/150‑1 (MQ=255)
                                                                                                  tGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACt                                                    >  2:4249714/1‑150 (MQ=255)
                                                                                                   ggAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACtt                                                   <  2:2771568/150‑1 (MQ=255)
                                                                                                     aaCAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTa                                                 <  2:4112025/150‑1 (MQ=255)
                                                                                                      aCAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTAc                                                >  1:1100860/1‑150 (MQ=255)
                                                                                                      aCAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTAc                                                >  1:1459952/1‑150 (MQ=255)
                                                                                                           attattTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAgg                                           <  1:3116672/150‑1 (MQ=255)
                                                                                                            ttattTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAggg                                          >  2:2837479/1‑150 (MQ=255)
                                                                                                              attTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTg                                        >  1:4064942/1‑150 (MQ=255)
                                                                                                                 tGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGaaa                                     >  1:235433/1‑150 (MQ=255)
                                                                                                                  ggTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAAt                                    <  1:505477/150‑1 (MQ=255)
                                                                                                                  ggTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAAt                                    <  2:3872646/150‑1 (MQ=255)
                                                                                                                   gTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAAtt                                   <  1:3242893/150‑1 (MQ=255)
                                                                                                                       aGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCa                               >  1:3674187/1‑150 (MQ=255)
                                                                                                                        gggCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCaa                              >  2:4270607/1‑150 (MQ=255)
                                                                                                                        gggCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCaa                              >  2:2645522/1‑150 (MQ=255)
                                                                                                                          gCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATa                            <  2:3193569/150‑1 (MQ=255)
                                                                                                                             gTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCa                         >  1:743354/1‑150 (MQ=255)
                                                                                                                                 tAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTAGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGAt                     >  1:4538141/1‑150 (MQ=255)
                                                                                                                                  aaGCGCACTGACCTGAACTATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATg                    <  2:1506059/150‑1 (MQ=255)
                                                                                                                                  aaGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATg                    <  2:449477/150‑1 (MQ=255)
                                                                                                                                   agagcACTGACCTGAACGATACAGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATAACACGATGc                   >  2:1815180/4‑150 (MQ=255)
                                                                                                                                     cgcACTGACCTGAACTATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTg                 >  1:3627336/1‑150 (MQ=255)
                                                                                                                                     cgcACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTg                 >  2:3230067/1‑150 (MQ=255)
                                                                                                                                     cgcACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTg                 <  1:1288571/150‑1 (MQ=255)
                                                                                                                                      gcACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGt                <  2:373485/150‑1 (MQ=255)
                                                                                                                                       cACTGACCTGAACGATACCGAATCTGCCTCCACTCCGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGtt               <  2:2406253/150‑1 (MQ=255)
                                                                                                                                               tGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGTTTTATCCTg       >  2:77196/1‑150 (MQ=255)
                                                                                                                                               tGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGTTTTATCCTg       <  1:4316199/150‑1 (MQ=255)
                                                                                                                                               tGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGTTTTATCCTg       >  1:595145/1‑150 (MQ=255)
                                                                                                                                               tGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGTTTTATCCTg       >  2:1772620/1‑150 (MQ=255)
                                                                                                                                                 aaCGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGTTTTATCCTggg     >  1:3086651/1‑150 (MQ=255)
                                                                                                                                                  aCGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGTTTTATCCTgggg    <  2:1887052/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAACTCTTATGGTCACAGTCTGACCACCAGTACTAGTATTTCCGCGCCGGAACAGACCCTCGACTTCAGCTATAAAATGCCGCTGCTGAAGAATCCACTGGAACAATATTATTTGGTGCAGGGCGGTTTTAAGCGCACTGACCTGAACGATACCGAATCTGACTCCACTACGCTGGTGGCTTCTCGCTACTGGGATCTCTCCAGCGGCTGGCAGCGTGCCATTAACCTGCGCTGGAGTCTCGACCACTTTACTCAGGGTGAAATTACCAATACCACGATGCTGTTTTATCCTGGGGTG  >  CP009273/4433071‑4433368

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: