Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,057,706 A→T 7.4% intergenic (‑51/+17) ompL ← / ← yihO outer membrane porin L/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,057,7060AT 7.4% 937.3 / 50.8 285intergenic (‑51/+17)ompL/yihOouter membrane porin L/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (153/111);  new base T (8/13);  total (161/124)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 291 total aligned reads displayed.

GGCAAGATTATATGCCTCCCTGTTTTCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAAT  >  CP009273/4057557‑4057854
                                                                                                                                                     |                                                                                                                                                    
     gATTATATGCCTCCCTGTTTTCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCg                                                                                                                                                 <  2:3835768/150‑1 (MQ=255)
        tataTGCCTCCCTGTTTTCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCt                                                                                                                                              <  1:2856484/150‑1 (MQ=255)
        tataTGCCTCCCTGTTTTCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCt                                                                                                                                              <  1:4245021/150‑1 (MQ=255)
             gCCTCCCTGTTTTCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGa                                                                                                                                         >  2:4162168/1‑150 (MQ=255)
                tCCCTGTTTTCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACAt                                                                                                                                      >  1:4975877/1‑150 (MQ=255)
                   cTGTTTTCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTta                                                                                                                                   >  2:42981/1‑150 (MQ=255)
                       tttCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAAc                                                                                                                               <  1:2658800/150‑1 (MQ=255)
                        ttCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACg                                                                                                                              <  1:3322502/150‑1 (MQ=255)
                        ttCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACg                                                                                                                              >  1:4428230/1‑150 (MQ=255)
                          cGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGtt                                                                                                                            >  1:4079248/1‑150 (MQ=255)
                           gACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTa                                                                                                                           >  1:794159/1‑150 (MQ=255)
                              gTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACgg                                                                                                                        >  1:4996853/1‑150 (MQ=255)
                              gTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACgg                                                                                                                        <  1:669589/150‑1 (MQ=255)
                               tATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGa                                                                                                                       <  1:2617740/150‑1 (MQ=255)
                                   cgcCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGcc                                                                                                                   <  2:3107845/150‑1 (MQ=255)
                                   cgcCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGcc                                                                                                                   <  2:170927/150‑1 (MQ=255)
                                   cgcCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGcc                                                                                                                   <  1:57254/150‑1 (MQ=255)
                                   cgcCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGcc                                                                                                                   <  2:3380576/150‑1 (MQ=255)
                                    gcCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGAAGCCg                                                                                                                  <  2:1178865/150‑1 (MQ=255)
                                    gcCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCg                                                                                                                  >  2:3240816/1‑150 (MQ=255)
                                      cAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGtt                                                                                                                <  2:2586925/150‑1 (MQ=255)
                                      cAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGtt                                                                                                                >  1:1447963/1‑150 (MQ=255)
                                       aGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTTCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGttt                                                                                                               <  1:1940016/150‑1 (MQ=255)
                                        gCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGtttt                                                                                                              <  1:107066/150‑1 (MQ=255)
                                         cAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGttttt                                                                                                             >  2:3053629/1‑150 (MQ=255)
                                          aaaTACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGtttttt                                                                                                            >  1:1392426/1‑150 (MQ=255)
                                          aaaTACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGtttttt                                                                                                            >  1:5040543/1‑150 (MQ=255)
                                             tACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTttctt                                                                                                         >  2:3681122/1‑150 (MQ=255)
                                               ccGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTtcttct                                                                                                       <  1:1717620/150‑1 (MQ=255)
                                                 gAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGt                                                                                                     <  1:4189633/150‑1 (MQ=255)
                                                 gAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGt                                                                                                     >  2:3421615/1‑150 (MQ=255)
                                                  aGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTc                                                                                                    <  1:2633251/150‑1 (MQ=255)
                                                   ggCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCg                                                                                                   <  1:3184170/150‑1 (MQ=255)
                                                    gCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGc                                                                                                  >  1:3894574/1‑150 (MQ=255)
                                                      agagGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCa                                                                                                >  2:64562/1‑150 (MQ=255)
                                                             agagAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCttt                                                                                         >  2:1262445/1‑150 (MQ=255)
                                                                 aagaTAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTAttt                                                                                     >  2:2120923/1‑150 (MQ=255)
                                                                  agaTAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTAtttt                                                                                    <  1:4159954/150‑1 (MQ=255)
                                                                  agaTAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTAtttt                                                                                    <  2:2056106/150‑1 (MQ=255)
                                                                   gaTAATAAAATTATTGCATTAATCTTTTTCATAATAATTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGAGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTAttttt                                                                                   <  2:4729938/150‑1 (MQ=255)
                                                                     taataaAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGc                                                                                 >  2:2361430/1‑150 (MQ=255)
                                                                      aataaAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCg                                                                                >  2:2449901/1‑150 (MQ=255)
                                                                         aaaaTTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttg                                                                             <  1:3763530/150‑1 (MQ=255)
                                                                         aaaaTTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttg                                                                             <  1:1854748/150‑1 (MQ=255)
                                                                           aaTTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggt                                                                           <  1:3847163/150‑1 (MQ=255)
                                                                            attattGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggtt                                                                          <  2:1341118/150‑1 (MQ=255)
                                                                             ttattGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggttg                                                                         >  1:2286902/1‑150 (MQ=255)
                                                                              tattGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGa                                                                        <  2:3859034/150‑1 (MQ=255)
                                                                               attGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAt                                                                       <  1:1151051/150‑1 (MQ=255)
                                                                                 tGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTTGTTGAttt                                                                     >  1:4022690/1‑150 (MQ=255)
                                                                                  gCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTc                                                                    <  1:4277855/150‑1 (MQ=255)
                                                                                      tAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                                                >  2:5031225/1‑150 (MQ=255)
                                                                                        aTCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaat                                                              >  2:4522748/1‑150 (MQ=255)
                                                                                         tCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaata                                                             <  2:4496007/150‑1 (MQ=255)
                                                                                              ttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCTCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCaa                                                        >  2:2012732/1‑150 (MQ=255)
                                                                                              ttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCaa                                                        >  2:3114905/1‑150 (MQ=255)
                                                                                               tCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAAt                                                       >  1:3033503/1‑150 (MQ=255)
                                                                                               tCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAAt                                                       <  2:2293511/150‑1 (MQ=255)
                                                                                                  tcataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGGCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCa                                                    >  2:3061774/3‑150 (MQ=255)
                                                                                                  taataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCa                                                    >  2:2794985/1‑150 (MQ=255)
                                                                                                   aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaa                                                   >  1:4916247/1‑150 (MQ=255)
                                                                                                   aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaa                                                   >  2:4534333/1‑150 (MQ=255)
                                                                                                    ataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaaa                                                  >  2:4568520/1‑150 (MQ=255)
                                                                                                    ataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaaa                                                  >  1:3523609/1‑150 (MQ=255)
                                                                                                    ataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaaa                                                  >  2:2871435/1‑150 (MQ=255)
                                                                                                    ataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaaa                                                  >  1:1507907/1‑150 (MQ=255)
                                                                                                     taaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAATTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAc                                                 <  1:3039130/150‑1 (MQ=255)
                                                                                                      aaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                                                <  1:4361212/150‑1 (MQ=255)
                                                                                                       aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGAGTTGGTTGATTTCCTCAATAATCAATGCAAACCg                                               >  1:3010607/1‑150 (MQ=255)
                                                                                                           cATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                                           >  1:819016/1‑150 (MQ=255)
                                                                                                            aTTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAACGttttt                                          >  2:840724/1‑150 (MQ=255)
                                                                                                             ttttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc                                         >  1:2910255/1‑150 (MQ=255)
                                                                                                             ttttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc                                         >  2:743726/1‑150 (MQ=255)
                                                                                                               ttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAt                                       <  2:4740832/150‑1 (MQ=255)
                                                                                                               ttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAt                                       >  2:446064/1‑150 (MQ=255)
                                                                                                                  tAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGa                                    >  2:1608755/1‑150 (MQ=255)
                                                                                                                   aaaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGATATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAg                                   >  2:2206245/1‑150 (MQ=255)
                                                                                                                    aaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgt                                  >  2:3462528/1‑150 (MQ=255)
                                                                                                                      cGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgt                                <  2:3214016/150‑1 (MQ=255)
                                                                                                                      cGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgt                                >  1:2246384/1‑150 (MQ=255)
                                                                                                                       gTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgtg                               <  2:1100095/150‑1 (MQ=255)
                                                                                                                         tgggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTa                             >  2:3314270/2‑150 (MQ=255)
                                                                                                                         aGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTa                             >  1:4682946/1‑150 (MQ=255)
                                                                                                                            gAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTaaaa                          >  2:4209951/1‑150 (MQ=255)
                                                                                                                            gAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTaaaa                          <  1:2446814/150‑1 (MQ=255)
                                                                                                                             aGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAg                         <  1:3132742/150‑1 (MQ=255)
                                                                                                                                ccgccACCCACCCGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAAGCAATGCAAACCGTTTTGCCGTGAGTGTGTAAAAGAAt                      >  2:4991998/1‑150 (MQ=255)
                                                                                                                                  gccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATcc                    <  2:4120033/150‑1 (MQ=255)
                                                                                                                                  gccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATcc                    <  1:3421615/150‑1 (MQ=255)
                                                                                                                                   ccacccacCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATccc                   <  1:3601216/150‑1 (MQ=255)
                                                                                                                                    cacccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCa                  <  2:2692706/150‑1 (MQ=255)
                                                                                                                                      cccaccAGGTTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATc                <  1:4506165/150‑1 (MQ=255)
                                                                                                                                       ccaccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCg               >  1:2752540/1‑150 (MQ=255)
                                                                                                                                       ccaccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCg               <  2:3371699/150‑1 (MQ=255)
                                                                                                                                            aGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAAt          <  1:3366486/150‑1 (MQ=255)
                                                                                                                                                  gAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATgcagca    >  2:2781727/1‑150 (MQ=255)
                                                                                                                                                  gAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATgcagca    >  2:183447/1‑150 (MQ=255)
                                                                                                                                                  gAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATgcagca    <  2:952150/150‑1 (MQ=255)
                                                                                                                                                   aaaTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCaa   >  2:1650497/1‑150 (MQ=255)
                                                                                                                                                    aaTAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCaat  <  1:455028/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCAAGATTATATGCCTCCCTGTTTTCGACGTATGCGCCAGCAAATACCGAGGCAGAGGTCAGAGAAGATAATAAAATTATTGCATTAATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAAT  >  CP009273/4057557‑4057854

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: