Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,422 G→T 91.7% intergenic (‑110/‑269) citC ← / → citA [citrate [pro‑3S]‑lyase] ligase/sensory histidine kinase in two‑component regulatory system with CitB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273647,4220GT91.7% 841.3 / 64.5 287intergenic (‑110/‑269)citC/citA[citrate [pro‑3S]‑lyase] ligase/sensory histidine kinase in two‑component regulatory system with CitB
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base T (121/142);  total (133/154)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 288 total aligned reads displayed.

CTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTT  >  CP009273/647273‑647568
                                                                                                                                                     |                                                                                                                                                  
   aaCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGta                                                                                                                                                 <  2:3513389/150‑1 (MQ=255)
    aCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGtat                                                                                                                                                <  2:348321/150‑1 (MQ=255)
     cGTTTTACGCGTGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCTATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGtata                                                                                                                                               <  2:2031120/150‑1 (MQ=255)
        tttACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCt                                                                                                                                            <  1:3512198/150‑1 (MQ=255)
            cgcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAGGGGGGTGATCGATGTATATCttttt                                                                                                                                        >  1:532333/1‑150 (MQ=255)
            cgcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATtttttt                                                                                                                                        >  2:3162847/1‑150 (MQ=255)
            cgcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCttttt                                                                                                                                        <  2:4184534/150‑1 (MQ=255)
            cgcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCttttt                                                                                                                                        <  2:4502229/150‑1 (MQ=255)
             gcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTtc                                                                                                                                       >  1:1236088/1‑150 (MQ=255)
             gcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtc                                                                                                                                       <  1:1116208/150‑1 (MQ=255)
              cgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtct                                                                                                                                      <  1:3559776/150‑1 (MQ=255)
                ggTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTcc                                                                                                                                    >  1:2642539/1‑150 (MQ=255)
                  tGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCtt                                                                                                                                  <  1:2035990/150‑1 (MQ=255)
                   gAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCttt                                                                                                                                 >  1:1866992/1‑150 (MQ=255)
                     aaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCg                                                                                                                               <  2:2538806/150‑1 (MQ=255)
                     aaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCg                                                                                                                               >  1:2618245/1‑150 (MQ=255)
                       atatCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCt                                                                                                                             <  1:1924736/150‑1 (MQ=255)
                        tatCATTGCCGAACATAATAAATAGTATCCTGAAAGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCta                                                                                                                            >  2:1755338/1‑150 (MQ=255)
                         atCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCtat                                                                                                                           <  1:2266399/150‑1 (MQ=255)
                                 cGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTa                                                                                                                   <  2:3630114/150‑1 (MQ=255)
                                  gAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTaa                                                                                                                  <  2:3307778/150‑1 (MQ=255)
                                  gAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTaa                                                                                                                  <  1:3545004/150‑1 (MQ=255)
                                   aaCATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAg                                                                                                                 >  1:2736115/1‑150 (MQ=255)
                                   aaCATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAg                                                                                                                 >  2:2135630/1‑150 (MQ=255)
                                      ataataAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAg                                                                                                              <  2:591553/150‑1 (MQ=255)
                                            aataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTc                                                                                                        >  2:1147997/1‑150 (MQ=255)
                                             ataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCg                                                                                                       <  1:3654146/150‑1 (MQ=255)
                                             ataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCg                                                                                                       >  1:4003706/1‑150 (MQ=255)
                                               aGTATCCTGCAGGTGCATGTTGTTATCGATTTGCCACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGc                                                                                                     <  1:4268340/150‑1 (MQ=255)
                                                  aTCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTg                                                                                                  <  2:2868334/150‑1 (MQ=255)
                                                    ccTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGaa                                                                                                >  1:3250337/1‑150 (MQ=255)
                                                    ccTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGaa                                                                                                <  1:1621891/150‑1 (MQ=255)
                                                      tGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAt                                                                                              >  1:2961165/1‑150 (MQ=255)
                                                      tGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAt                                                                                              <  1:900416/150‑1 (MQ=255)
                                                       gAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCTCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtt                                                                                             <  2:1262227/150‑1 (MQ=255)
                                                        aaGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAttt                                                                                            >  2:4902742/1‑150 (MQ=255)
                                                        aaGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAttt                                                                                            <  2:940312/150‑1 (MQ=255)
                                                          ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTg                                                                                          <  2:3948913/150‑1 (MQ=255)
                                                          ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTg                                                                                          <  2:771893/150‑1 (MQ=255)
                                                            tGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGat                                                                                        >  1:2296782/1‑150 (MQ=255)
                                                            tGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGat                                                                                        <  2:2673913/150‑1 (MQ=255)
                                                                   tgttATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTc                                                                                 <  1:625148/150‑1 (MQ=255)
                                                                    gttATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTca                                                                                >  1:2076069/1‑150 (MQ=255)
                                                                    gttATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTca                                                                                >  2:4461448/1‑150 (MQ=255)
                                                                       aTCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACAt                                                                             >  2:2224629/1‑150 (MQ=255)
                                                                              tGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTa                                                                      >  1:1578715/1‑150 (MQ=255)
                                                                                  aCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTgg                                                                  >  2:2039044/1‑150 (MQ=255)
                                                                                     aaTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttt                                                               >  2:2778215/1‑150 (MQ=255)
                                                                                      aTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGtttt                                                              <  1:3792443/150‑1 (MQ=255)
                                                                                       tgttgtTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttttc                                                             <  2:226663/150‑2 (MQ=255)
                                                                                       tgttgtTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttttt                                                             <  2:14562/150‑1 (MQ=255)
                                                                                        gttgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttttct                                                            <  1:487201/150‑3 (MQ=255)
                                                                                         ttgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTACGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtcta                                                           <  2:475564/150‑1 (MQ=255)
                                                                                           gttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatt                                                         <  2:1163865/150‑1 (MQ=255)
                                                                                               taTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGATAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttatt                                                     <  2:4978010/149‑1 (MQ=255)
                                                                                               aaTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttatt                                                     >  2:667549/1‑150 (MQ=255)
                                                                                                aTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattt                                                    <  1:2343482/150‑1 (MQ=255)
                                                                                                 tGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattta                                                   <  2:4816918/150‑1 (MQ=255)
                                                                                                 tGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattta                                                   <  2:1610179/150‑1 (MQ=255)
                                                                                                 tGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattta                                                   >  2:1693514/1‑150 (MQ=255)
                                                                                                  gTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTaa                                                  <  1:3272326/150‑1 (MQ=255)
                                                                                                   ttGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAAt                                                 <  2:2322152/150‑1 (MQ=255)
                                                                                                     gCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGa                                               <  2:5084828/150‑1 (MQ=255)
                                                                                                     gCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGa                                               <  1:3845472/150‑1 (MQ=255)
                                                                                                     gCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGa                                               <  2:1084917/150‑1 (MQ=255)
                                                                                                      cAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAt                                              >  2:4663987/1‑150 (MQ=255)
                                                                                                         aCTGATAACCTTTTATTTTCACTTGGGAGAATGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGATTTAGTGGTTTTCTATTTATTTAATGAtttt                                           <  2:532333/150‑1 (MQ=255)
                                                                                                           tGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTaa                                         <  2:3601185/150‑1 (MQ=255)
                                                                                                            gATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAg                                        <  1:2806277/150‑1 (MQ=255)
                                                                                                              tAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtt                                      <  1:3737857/150‑1 (MQ=255)
                                                                                                              tAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGtt                                      >  1:4133826/1‑150 (MQ=255)
                                                                                                                aCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttt                                    <  2:2961165/150‑1 (MQ=255)
                                                                                                                 ccTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGttttt                                   >  2:2103755/1‑150 (MQ=255)
                                                                                                                   ttttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTtta                                 >  2:2568971/1‑150 (MQ=255)
                                                                                                                     ttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaat                               <  1:1690501/150‑1 (MQ=255)
                                                                                                                         tttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTtaattaat                           >  2:1402406/1‑150 (MQ=255)
                                                                                                                             aCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTaa                       >  1:3465486/1‑150 (MQ=255)
                                                                                                                             aCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTaa                       <  2:2988572/150‑1 (MQ=255)
                                                                                                                             aCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTaa                       >  1:4375924/1‑150 (MQ=255)
                                                                                                                               ttGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAAtt                     <  2:4141211/150‑1 (MQ=255)
                                                                                                                                tGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTa                    >  1:1050224/1‑150 (MQ=255)
                                                                                                                                tGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTa                    <  2:386314/150‑1 (MQ=255)
                                                                                                                                 gggAGAAAGGGTGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTAc                   <  1:4620247/150‑1 (MQ=255)
                                                                                                                                   gagaAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGa                 >  1:2439143/1‑150 (MQ=255)
                                                                                                                                       aaGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg             >  1:1702355/1‑150 (MQ=255)
                                                                                                                                       aaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg             <  1:1333977/150‑1 (MQ=255)
                                                                                                                                       aaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg             >  2:2474740/1‑150 (MQ=255)
                                                                                                                                        aGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGa            <  2:3080069/150‑1 (MQ=255)
                                                                                                                                         gggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGAc           <  2:4053887/150‑1 (MQ=255)
                                                                                                                                         gggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGAc           <  1:1749591/150‑1 (MQ=255)
                                                                                                                                          ggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACt          >  1:4151750/1‑150 (MQ=255)
                                                                                                                                          ggggTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACt          >  1:1205412/1‑150 (MQ=255)
                                                                                                                                           gggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTc         >  1:663932/1‑150 (MQ=255)
                                                                                                                                            ggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCg        >  1:1016284/1‑150 (MQ=255)
                                                                                                                                             gTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGc       >  1:4907265/1‑150 (MQ=255)
                                                                                                                                              tGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCa      <  1:3351180/150‑1 (MQ=255)
                                                                                                                                              tGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCa      <  2:1891346/150‑1 (MQ=255)
                                                                                                                                              tGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCa      <  2:2444602/150‑1 (MQ=255)
                                                                                                                                                aTCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAgg    >  1:4524871/1‑150 (MQ=255)
                                                                                                                                                 tCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGt   >  1:5040019/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTT  >  CP009273/647273‑647568

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: