Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 287,084 C→A 6.7% intergenic (‑27/‑31) insA ← / → yagJ IS1 repressor TnpA/CP4‑6 prophage; uncharacterized protein;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273287,0840CA 6.7% 779.7 / 23.5 206intergenic (‑27/‑31)insA/yagJIS1 repressor TnpA/CP4‑6 prophage; uncharacterized protein;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (99/93);  new base A (10/4);  total (109/97)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

GATAGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTACCGTTGC  >  CP009273/286967‑287233
                                                                                                                     |                                                                                                                                                     
gaTAGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAt                                                                                                                       >  1:597888/1‑150 (MQ=31)
gaTAGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATAAACTAAATCAGTAAGTTGGCAGCATCACCGAAt                                                                                                                       <  2:3054081/150‑1 (MQ=14)
 aTAGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAtt                                                                                                                      >  2:3830296/1‑150 (MQ=32)
   aGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAtttt                                                                                                                    <  1:4212336/150‑1 (MQ=33)
     cgcTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCa                                                                                                                  >  2:1926675/1‑150 (MQ=35)
        tGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGcc                                                                                                               >  2:441796/1‑150 (MQ=38)
        tGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGcc                                                                                                               <  2:1548651/150‑1 (MQ=38)
             tCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAAAAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGAt                                                                                                          <  2:1491872/150‑1 (MQ=255)
               cGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGt                                                                                                        >  1:1513012/1‑150 (MQ=255)
                 gCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGttt                                                                                                      >  2:1661040/1‑150 (MQ=255)
                      gCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGt                                                                                                 >  1:433910/1‑150 (MQ=255)
                      gCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGt                                                                                                 >  2:1068326/1‑150 (MQ=255)
                        ttttGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACCCTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTAc                                                                                               <  2:229283/150‑1 (MQ=255)
                           tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                                                                                            <  2:357365/150‑1 (MQ=255)
                           tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                                                                                            >  2:3385587/1‑150 (MQ=255)
                           tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATAAACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                                                                                            >  2:881738/1‑150 (MQ=255)
                           tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATAAACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                                                                                            >  2:1797224/1‑150 (MQ=255)
                           tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATAAACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                                                                                            >  1:1812072/1‑150 (MQ=255)
                             ccGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATg                                                                                          <  2:2034203/150‑1 (MQ=255)
                             ccGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATg                                                                                          <  1:2678468/150‑1 (MQ=255)
                             ccGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATg                                                                                          <  1:504141/150‑1 (MQ=255)
                               gTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTc                                                                                        >  1:287904/1‑150 (MQ=255)
                                 tACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                                                                                      <  2:1001801/150‑1 (MQ=255)
                                   cGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAat                                                                                    >  2:2527569/1‑150 (MQ=255)
                                     caccacCCCGTCAGTAGCTGAACTGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAatat                                                                                  <  2:1483871/150‑1 (MQ=255)
                                     caccacCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAatat                                                                                  >  1:2547044/1‑150 (MQ=255)
                                        caccCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTg                                                                               >  2:3526732/1‑150 (MQ=255)
                                         accCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTgg                                                                              >  1:2146182/1‑150 (MQ=255)
                                          ccccGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGc                                                                             <  2:1988901/150‑1 (MQ=255)
                                           cccGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTCAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGcc                                                                            >  2:4258762/1‑150 (MQ=255)
                                           cccGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGcc                                                                            >  2:3112994/1‑150 (MQ=255)
                                           cccGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGcc                                                                            >  2:3450178/1‑150 (MQ=255)
                                              gTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGAt                                                                         <  2:1217372/150‑1 (MQ=255)
                                               tCAGTAGCTGAACTGGAGGGACAGCTGATAGAAACAGTAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa                                                                        <  2:3982386/150‑1 (MQ=255)
                                                cAGTAGCTTAACTGGTGGGACAGCTGATAGAAACAGACGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATaa                                                                       <  1:1462732/150‑1 (MQ=255)
                                                 aGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAg                                                                      <  2:2813964/150‑1 (MQ=255)
                                                 aGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAg                                                                      >  1:3663202/1‑150 (MQ=255)
                                                    aGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGtt                                                                   <  1:2408539/150‑1 (MQ=255)
                                                     gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTg                                                                  <  2:1561011/150‑1 (MQ=255)
                                                      cTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGc                                                                 >  2:1150409/1‑150 (MQ=255)
                                                      cTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGc                                                                 >  2:1150630/1‑150 (MQ=255)
                                                          aCAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGAt                                                             <  1:1554153/150‑1 (MQ=255)
                                                              gaggGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCgag                                                         <  2:2247360/150‑1 (MQ=255)
                                                                 ggACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGagaag                                                      <  2:168348/150‑1 (MQ=255)
                                                                 ggACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGagaag                                                      <  2:348245/150‑1 (MQ=255)
                                                                  gACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGagaaga                                                     >  2:1840295/1‑150 (MQ=255)
                                                                       cTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATg                                                <  1:4306353/150‑1 (MQ=255)
                                                                        tGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGc                                               >  2:775197/1‑150 (MQ=255)
                                                                        tGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGc                                               >  1:1341023/1‑150 (MQ=255)
                                                                          aTAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCag                                             >  1:1334469/1‑150 (MQ=255)
                                                                          aTAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATAAACTAAATCCGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTCTGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGCTGCag                                             >  2:4292720/1‑150 (MQ=255)
                                                                          aTAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATAAACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCag                                             >  1:1686290/1‑150 (MQ=255)
                                                                           tAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCaga                                            <  2:3185740/150‑1 (MQ=255)
                                                                            aGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCagag                                           <  2:1404978/150‑1 (MQ=255)
                                                                             taaaCAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAgg                                          <  1:1486566/149‑1 (MQ=255)
                                                                             gAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAgg                                          <  2:1029705/150‑1 (MQ=255)
                                                                                  aGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAatat                                     <  2:2132734/150‑1 (MQ=255)
                                                                                   gAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATc                                    <  2:909503/150‑1 (MQ=255)
                                                                                     aGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTg                                  <  2:817521/150‑1 (MQ=255)
                                                                                     aGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTg                                  >  1:2742207/1‑150 (MQ=255)
                                                                                        cACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGaa                               <  2:875677/150‑1 (MQ=255)
                                                                                        cACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGCTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGaa                               <  2:330649/150‑1 (MQ=255)
                                                                                         cactgAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAg                              <  2:2726592/146‑1 (MQ=255)
                                                                                           tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCa                            <  1:1907428/150‑1 (MQ=255)
                                                                                             gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAAt                          >  1:3795362/1‑150 (MQ=255)
                                                                                               gCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCa                        <  1:1890425/150‑1 (MQ=255)
                                                                                               gCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCa                        >  1:1714674/1‑150 (MQ=255)
                                                                                                cACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCAt                       >  1:162016/1‑150 (MQ=255)
                                                                                                cACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCAt                       <  2:1241714/150‑1 (MQ=255)
                                                                                                cACCTCAAAAACACCATCAGACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCAt                       <  2:3800495/150‑1 (MQ=255)
                                                                                                 aCCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATg                      <  1:3878482/150‑1 (MQ=255)
                                                                                                 aCCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATg                      <  1:3150671/150‑1 (MQ=255)
                                                                                                 aCCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATg                      <  2:3877021/150‑1 (MQ=255)
                                                                                                      aaacaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGc                 <  2:2457850/150‑1 (MQ=255)
                                                                                                      aaaaaCACCATCATAAACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGc                 >  2:3741393/1‑150 (MQ=255)
                                                                                                        aaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTc               >  2:4140046/1‑150 (MQ=255)
                                                                                                        aaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTc               >  2:4103077/1‑150 (MQ=255)
                                                                                                        aaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTc               >  2:2457215/1‑150 (MQ=255)
                                                                                                        aaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGAGAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTc               >  1:1101947/1‑150 (MQ=255)
                                                                                                         aaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCg              <  2:4062351/150‑1 (MQ=255)
                                                                                                         aaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCg              >  2:210655/1‑150 (MQ=255)
                                                                                                         aaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCg              >  2:3542471/1‑150 (MQ=255)
                                                                                                             ccATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGtt          >  1:2053664/1‑150 (MQ=255)
                                                                                                             ccATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGtt          >  1:2456434/1‑150 (MQ=255)
                                                                                                             ccATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGtt          >  2:3344207/1‑150 (MQ=255)
                                                                                                             ccATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGtt          >  2:681151/1‑150 (MQ=255)
                                                                                                             ccATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGtt          >  2:1169613/1‑150 (MQ=255)
                                                                                                              catcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTa         <  2:1114307/150‑1 (MQ=255)
                                                                                                              catcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTa         >  2:3309399/1‑150 (MQ=255)
                                                                                                              catcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTa         >  1:4001243/1‑150 (MQ=255)
                                                                                                               atcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTAc        >  1:13721/1‑150 (MQ=255)
                                                                                                               atcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTAc        <  2:625216/150‑1 (MQ=255)
                                                                                                                tcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTAcc       >  2:2863336/1‑150 (MQ=255)
                                                                                                                 catACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTACCg      >  1:2085827/1‑150 (MQ=255)
                                                                                                                 catACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTACCg      >  2:2191786/1‑150 (MQ=255)
                                                                                                                  atACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTACCGt     >  2:1935231/1‑150 (MQ=255)
                                                                                                                   tACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTACCGtt    <  2:2533810/150‑1 (MQ=255)
                                                                                                                    acacTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTACCGTTg   >  2:1475647/1‑150 (MQ=255)
                                                                                                                    acacTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTACCGTTg   <  1:937007/150‑1 (MQ=255)
                                                                                                                     cacTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCACGTGTTACCGTTGc  >  1:422548/1‑150 (MQ=255)
                                                                                                                     |                                                                                                                                                     
GATAGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCAGAGGAATATCTGGAAGCAATCATGGAGGCTCGTGTTACCGTTGC  >  CP009273/286967‑287233

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: