Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,509 T→C 94.2% intergenic (‑197/‑182) 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,5090TC94.2% 880.0 / 30.6 242intergenic (‑197/‑182)citC/citA[citrate [pro‑3S]‑lyase] ligase/sensory histidine kinase in two‑component regulatory system with CitB
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base C (109/119);  total (117/125)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 267 total aligned reads displayed.

TGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTAATTACATCCG  >  CP009273/647360‑647656
                                                                                                                                                     |                                                                                                                                                   
tgttgtTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttttt                                                                                                                                                     <  2:14562/150‑1 (MQ=255)
  ttgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTACGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtcta                                                                                                                                                   <  2:475564/150‑1 (MQ=255)
   tgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctat                                                                                                                                                  >  2:2546056/1‑150 (MQ=255)
    gttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatt                                                                                                                                                 <  1:2263023/150‑1 (MQ=255)
    gttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatt                                                                                                                                                 <  2:1163865/150‑1 (MQ=255)
        aaTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttatt                                                                                                                                             >  2:667549/1‑150 (MQ=255)
         aTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattt                                                                                                                                            >  2:677134/1‑150 (MQ=255)
         aTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattt                                                                                                                                            <  1:2343482/150‑1 (MQ=255)
          tGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattta                                                                                                                                           <  2:4816624/150‑1 (MQ=255)
          tGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattta                                                                                                                                           <  2:1610179/150‑1 (MQ=255)
           gTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTaa                                                                                                                                          <  2:3621943/150‑1 (MQ=255)
           gTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTaa                                                                                                                                          <  1:3272326/150‑1 (MQ=255)
             tGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATg                                                                                                                                        >  1:1911147/1‑150 (MQ=255)
              gCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGa                                                                                                                                       <  1:3845472/150‑1 (MQ=255)
                aaaCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGAtt                                                                                                                                     >  2:2913366/1‑150 (MQ=255)
                  aCTGATAACCTTTTATTTTCACTTGGGAGAATGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGATTTAGTGGTTTTCTATTTATTTAATGAtttt                                                                                                                                   <  2:532333/150‑1 (MQ=255)
                       tAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtt                                                                                                                              <  1:3737864/150‑1 (MQ=255)
                          ccTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGttttt                                                                                                                           >  1:4974728/1‑150 (MQ=255)
                          ccTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGttttt                                                                                                                           >  2:2103755/1‑150 (MQ=255)
                           cTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttttt                                                                                                                          <  1:3545206/150‑1 (MQ=255)
                            ttttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTtta                                                                                                                         >  2:2568971/1‑150 (MQ=255)
                              ttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaat                                                                                                                       <  1:1690501/150‑1 (MQ=255)
                               tATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTGTTCTATTTATTTAATGATTTTAAGTTTTttaatt                                                                                                                      >  2:4085858/1‑150 (MQ=255)
                                aTTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaatta                                                                                                                     >  1:663285/1‑150 (MQ=255)
                                  tttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTtaattaat                                                                                                                   >  2:1402406/1‑150 (MQ=255)
                                    tCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGt                                                                                                                 >  1:1567791/1‑150 (MQ=255)
                                    tCACTTGGGAGAAAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGt                                                                                                                 >  1:4129092/1‑150 (MQ=255)
                                      aCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTaa                                                                                                               <  2:2988526/150‑1 (MQ=255)
                                      aCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTaa                                                                                                               >  1:4375924/1‑150 (MQ=255)
                                      aCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTaa                                                                                                               >  1:3465486/1‑150 (MQ=255)
                                        ttGGGAGATAGGGGGTGATCGATGTATATTTTTTTCTCCTTTCGATATACATCATAAGGAGTATTTCTGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGATTTTTAATTAATGTAttt                                                                                                             <  2:4141257/150‑1 (MQ=255)
                                         tGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTa                                                                                                            <  1:4902742/150‑1 (MQ=255)
                                         tGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTa                                                                                                            >  1:1050224/1‑150 (MQ=255)
                                               aaaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAAt                                                                                                      <  2:2759463/150‑1 (MQ=255)
                                                aaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg                                                                                                     >  2:2474740/1‑150 (MQ=255)
                                                aaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg                                                                                                     <  1:1333977/150‑1 (MQ=255)
                                                 aGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGa                                                                                                    <  2:3080069/150‑1 (MQ=255)
                                                  gggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGAc                                                                                                   <  1:1749591/150‑1 (MQ=255)
                                                  gggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGAc                                                                                                   <  2:4591614/150‑1 (MQ=255)
                                                     ggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCg                                                                                                >  1:1016284/1‑150 (MQ=255)
                                                       tGATCGATGTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCa                                                                                              <  1:3351180/150‑1 (MQ=255)
                                                       tGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCa                                                                                              <  2:1891346/150‑1 (MQ=255)
                                                       tGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCa                                                                                              >  1:3142639/1‑150 (MQ=255)
                                                       tGATCGATGTATATCTTTTTCTCCTTTCGCTATACACCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCa                                                                                              <  1:2879713/150‑1 (MQ=255)
                                                           cGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGtt                                                                                          <  1:3622709/150‑1 (MQ=255)
                                                           cGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGtt                                                                                          >  2:3292883/1‑150 (MQ=255)
                                                            gatgTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGttt                                                                                         <  1:892659/147‑1 (MQ=255)
                                                            gatgTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGttt                                                                                         >  2:4970638/4‑150 (MQ=255)
                                                             atgTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTa                                                                                        >  1:2748537/3‑150 (MQ=255)
                                                             atgTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTa                                                                                        >  2:1051059/3‑150 (MQ=255)
                                                              tgTATATTTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTaa                                                                                       >  2:4344764/2‑150 (MQ=255)
                                                               ctataTCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAg                                                                                      <  2:2346246/149‑1 (MQ=255)
                                                                 atatCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGtg                                                                                    >  2:2451348/1‑150 (MQ=255)
                                                                  tatCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGtga                                                                                   >  2:426744/1‑150 (MQ=255)
                                                                   atCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGtgat                                                                                  >  2:2165664/1‑150 (MQ=255)
                                                                        tttCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGtgatttaat                                                                             >  2:1585586/1‑150 (MQ=255)
                                                                        tttCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGtgatttaat                                                                             <  2:812310/150‑1 (MQ=255)
                                                                         ttCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGtgatttaatt                                                                            <  1:356064/150‑1 (MQ=255)
                                                                               tttCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGtgatttaattgattta                                                                      >  2:3272834/1‑150 (MQ=255)
                                                                                   gCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGa                                                                  >  1:1764014/1‑150 (MQ=255)
                                                                                    cTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGaa                                                                 >  2:3498965/1‑150 (MQ=255)
                                                                                     tataCATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAAt                                                                <  2:595104/150‑1 (MQ=255)
                                                                                     tataCATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAAt                                                                <  1:5080120/150‑1 (MQ=255)
                                                                                     tataCATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAAt                                                                >  1:3177562/1‑150 (MQ=255)
                                                                                      ataCATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATa                                                               <  2:1236033/150‑1 (MQ=255)
                                                                                      ataCATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATa                                                               <  2:1769601/150‑1 (MQ=255)
                                                                                          aTCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAAt                                                           <  1:1264362/150‑1 (MQ=255)
                                                                                              tAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGc                                                       >  2:2075048/1‑150 (MQ=255)
                                                                                                  gAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACg                                                   <  1:3687957/150‑1 (MQ=255)
                                                                                                  gAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACg                                                   <  1:2124742/150‑1 (MQ=255)
                                                                                                           ggCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCAtaatt                                          <  1:2955524/150‑1 (MQ=255)
                                                                                                             cGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCAtaattaa                                        <  2:4464184/150‑1 (MQ=255)
                                                                                                                gAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAatat                                     <  2:2642539/150‑1 (MQ=255)
                                                                                                                  aaTTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCt                                   <  2:4436127/150‑1 (MQ=255)
                                                                                                                  aaTTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTAAGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCt                                   >  2:2717349/1‑150 (MQ=255)
                                                                                                                   aTTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTa                                  >  1:1604279/1‑150 (MQ=255)
                                                                                                                       tGATTTATTTCACATAGAGTTAGTTGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTa                              <  1:2201830/150‑1 (MQ=255)
                                                                                                                        gATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTAt                             <  1:4538009/150‑1 (MQ=255)
                                                                                                                            tatttCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTg                         >  2:5098861/1‑150 (MQ=255)
                                                                                                                             atttCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGa                        <  2:36503/150‑1 (MQ=255)
                                                                                                                              tttCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGAt                       <  1:4493627/150‑1 (MQ=255)
                                                                                                                                tCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTc                     >  2:4622305/1‑150 (MQ=255)
                                                                                                                                tCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTc                     <  1:2568971/150‑1 (MQ=255)
                                                                                                                                 cacaTAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCa                    >  2:2985221/1‑150 (MQ=255)
                                                                                                                                  acaTAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCaa                   <  1:1174807/150‑1 (MQ=255)
                                                                                                                                    aTAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACa                 >  2:3984382/1‑150 (MQ=255)
                                                                                                                                     tAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACAt                <  2:2439143/150‑1 (MQ=255)
                                                                                                                                     tAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACAt                >  2:165394/1‑150 (MQ=255)
                                                                                                                                      atagTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCTCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCAACTATAATAATTAATATATATGTATTTTTATTCAAAAtt               >  2:895201/3‑150 (MQ=255)
                                                                                                                                        agTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACAtttt             <  1:251571/150‑1 (MQ=255)
                                                                                                                                        agTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACAtttt             >  1:4621243/1‑150 (MQ=255)
                                                                                                                                         gTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTa            >  1:3751764/1‑150 (MQ=255)
                                                                                                                                          ttAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCATGTTTAAGTGATTTAATTTATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTaa           <  1:557659/150‑1 (MQ=255)
                                                                                                                                          ttAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTaa           <  1:2397576/150‑1 (MQ=255)
                                                                                                                                           tAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTAAt          <  1:2073561/150‑1 (MQ=255)
                                                                                                                                            aGTGGTTTTCTATTTATTTAATAATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTAAtt         >  2:4715763/1‑150 (MQ=255)
                                                                                                                                              tGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTAATTAc       >  2:1186655/1‑150 (MQ=255)
                                                                                                                                               ggTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTAATAACa      >  2:1531149/1‑150 (MQ=255)
                                                                                                                                                gTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTAATTACAt     <  1:1613551/150‑1 (MQ=255)
                                                                                                                                                  tttctATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTAATTACATcc   >  1:4293625/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTAATTACATCCG  >  CP009273/647360‑647656

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: