Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,792,473 G→A 11.8% intergenic (‑619/‑50) stpA ← / → alaE DNA binding protein, nucleoid‑associated/alanine exporter, alanine‑inducible, stress‑responsive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,792,4730GA11.8% 539.3 / 47.8 177intergenic (‑619/‑50)stpA/alaEDNA binding protein, nucleoid‑associated/alanine exporter, alanine‑inducible, stress‑responsive
Reads supporting (aligned to +/- strand):  ref base G (103/53);  new base A (13/8);  total (116/61)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

CATGTAAAAATCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGTTTTCCTC  >  CP009273/2792326‑2792621
                                                                                                                                                   |                                                                                                                                                    
cATGTAAAAATCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCt                                                                                                                                                    >  1:4503402/1‑150 (MQ=255)
 aTGTAAAAATCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCtt                                                                                                                                                   <  2:11354/150‑1 (MQ=255)
 aTGTAAAAATCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCtt                                                                                                                                                   <  2:11511/150‑1 (MQ=255)
  tGTAAAAATCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCttt                                                                                                                                                  >  2:1034885/1‑150 (MQ=255)
  tGTAAAAATCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCttt                                                                                                                                                  <  1:2992855/150‑1 (MQ=255)
  tGTAAAAATCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCttt                                                                                                                                                  >  2:2327332/1‑150 (MQ=255)
     aaaaaTCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCAt                                                                                                                                               >  1:3036735/1‑150 (MQ=255)
     aaaaaTCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCAt                                                                                                                                               >  2:1818156/1‑150 (MQ=255)
         aTCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTcc                                                                                                                                           >  2:1970132/1‑150 (MQ=255)
         aTCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTcc                                                                                                                                           >  2:1970310/1‑150 (MQ=255)
           ccACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTccct                                                                                                                                         >  2:3420426/1‑148 (MQ=255)
              caATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTaa                                                                                                                                      <  1:2675728/150‑1 (MQ=255)
               aaTGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAAc                                                                                                                                     >  1:1411953/1‑150 (MQ=255)
                  gTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATc                                                                                                                                  <  1:1413646/150‑1 (MQ=255)
                  gTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATc                                                                                                                                  >  1:327358/1‑150 (MQ=255)
                        aaaaCGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTa                                                                                                                            <  1:4861263/150‑1 (MQ=255)
                        aaaaCGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTACTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTa                                                                                                                            >  1:22128/1‑150 (MQ=255)
                             gAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCtt                                                                                                                       <  1:2663952/150‑1 (MQ=255)
                              aGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCttt                                                                                                                      <  2:1155294/150‑1 (MQ=255)
                                  gTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTttatt                                                                                                                  <  1:1503581/150‑1 (MQ=255)
                                              tttATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAg                                                                                                      >  2:4233578/1‑150 (MQ=255)
                                              tttATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAg                                                                                                      >  2:2203875/1‑150 (MQ=255)
                                               ttATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGc                                                                                                     <  1:1761307/150‑1 (MQ=255)
                                                   cTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGt                                                                                                 >  2:746701/1‑150 (MQ=255)
                                                        aaaaCTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCa                                                                                            <  2:1902106/150‑1 (MQ=255)
                                                        aaaaCTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCa                                                                                            <  2:4048043/150‑1 (MQ=255)
                                                        aaaaCTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCa                                                                                            >  2:2276283/1‑150 (MQ=255)
                                                         aaaCTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATACACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCAc                                                                                           <  1:1552433/150‑1 (MQ=255)
                                                           aCTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCg                                                                                         >  2:1476922/1‑150 (MQ=255)
                                                           aCTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCg                                                                                         >  2:1476855/1‑150 (MQ=255)
                                                           aCTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCg                                                                                         >  2:775490/1‑150 (MQ=255)
                                                            cTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGc                                                                                        >  2:233337/1‑150 (MQ=255)
                                                              gACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAg                                                                                      >  2:2441450/1‑150 (MQ=255)
                                                                cTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTc                                                                                    >  2:3766017/1‑150 (MQ=255)
                                                                 taaataaaAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCa                                                                                   >  1:1991620/1‑150 (MQ=255)
                                                                  aaataaaAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCAc                                                                                  >  1:4836246/1‑150 (MQ=255)
                                                                    ataaaAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGc                                                                                >  1:4196752/1‑150 (MQ=255)
                                                                    ataaaAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGc                                                                                >  2:2446660/1‑150 (MQ=255)
                                                                      aaaaaTTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCtt                                                                              >  1:3596010/1‑150 (MQ=255)
                                                                           tttttCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTc                                                                         >  1:1824068/1‑150 (MQ=255)
                                                                               tCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGc                                                                     <  1:4220723/150‑1 (MQ=255)
                                                                                cACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCa                                                                    >  1:796195/1‑150 (MQ=255)
                                                                                   tAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGtt                                                                 >  1:1099406/1‑150 (MQ=255)
                                                                                   tAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGtt                                                                 >  1:906418/1‑150 (MQ=255)
                                                                                    aaTTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTg                                                                >  2:4477358/1‑150 (MQ=255)
                                                                                    aaTTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTg                                                                >  2:4066666/1‑150 (MQ=255)
                                                                                    aaTTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTg                                                                >  2:84935/1‑150 (MQ=255)
                                                                                    aaTTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTg                                                                <  2:3157373/150‑1 (MQ=255)
                                                                                     aTTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGc                                                               >  2:4756039/1‑150 (MQ=255)
                                                                                       tGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAg                                                             >  2:4168345/1‑150 (MQ=255)
                                                                                       tGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAg                                                             >  2:1673714/1‑150 (MQ=255)
                                                                                       tGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAg                                                             >  2:1673719/1‑150 (MQ=255)
                                                                                        gATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGa                                                            <  2:2620721/150‑1 (MQ=255)
                                                                                        gATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGa                                                            <  2:2620619/150‑1 (MQ=255)
                                                                                         aTTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGAt                                                           >  2:2220113/1‑150 (MQ=255)
                                                                                            aGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACg                                                        >  2:3887006/1‑150 (MQ=255)
                                                                                             gTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGt                                                       >  2:70952/1‑150 (MQ=255)
                                                                                                aTAGCCAGTGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTcg                                                    >  1:4132673/1‑150 (MQ=255)
                                                                                                aTAGCCAGTGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCAATTAACATCCAATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTcg                                                    >  1:1674899/1‑150 (MQ=255)
                                                                                                 tAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTcgc                                                   >  2:2429460/1‑150 (MQ=255)
                                                                                                  aGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTcgcg                                                  >  2:394266/1‑150 (MQ=255)
                                                                                                    ccAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGAt                                                >  1:1612555/1‑150 (MQ=255)
                                                                                                     cAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCAGTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATg                                               >  2:2635966/1‑150 (MQ=255)
                                                                                                       gCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGgt                                             <  1:3582345/150‑1 (MQ=255)
                                                                                                       gCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGgt                                             <  1:2465439/150‑1 (MQ=255)
                                                                                                         gATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGgttg                                           >  2:4204759/1‑150 (MQ=255)
                                                                                                            ataCGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGgttgtta                                        >  2:317455/1‑149 (MQ=255)
                                                                                                                gCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTg                                    >  2:4693877/1‑150 (MQ=255)
                                                                                                                gCTATGCGAAAATGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTg                                    >  1:2821516/1‑150 (MQ=255)
                                                                                                                 cTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGt                                   >  1:3020319/1‑150 (MQ=255)
                                                                                                                    tGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCt                                >  1:2248419/1‑150 (MQ=255)
                                                                                                                    tGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGAGGGTTGTTTACTGTTCt                                >  1:1938164/1‑150 (MQ=255)
                                                                                                                      cGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGt                              <  1:2970056/150‑1 (MQ=255)
                                                                                                                      cGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGt                              >  2:1360773/1‑150 (MQ=255)
                                                                                                                        aaaaTGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCg                            >  1:4636751/1‑150 (MQ=255)
                                                                                                                        aaaaTGCAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCg                            >  1:4194974/1‑150 (MQ=255)
                                                                                                                          aaTGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTg                          <  2:3097563/150‑1 (MQ=255)
                                                                                                                           aTGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCGTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGa                         <  1:4805120/150‑1 (MQ=255)
                                                                                                                            tGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGaa                        <  2:3154881/150‑1 (MQ=255)
                                                                                                                              cAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACa                      <  2:849760/150‑1 (MQ=255)
                                                                                                                              cAGATGGCAATGAGATCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACa                      >  1:1143638/1‑150 (MQ=255)
                                                                                                                               aGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACAt                     <  1:2540126/150‑1 (MQ=255)
                                                                                                                               aGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACAt                     >  1:3865611/1‑150 (MQ=255)
                                                                                                                                gATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACAtg                    >  2:791926/1‑150 (MQ=255)
                                                                                                                                 aTGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACAtgt                   <  2:2020335/150‑1 (MQ=255)
                                                                                                                                  tGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACAtgtg                  <  2:1972794/150‑1 (MQ=255)
                                                                                                                                    gCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTa                <  1:412528/150‑1 (MQ=255)
                                                                                                                                     cAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTAt               >  2:4258408/1‑150 (MQ=255)
                                                                                                                                     cAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTAt               >  2:4259600/1‑150 (MQ=255)
                                                                                                                                     cAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTAt               <  2:373411/150‑1 (MQ=255)
                                                                                                                                        tGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGa            >  2:769848/1‑150 (MQ=255)
                                                                                                                                         gagaTCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGaa           <  2:3125506/150‑1 (MQ=255)
                                                                                                                                           gaTCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGt         >  1:433416/1‑150 (MQ=255)
                                                                                                                                           gaTCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGt         >  1:6465/1‑150 (MQ=255)
                                                                                                                                             tCCACTACTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGttt       >  2:57153/1‑150 (MQ=255)
                                                                                                                                              ccACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGtttt      >  1:3854778/1‑150 (MQ=255)
                                                                                                                                               cACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGTTTTc     >  1:715489/1‑150 (MQ=255)
                                                                                                                                                 cTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGTTTTCct   <  2:3036735/150‑1 (MQ=255)
                                                                                                                                                  tGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGTTTTCctc  <  1:3420426/150‑1 (MQ=255)
                                                                                                                                                  tGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGTTTTCctc  <  1:2203875/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATGTAAAAATCCACAATGTACAAAAAACGAGCCGTTACGGAATATTTTATCTACAAAAACTGACTAAATAAAAATTTTTCACTAATTGATTAGTCATAGCCAGCGATATACGCTATGCGAAAATGCAGATGGCAATGAGATCCACTGCTTTCATCTCCATTAACATCCCATTACGCTTTTATTAAGGAGCATTAGCATGTTCTCACCGCAGTCACGCTTGCGTCATGCAGTTGCAGATACGTTCGCGATGGTTGTTTACTGTTCTGTCGTGAACATGTGTATTGAAGTTTTCCTC  >  CP009273/2792326‑2792621

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: