Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,323,686 C→A 21.6% D28Y (GAT→TAT)  basR ← response regulator in two‑component regulatory system with BasS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,323,6860CA21.6% 829.3 / 183.8 309D28Y (GAT→TAT) basRresponse regulator in two‑component regulatory system with BasS
Reads supporting (aligned to +/- strand):  ref base C (111/131);  new base A (26/41);  total (137/172)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

GATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGT  >  CP009273/4323537‑4323835
                                                                                                                                                     |                                                                                                                                                     
 aTCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCg                                                                                                                                                      <  1:1450077/150‑1 (MQ=255)
  tCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGc                                                                                                                                                     >  2:962947/1‑150 (MQ=255)
         cGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTa                                                                                                                                              >  2:2965155/1‑150 (MQ=255)
           gCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGc                                                                                                                                            >  2:3122171/1‑150 (MQ=255)
              gggTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCtt                                                                                                                                         >  2:1231089/1‑150 (MQ=255)
                gTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCg                                                                                                                                       <  1:871309/150‑1 (MQ=255)
                 tataTTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCgg                                                                                                                                      >  1:594143/1‑150 (MQ=255)
                    atTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGttt                                                                                                                                   >  1:4368146/1‑150 (MQ=255)
                     ttttttCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTg                                                                                                                                  >  2:3247823/1‑150 (MQ=255)
                      tttttCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTgc                                                                                                                                 >  1:4459966/1‑150 (MQ=255)
                      tttttCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTgc                                                                                                                                 <  1:58851/150‑1 (MQ=255)
                       ttttCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTgcg                                                                                                                                <  1:4590248/150‑1 (MQ=255)
                           cTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGcgccgc                                                                                                                            <  2:4059968/150‑1 (MQ=255)
                           cTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGcgccgc                                                                                                                            >  1:2355598/1‑150 (MQ=255)
                           cTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGcgccgc                                                                                                                            >  1:2356432/1‑150 (MQ=255)
                            tGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCgccgcc                                                                                                                           <  2:731516/150‑1 (MQ=255)
                             gCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCa                                                                                                                          <  1:448953/150‑1 (MQ=255)
                               cGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGa                                                                                                                        <  2:1108265/150‑1 (MQ=255)
                                ggATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGaa                                                                                                                       >  2:767410/1‑150 (MQ=255)
                                ggATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGaa                                                                                                                       >  2:817808/1‑150 (MQ=255)
                                 gATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAAt                                                                                                                      >  1:3304479/1‑150 (MQ=255)
                                 gATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAAt                                                                                                                      <  2:4669002/150‑1 (MQ=255)
                                     cGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGt                                                                                                                  >  2:3335786/1‑150 (MQ=255)
                                      gggCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTc                                                                                                                 >  1:1016477/1‑150 (MQ=255)
                                      gggCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTc                                                                                                                 >  2:4777350/1‑150 (MQ=255)
                                         cGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCt                                                                                                              <  1:4158370/150‑1 (MQ=255)
                                         cGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCt                                                                                                              <  1:4304300/150‑1 (MQ=255)
                                          gagaAAATGCAGTCCATCTTCGTCGGGTAACCCTACATCCAGTTCCACGAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCTTCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTg                                                                                                             <  2:2286438/150‑1 (MQ=255)
                                          gagaAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTg                                                                                                             <  1:4093994/150‑1 (MQ=255)
                                              aaaTGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAAt                                                                                                         <  1:1556430/150‑1 (MQ=255)
                                                aTGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATaa                                                                                                       >  1:2963228/1‑150 (MQ=255)
                                                aTGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATaa                                                                                                       <  2:3935507/150‑1 (MQ=255)
                                                  gCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACa                                                                                                     <  1:2082401/150‑1 (MQ=255)
                                                   cAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTATTGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAg                                                                                                    <  2:3155435/150‑1 (MQ=255)
                                                      tCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGt                                                                                                 <  1:3213216/150‑1 (MQ=255)
                                                      tCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGt                                                                                                 >  1:2053237/1‑150 (MQ=255)
                                                         aTCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATc                                                                                              <  2:3593202/150‑1 (MQ=255)
                                                          tCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCg                                                                                             <  1:983102/150‑1 (MQ=255)
                                                          tCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCg                                                                                             >  2:1551661/1‑150 (MQ=255)
                                                          tCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCg                                                                                             >  2:2406546/1‑150 (MQ=255)
                                                            ttCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTc                                                                                           <  2:1827197/150‑1 (MQ=255)
                                                            ttCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTc                                                                                           <  2:2503802/150‑1 (MQ=255)
                                                             tcgtcgGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCt                                                                                          <  2:826056/150‑1 (MQ=255)
                                                               gtcgGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTc                                                                                        >  1:3347155/1‑150 (MQ=255)
                                                                tcgGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTca                                                                                       >  2:3782626/1‑150 (MQ=255)
                                                                 cgGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTcaa                                                                                      >  1:4573255/1‑150 (MQ=255)
                                                                  gggTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTcaac                                                                                     >  1:2321577/1‑150 (MQ=255)
                                                                     tAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAAt                                                                                  <  2:4388830/150‑1 (MQ=255)
                                                                     tAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAAt                                                                                  >  1:946253/1‑150 (MQ=255)
                                                                     tAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAAt                                                                                  <  2:1607739/150‑1 (MQ=255)
                                                                      aaCCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATc                                                                                 <  2:2182650/150‑1 (MQ=255)
                                                                       aCCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCa                                                                                <  2:4481252/150‑1 (MQ=255)
                                                                        cccTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAg                                                                               >  1:2173976/1‑150 (MQ=255)
                                                                          cTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGaa                                                                             <  1:4800285/150‑1 (MQ=255)
                                                                           tAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAt                                                                            >  2:273211/1‑150 (MQ=255)
                                                                           tAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAt                                                                            <  2:640091/150‑1 (MQ=255)
                                                                            aaaTCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAtt                                                                           <  1:2639423/150‑1 (MQ=255)
                                                                             aaTCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAttt                                                                          <  1:2933661/150‑1 (MQ=255)
                                                                              aTCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAtttt                                                                         <  2:3520222/150‑1 (MQ=255)
                                                                              aTCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAtttt                                                                         >  2:510657/1‑150 (MQ=255)
                                                                               tCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTttc                                                                        <  1:3757218/150‑1 (MQ=255)
                                                                                   gTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTttcattc                                                                    >  1:4689517/1‑150 (MQ=255)
                                                                                    tACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTttcattca                                                                   <  1:273211/150‑1 (MQ=255)
                                                                                     accaccAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATtcac                                                                  <  2:3006545/150‑1 (MQ=255)
                                                                                      ccaccaGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATtcact                                                                 <  1:3576007/150‑1 (MQ=255)
                                                                                      ccaccaGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATtcact                                                                 <  1:3494357/150‑1 (MQ=255)
                                                                                         ccaGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATtcactcac                                                              >  1:813337/1‑150 (MQ=255)
                                                                                          caGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTcactcact                                                             <  1:3167915/150‑1 (MQ=255)
                                                                                            ggCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCActct                                                           >  1:66241/1‑150 (MQ=255)
                                                                                              cTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTcc                                                         <  2:1682397/150‑1 (MQ=255)
                                                                                              cTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTcc                                                         >  1:2906625/1‑150 (MQ=255)
                                                                                               tGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCt                                                        >  1:1135370/1‑150 (MQ=255)
                                                                                                gTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTg                                                       <  2:442911/150‑1 (MQ=255)
                                                                                                 tAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGc                                                      >  2:2002505/1‑150 (MQ=255)
                                                                                                 tAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGc                                                      >  2:1465788/1‑150 (MQ=255)
                                                                                                 tAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGc                                                      >  2:1464437/1‑150 (MQ=255)
                                                                                                    tGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAg                                                   >  2:1016290/1‑150 (MQ=255)
                                                                                                           ccTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTTCAg                                            >  2:2546003/1‑150 (MQ=255)
                                                                                                            cTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGa                                           <  1:3866841/150‑1 (MQ=255)
                                                                                                             tCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGaa                                          <  1:291174/150‑1 (MQ=255)
                                                                                                              cAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAat                                         >  1:1986675/1‑150 (MQ=255)
                                                                                                                      tGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCg                                 <  2:663122/150‑1 (MQ=255)
                                                                                                                      tGTTCCGCCATCCGCGCGGTTGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCg                                 <  2:751591/150‑1 (MQ=255)
                                                                                                                       gTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGc                                >  2:88708/1‑150 (MQ=255)
                                                                                                                         tCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAg                              >  1:2518619/1‑150 (MQ=255)
                                                                                                                            gCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCt                           <  1:3947970/150‑1 (MQ=255)
                                                                                                                               aTCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGt                        <  1:3251022/150‑1 (MQ=255)
                                                                                                                                tCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTa                       >  1:2251005/1‑150 (MQ=255)
                                                                                                                                tCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTa                       <  1:1762741/150‑1 (MQ=255)
                                                                                                                                tCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTa                       >  2:4600407/1‑150 (MQ=255)
                                                                                                                                tCCGCGCGGTTGTCACGCTATCGCACGCGTAGACTTCGGTTTGCGCCGCCAGAATCAGTACCTGCAATAACAGAGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTa                       <  1:1334069/150‑1 (MQ=255)
                                                                                                                                       ggTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTc                <  1:2263862/150‑1 (MQ=255)
                                                                                                                                          tGTCACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGtc             >  1:742396/1‑150 (MQ=255)
                                                                                                                                            tCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACCGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGtctc           <  2:4579037/150‑1 (MQ=255)
                                                                                                                                            tCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGtctc           <  1:960334/150‑1 (MQ=255)
                                                                                                                                            tCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGtctc           <  2:1384556/150‑1 (MQ=255)
                                                                                                                                             cACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCa          >  2:2630814/1‑150 (MQ=255)
                                                                                                                                             cACGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCa          >  2:1874834/1‑150 (MQ=255)
                                                                                                                                              aCGCTATAGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCaa         <  1:1596253/150‑1 (MQ=255)
                                                                                                                                                    tCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAg   >  2:2239459/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGT  >  CP009273/4323537‑4323835

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: