Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 549,397 G→T 85.8% intergenic (‑1/+2) lpxH ← / ← ppiB UDP‑2,3‑diacylglucosamine pyrophosphohydrolase/peptidyl‑prolyl cis‑trans isomerase B (rotamase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273549,3970GT85.8% 734.1 / 125.6 296intergenic (‑1/+2)lpxH/ppiBUDP‑2,3‑diacylglucosamine pyrophosphohydrolase/peptidyl‑prolyl cis‑trans isomerase B (rotamase B)
Reads supporting (aligned to +/- strand):  ref base G (16/26);  new base T (126/128);  total (142/154)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

Only 100 of 299 total aligned reads displayed.

TCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCT  >  CP009273/549248‑549546
                                                                                                                                                     |                                                                                                                                                     
tcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACg                                                                                                                                                       <  2:4403483/150‑1 (MQ=255)
  gtcgCCAATCCTTGCTTAAAACTGATCGCCAATATTATACCGCGGGTCGGCCTGGGGGGCTGCCCCCGCTAAAAAACGAAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTAt                                                                                                                                                     <  1:1203440/150‑1 (MQ=255)
  gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCCGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTAt                                                                                                                                                     >  2:4025864/1‑150 (MQ=255)
  gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGAt                                                                                                                                                     >  2:3958425/1‑150 (MQ=255)
    cgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTa                                                                                                                                                   >  2:874695/1‑150 (MQ=255)
    cgCAAATCAATGCTACAAACAGATCGCCAATAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGATAAAAAACGCATAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTa                                                                                                                                                   <  1:1949838/150‑1 (MQ=255)
     gCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAGCTGCAATAAAGAGAGTCGCCACGAGTAc                                                                                                                                                  >  2:227630/1‑150 (MQ=255)
     gCCAATCCATGCTTCAAACAGATCGCCAAGAATAAACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTAc                                                                                                                                                  <  1:822112/150‑1 (MQ=255)
       cAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTc                                                                                                                                                >  1:3076053/1‑150 (MQ=255)
           ccATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTa                                                                                                                                            <  2:3194250/150‑1 (MQ=255)
            cATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTaa                                                                                                                                           >  2:1836652/1‑150 (MQ=255)
            cATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTaa                                                                                                                                           >  2:2288201/1‑150 (MQ=255)
            cATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTaa                                                                                                                                           >  2:4258790/1‑150 (MQ=255)
            cATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTaa                                                                                                                                           >  1:4371374/1‑150 (MQ=255)
             aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAAc                                                                                                                                          <  2:1717554/150‑1 (MQ=255)
             aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAAc                                                                                                                                          >  2:1806779/1‑150 (MQ=255)
              tGCTTCAAAAAGATCGCCAAGAATATACAGCGCGTCGTCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACg                                                                                                                                         <  2:2752655/150‑1 (MQ=255)
                cTTCAAACAGATCGCCAATAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGt                                                                                                                                       <  1:3561404/150‑1 (MQ=255)
                 ttCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTc                                                                                                                                      <  2:1124801/150‑1 (MQ=255)
                     aaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGc                                                                                                                                  >  2:634112/1‑150 (MQ=255)
                       cAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCtt                                                                                                                                <  2:1241273/150‑1 (MQ=255)
                       cAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCtt                                                                                                                                <  2:4154200/150‑1 (MQ=255)
                         gATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTc                                                                                                                              <  1:911335/150‑1 (MQ=255)
                          aTCGCCAAGAATATACCGCGGGTCGGCCTTGCGGGCTTCCCCCGATAAAAAACGCATAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCa                                                                                                                             <  1:865606/150‑1 (MQ=255)
                            cGCCAAGAATATCCCGCGCGTCGTCGTTGCGGGCTTCCCCCCCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAAt                                                                                                                           <  1:2909593/150‑1 (MQ=255)
                               cAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGAt                                                                                                                        >  2:2339501/1‑150 (MQ=255)
                                 aGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATaa                                                                                                                      >  2:1536429/1‑150 (MQ=255)
                                  gAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAAc                                                                                                                     <  2:2172538/150‑1 (MQ=255)
                                  gAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAAc                                                                                                                     >  1:2111713/1‑150 (MQ=255)
                                     tataCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGtc                                                                                                                  >  2:2983728/1‑150 (MQ=255)
                                      ataCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGtct                                                                                                                 <  1:209073/150‑1 (MQ=255)
                                      ataCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGtct                                                                                                                 <  2:2909431/150‑1 (MQ=255)
                                      ataCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGtct                                                                                                                 <  2:699144/150‑1 (MQ=255)
                                        cccgcgcGTCGGCGTGGGGGGGTCCCCCCGATAAAAAACGAATAAACCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGtcttc                                                                                                               <  1:2663347/147‑1 (MQ=255)
                                         cAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGtcttct                                                                                                              >  1:4029742/1‑150 (MQ=255)
                                         cAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGtcttct                                                                                                              >  1:604965/1‑150 (MQ=255)
                                         cAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGtcttct                                                                                                              <  1:1800207/150‑1 (MQ=255)
                                          aGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTcttctt                                                                                                             <  1:654803/150‑1 (MQ=255)
                                              cgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTgg                                                                                                         <  2:474882/150‑1 (MQ=255)
                                              cgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTgg                                                                                                         <  2:3721934/150‑1 (MQ=255)
                                                    ccTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAACGAGTGTCGCCACTAGTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTc                                                                                                   <  2:522134/150‑1 (MQ=255)
                                                       tGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTg                                                                                                >  1:21943/1‑150 (MQ=255)
                                                          gggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTg                                                                                             >  1:3694893/1‑150 (MQ=255)
                                                           ggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGc                                                                                            >  2:198955/1‑150 (MQ=255)
                                                            gcgtCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCa                                                                                           <  2:2232064/147‑1 (MQ=255)
                                                              ttCCCCCGCTAAAAAACGCAGACAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATa                                                                                         <  2:3546263/150‑1 (MQ=255)
                                                              ttCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATa                                                                                         >  2:1382459/1‑150 (MQ=255)
                                                                   ccGCTAAAAAACGCAGAACACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCt                                                                                    <  2:4527769/150‑1 (MQ=255)
                                                                   ccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCt                                                                                    >  2:159909/1‑150 (MQ=255)
                                                                     gCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTAc                                                                                  <  2:1270908/150‑1 (MQ=255)
                                                                        aaaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAc                                                                               >  2:3830199/1‑150 (MQ=255)
                                                                         aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGACGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCTTACCGCTACGAcc                                                                              <  1:2988097/150‑1 (MQ=255)
                                                                          aaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCg                                                                             <  2:566441/150‑1 (MQ=255)
                                                                           aaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCgg                                                                            >  1:267315/1‑150 (MQ=255)
                                                                           aaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCgg                                                                            <  1:632725/150‑1 (MQ=255)
                                                                           aaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAACAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAcccg                                                                            <  2:368958/150‑3 (MQ=255)
                                                                            aaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATCCCGCTACGACCGGt                                                                           >  2:3474777/1‑150 (MQ=255)
                                                                             aCGCAGAAACCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGtt                                                                          <  1:4245188/150‑1 (MQ=255)
                                                                             aCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGtt                                                                          <  1:2740466/150‑1 (MQ=255)
                                                                              cGCAGAAACCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTg                                                                         <  1:2288201/150‑1 (MQ=255)
                                                                               gCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGc                                                                        >  1:4520788/1‑150 (MQ=255)
                                                                                cAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCa                                                                       >  2:1448483/1‑150 (MQ=255)
                                                                                     aaCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCaccacc                                                                  >  2:3534610/1‑150 (MQ=255)
                                                                                        cggcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCttt                                                               >  1:2486723/1‑150 (MQ=255)
                                                                                         ggcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTg                                                              <  1:4459610/150‑1 (MQ=255)
                                                                                         ggcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTg                                                              <  1:1527985/150‑1 (MQ=255)
                                                                                         ggcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTg                                                              >  2:1987827/1‑150 (MQ=255)
                                                                                         ggcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTg                                                              >  2:3130219/1‑150 (MQ=255)
                                                                                         ggcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTg                                                              >  1:2883879/1‑150 (MQ=255)
                                                                                          gcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGa                                                             <  1:3376911/150‑1 (MQ=255)
                                                                                             gTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAttt                                                          <  2:2734417/150‑1 (MQ=255)
                                                                                              tGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAtttt                                                         >  1:1040516/1‑150 (MQ=255)
                                                                                               gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTg                                                        >  1:1990056/1‑150 (MQ=255)
                                                                                               gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCGTTGATTTTg                                                        >  1:744/1‑150 (MQ=255)
                                                                                                  cGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCt                                                     >  1:1292548/1‑150 (MQ=255)
                                                                                                   gCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTa                                                    <  1:1409298/150‑1 (MQ=255)
                                                                                                       gTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCAc                                                >  1:151002/1‑150 (MQ=255)
                                                                                                       gTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTAGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCAc                                                <  1:956857/150‑1 (MQ=255)
                                                                                                        ttcttcCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACg                                               >  2:2546036/1‑150 (MQ=255)
                                                                                                        ttcttcCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACg                                               <  2:4098601/150‑1 (MQ=255)
                                                                                                         tcttcCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGt                                              >  2:1215907/1‑150 (MQ=255)
                                                                                                                      gaTGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAAc                                 <  1:1622001/150‑1 (MQ=255)
                                                                                                                         gAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCAc                              >  1:1188469/1‑150 (MQ=255)
                                                                                                                          aaGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACt                             >  1:2684088/1‑150 (MQ=255)
                                                                                                                          aaGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACt                             >  1:4139440/1‑150 (MQ=255)
                                                                                                                           aGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACtt                            <  1:3212139/150‑1 (MQ=255)
                                                                                                                              tCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAg                         <  2:1684012/150‑1 (MQ=255)
                                                                                                                                tGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCa                       <  2:3400098/150‑1 (MQ=255)
                                                                                                                                     tAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAacac                  <  2:3794633/150‑1 (MQ=255)
                                                                                                                                      aaaGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACg                 >  1:536890/1‑150 (MQ=255)
                                                                                                                                          attgtCGCCACTATTACTCGCTAAAGGTCACGCTTTCAAAGATAACGTCTTATTATGGCACGTCCTGGTGCATACCGCTACTAACGGTTGCAACACCTATCATTTAGTCTACAACGTCCATGCCGTCAACCACTTCAGCAAACACTCAGt             >  2:2375681/3‑150 (MQ=255)
                                                                                                                                           gtgtCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTa            <  2:2331526/150‑1 (MQ=255)
                                                                                                                                            tgtCGCCACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAg           >  1:701835/1‑150 (MQ=255)
                                                                                                                                                gCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTCCCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGcccc       >  2:2867667/1‑150 (MQ=255)
                                                                                                                                                gCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGcccc       >  2:189752/1‑150 (MQ=255)
                                                                                                                                                 ccACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCa      <  1:3888257/150‑1 (MQ=255)
                                                                                                                                                 ccACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCa      >  2:770437/1‑150 (MQ=255)
                                                                                                                                                 ccACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCa      >  2:3442741/1‑150 (MQ=255)
                                                                                                                                                  cACTATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCaa     >  1:3815047/1‑150 (MQ=255)
                                                                                                                                                     gATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCt  >  2:1194385/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCT  >  CP009273/549248‑549546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: