Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,422 G→T 76.5% 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,4220GT76.5% 320.9 / 143.7 183intergenic (‑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 (23/20);  new base T (71/69);  total (94/89)
Fisher's exact test for biased strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

Only 100 of 184 total aligned reads displayed.

CTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTT  >  CP009273/647273‑647568
                                                                                                                                                     |                                                                                                                                                  
cTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAg                                                                                                                                                    >  2:421845/1‑150 (MQ=255)
   aaCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGta                                                                                                                                                 <  1:230868/150‑1 (MQ=255)
   aaCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGta                                                                                                                                                 >  1:2890080/1‑150 (MQ=255)
          tACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCttt                                                                                                                                          <  1:2496281/150‑1 (MQ=255)
          tACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCttt                                                                                                                                          <  2:2965505/150‑1 (MQ=255)
             gcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtc                                                                                                                                       <  1:138168/150‑1 (MQ=255)
              cgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtct                                                                                                                                      <  2:894780/150‑1 (MQ=255)
               gggTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTtctc                                                                                                                                     >  2:1402238/1‑150 (MQ=255)
               gggTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtctc                                                                                                                                     <  1:4374916/150‑1 (MQ=255)
               gggTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTtctc                                                                                                                                     >  2:1535525/1‑150 (MQ=255)
                 gTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCt                                                                                                                                   <  1:696746/150‑1 (MQ=255)
                 gTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCt                                                                                                                                   <  2:38154/150‑1 (MQ=255)
                 gTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCt                                                                                                                                   <  1:528591/150‑1 (MQ=255)
                   gAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCttt                                                                                                                                 <  1:383460/150‑1 (MQ=255)
                     aaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCg                                                                                                                               <  1:2215846/150‑1 (MQ=255)
                     aaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCg                                                                                                                               <  1:582789/150‑1 (MQ=255)
                        tatCATTGCCGTACATAATAAATAGTATCCTGAAGGTGCATGTTGTTTTCGATTTGCAAATAATGTTGTTCAATGTTGCAAAATTATAAACTTTTATTTTCACTTGGTAGTATGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCta                                                                                                                            <  1:366355/150‑1 (MQ=255)
                         atCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCtat                                                                                                                           >  1:2174995/1‑150 (MQ=255)
                         atCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCtat                                                                                                                           <  2:1682726/150‑1 (MQ=255)
                           cATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATAc                                                                                                                         >  2:2426285/1‑150 (MQ=255)
                            aTTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACa                                                                                                                        >  2:607819/1‑150 (MQ=255)
                             ttGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACAt                                                                                                                       <  1:2358071/150‑1 (MQ=255)
                                ccGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCt                                                                                                                    >  1:721958/1‑150 (MQ=255)
                                 cGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTa                                                                                                                   <  1:3458270/150‑1 (MQ=255)
                                  gAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTaa                                                                                                                  >  2:2513693/1‑150 (MQ=255)
                                      ataataAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAg                                                                                                              <  2:2605853/150‑1 (MQ=255)
                                        aataaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTa                                                                                                            >  2:1690193/1‑150 (MQ=255)
                                         ataaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAt                                                                                                           <  2:682253/150‑1 (MQ=255)
                                         ataaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAt                                                                                                           <  2:708065/150‑1 (MQ=255)
                                          taaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAtt                                                                                                          >  2:1687821/1‑150 (MQ=255)
                                           aaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAttt                                                                                                         <  2:3080255/150‑1 (MQ=255)
                                           aaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAttt                                                                                                         >  1:2948654/1‑150 (MQ=255)
                                            aataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTc                                                                                                        >  1:252828/1‑150 (MQ=255)
                                            aataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTc                                                                                                        <  2:94826/150‑1 (MQ=255)
                                             ataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCg                                                                                                       >  2:276453/1‑150 (MQ=255)
                                              taGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCgg                                                                                                      <  1:1131334/150‑1 (MQ=255)
                                                  aTCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTg                                                                                                  >  2:3272476/1‑150 (MQ=255)
                                                       gAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtt                                                                                             >  1:800927/1‑150 (MQ=255)
                                                        aaGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAttt                                                                                            >  1:2292799/1‑150 (MQ=255)
                                                        aaGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAttt                                                                                            >  2:626131/1‑150 (MQ=255)
                                                        aaGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAttt                                                                                            >  1:576725/1‑150 (MQ=255)
                                                        aaGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAttt                                                                                            >  1:2836145/1‑150 (MQ=255)
                                                         aGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtttt                                                                                           <  1:3691345/150‑1 (MQ=255)
                                                          ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTg                                                                                          <  2:2935626/150‑1 (MQ=255)
                                                          ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTg                                                                                          <  1:3490457/150‑1 (MQ=255)
                                                          ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTg                                                                                          <  2:2783826/150‑1 (MQ=255)
                                                              cATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATAAATCCTAAGGAGTATTTCGGCGTGAAATTTTGattt                                                                                      >  1:1184741/1‑150 (MQ=255)
                                                               aTGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGattta                                                                                     >  1:1702818/1‑150 (MQ=255)
                                                               aTGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGattta                                                                                     <  2:3364389/150‑1 (MQ=255)
                                                                tgttgtTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGatttat                                                                                    <  1:3823636/150‑1 (MQ=255)
                                                                 gttgttATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGatttatt                                                                                   >  1:2036303/1‑150 (MQ=255)
                                                                   tgttATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAGTTTTGATTTATTTc                                                                                 >  1:3758488/1‑150 (MQ=255)
                                                                      tATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTcaca                                                                              >  2:2783200/1‑150 (MQ=255)
                                                                         cGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATag                                                                           >  1:40336/1‑150 (MQ=255)
                                                                           aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATagag                                                                         <  2:3763261/150‑1 (MQ=255)
                                                                           aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATagag                                                                         >  1:4114675/1‑150 (MQ=255)
                                                                            tttGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGt                                                                        >  2:3750066/1‑150 (MQ=255)
                                                                               gCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAg                                                                     >  2:3358658/1‑150 (MQ=255)
                                                                                  aCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTgg                                                                  >  2:2337758/1‑150 (MQ=255)
                                                                                  aCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTgg                                                                  >  2:501921/1‑150 (MQ=255)
                                                                                  aCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTgg                                                                  >  2:501716/1‑150 (MQ=255)
                                                                                   cGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGt                                                                 <  2:1637339/150‑1 (MQ=255)
                                                                                   cGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGt                                                                 <  2:3916630/150‑1 (MQ=255)
                                                                                   cGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGt                                                                 <  1:421845/150‑1 (MQ=255)
                                                                                     aaTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATGGAGTTAGTGGttt                                                               >  2:1177241/1‑150 (MQ=255)
                                                                                      aTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGtttt                                                              >  1:4495397/1‑150 (MQ=255)
                                                                                      aTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGtttt                                                              >  2:3822473/1‑150 (MQ=255)
                                                                                          tgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctat                                                          >  2:2822/1‑150 (MQ=255)
                                                                                          tgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctat                                                          >  2:1215122/1‑150 (MQ=255)
                                                                                            ttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctattt                                                        >  1:3009759/1‑150 (MQ=255)
                                                                                              cAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtttatttat                                                      <  2:864504/150‑1 (MQ=255)
                                                                                                aTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattt                                                    >  1:4159226/1‑150 (MQ=255)
                                                                                                 tGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtttatttattta                                                   <  1:4053052/150‑1 (MQ=255)
                                                                                                  gTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTaa                                                  <  2:2174995/150‑1 (MQ=255)
                                                                                                    tGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTCCTATTTATTTAATg                                                <  2:721958/150‑1 (MQ=255)
                                                                                                      cAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGAt                                              >  1:115474/1‑150 (MQ=255)
                                                                                                       aaaCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAtt                                             <  1:4043428/150‑1 (MQ=255)
                                                                                                          cTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTa                                          >  2:3848937/1‑150 (MQ=255)
                                                                                                            gATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAg                                        >  1:4207184/1‑150 (MQ=255)
                                                                                                             aTAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGt                                       <  2:1868119/150‑1 (MQ=255)
                                                                                                             aTAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGt                                       <  1:4010150/150‑1 (MQ=255)
                                                                                                             aTAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGt                                       <  2:211546/150‑1 (MQ=255)
                                                                                                              tAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtt                                      >  2:1118962/1‑150 (MQ=255)
                                                                                                                aCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttt                                    <  1:1294316/150‑1 (MQ=255)
                                                                                                                aCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGtttt                                    >  2:3380093/1‑150 (MQ=255)
                                                                                                                aCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGtttt                                    >  2:4471704/1‑150 (MQ=255)
                                                                                                                  cTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttttt                                  >  2:1261569/1‑150 (MQ=255)
                                                                                                                   ttttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTtta                                 <  2:3997508/150‑1 (MQ=255)
                                                                                                                        ttttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaattaa                            >  2:3560425/1‑150 (MQ=255)
                                                                                                                          ttCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATg                          <  1:3947034/150‑1 (MQ=255)
                                                                                                                           tCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGt                         >  1:426495/1‑150 (MQ=255)
                                                                                                                            cACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATATCGGCGTGAAAATTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAATTTAAGTTTTTTAATTAATGTa                        >  1:1709238/1‑150 (MQ=255)
                                                                                                                            cACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTa                        >  1:2070915/1‑150 (MQ=255)
                                                                                                                                     gaAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGaaa               >  1:3707913/1‑150 (MQ=255)
                                                                                                                                         gggggTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGAc           <  2:4103585/150‑1 (MQ=255)
                                                                                                                                         gggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGAc           <  1:1594345/150‑1 (MQ=255)
                                                                                                                                          ggggTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACt          <  2:2036303/150‑1 (MQ=255)
                                                                                                                                             gTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATTACTCGc       >  1:3028206/1‑150 (MQ=255)
                                                                                                                                              tGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTAACATAGAGTTAGTGGTTTTCTATTTATTTAATTATTTTAAGTTTTTTAATTAATGTAATTAATAAATGACTCgaa      >  2:2736520/1‑148 (MQ=255)
                                                                                                                                                  cGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATTTAATTACGAAATTACTCGCAGGtt  >  2:4224441/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: