Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,633,239 A→G 16.9% intergenic (‑315/+48) cspI ← / ← rzpQ Qin prophage; cold shock protein/Rz‑like protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,633,2390AG16.9% 371.8 / 61.8 148intergenic (‑315/+48)cspI/rzpQQin prophage; cold shock protein/Rz‑like protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base A (71/51);  new base G (16/8);  total (89/59)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

ACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTT  >  CP009273/1633091‑1633387
                                                                                                                                                    |                                                                                                                                                    
accaccCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTtttat                                                                                                                                                     <  1:1997422/150‑1 (MQ=255)
  caccCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTtttattt                                                                                                                                                   >  2:354269/1‑150 (MQ=255)
  caccCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTtttattt                                                                                                                                                   <  1:663387/150‑1 (MQ=255)
      cATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAAt                                                                                                                                               <  2:5298/150‑1 (MQ=255)
       aTAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAAtt                                                                                                                                              <  2:1886682/150‑1 (MQ=255)
       aTAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAAtt                                                                                                                                              <  2:46626/150‑1 (MQ=255)
        tAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAAttt                                                                                                                                             <  2:4445784/150‑1 (MQ=255)
           cAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATATCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTAc                                                                                                                                          <  2:2485233/150‑1 (MQ=255)
           cAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTAc                                                                                                                                          >  1:287883/1‑150 (MQ=255)
           cAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTAc                                                                                                                                          <  2:58438/150‑1 (MQ=255)
            aaGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGTCTCATGGAATAGTCCTTTTtatttttatttttaa                                                                                                                                         >  2:2449863/1‑147 (MQ=255)
            aaGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACa                                                                                                                                         >  1:3047950/1‑150 (MQ=255)
             aGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACAt                                                                                                                                        <  1:4257917/150‑1 (MQ=255)
             aGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTAATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTTtttaaattttaaat                                                                                                                                        >  2:1765506/1‑150 (MQ=255)
                aTTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTcc                                                                                                                                     >  1:3050269/1‑150 (MQ=255)
                aTTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTcc                                                                                                                                     >  1:3050249/1‑150 (MQ=255)
                 ttACTTTTGTAAAGTATAACAAAGCTATGAAAGTATATCTTTATTTATTAATAAAAAGTGTCATTCTGATTAAGACCTTTTATCTTACCATTAAGATTTCAGGTATTTTGGATCATTGAAGAGTCCTTTTTATTTAAATTTTACATTcct                                                                                                                                    <  2:4552034/150‑2 (MQ=255)
                     ttttGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGAt                                                                                                                                <  1:1626469/150‑1 (MQ=255)
                       ttGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGt                                                                                                                              >  2:3093814/1‑150 (MQ=255)
                       ttGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTGTGTATTTAAATTTTACATTCCGCGATGt                                                                                                                              >  1:2459625/1‑150 (MQ=255)
                          tAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTaaa                                                                                                                           <  2:4792601/150‑1 (MQ=255)
                          tAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTaaa                                                                                                                           >  1:756966/1‑150 (MQ=255)
                           aaaGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAAt                                                                                                                          <  1:3580123/150‑1 (MQ=255)
                                aaaaCACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTcc                                                                                                                     <  1:1369559/150‑1 (MQ=255)
                                 aaaCACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCg                                                                                                                    <  2:4336822/150‑1 (MQ=255)
                                     acaGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGAttt                                                                                                                <  1:616999/150‑1 (MQ=255)
                                     acaGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGAttt                                                                                                                <  1:2238698/150‑1 (MQ=255)
                                     acaGCAATGAAAGAAGAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGAttt                                                                                                                <  1:616933/150‑1 (MQ=255)
                                        gCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAat                                                                                                             <  1:1745794/150‑1 (MQ=255)
                                        gCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAat                                                                                                             <  1:1745763/150‑1 (MQ=255)
                                        gCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAat                                                                                                             >  2:1758176/1‑150 (MQ=255)
                                        gCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACATTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAat                                                                                                             <  1:2782503/150‑1 (MQ=255)
                                         cAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAata                                                                                                            >  1:4433163/1‑150 (MQ=255)
                                          aaTGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAatat                                                                                                           >  1:756681/1‑150 (MQ=255)
                                           aTGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATAtt                                                                                                          >  1:2272248/1‑150 (MQ=255)
                                            tGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTa                                                                                                         <  2:834599/150‑1 (MQ=255)
                                             gAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTAc                                                                                                        >  2:4233144/1‑150 (MQ=255)
                                              aaaGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATGTAATATTAcc                                                                                                       >  1:4137569/1‑150 (MQ=255)
                                               aaGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTAccc                                                                                                      >  2:3871161/1‑150 (MQ=255)
                                                 gAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACAATCCGCGAGGTAAATGTTCCCATTTAATATTAcccaa                                                                                                    >  1:84799/1‑148 (MQ=255)
                                                  aaTAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTAc                                                                                                   >  1:4365975/1‑150 (MQ=255)
                                                    tAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACAt                                                                                                 >  1:2600287/1‑150 (MQ=255)
                                                    tAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACAt                                                                                                 >  1:3181642/1‑150 (MQ=255)
                                                      gCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACAttt                                                                                               <  2:3206690/150‑1 (MQ=255)
                                                       cTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTg                                                                                              >  2:680401/1‑150 (MQ=255)
                                                       cTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTg                                                                                              >  2:3564702/1‑150 (MQ=255)
                                                         ttatttatTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGAt                                                                                            <  1:2052665/150‑1 (MQ=255)
                                                          tatttattAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATg                                                                                           >  1:3060279/1‑150 (MQ=255)
                                                           atttattAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGc                                                                                          >  2:4305048/1‑150 (MQ=255)
                                                            tttattAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTAAATTACGCGATGTAAATTTTGCGATTTAATATTAACCTACATTTGTTGCt                                                                                         >  1:418664/1‑150 (MQ=255)
                                                               attaATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTGCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCtttt                                                                                      >  2:858890/1‑150 (MQ=255)
                                                                 taataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATTCTTTTTa                                                                                    >  1:4662996/1‑150 (MQ=255)
                                                                 taataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTAAACTACATTTTATGCtttttt                                                                                    >  1:2727409/1‑149 (MQ=255)
                                                                 caataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTa                                                                                    >  2:4142741/2‑150 (MQ=255)
                                                                 caataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTa                                                                                    >  2:4142816/2‑150 (MQ=255)
                                                                  aataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTAt                                                                                   <  1:579615/150‑1 (MQ=255)
                                                                   ataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTAtc                                                                                  <  1:2251591/150‑1 (MQ=255)
                                                                   ataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTAtc                                                                                  <  1:2251693/150‑1 (MQ=255)
                                                                      aaaCGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTTCCCTACATTTGATGATTTTTATctct                                                                               >  1:4131927/1‑150 (MQ=255)
                                                                      aaaCGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTACGATTTAATATTACACTACATTTGATGCTTTTTATatct                                                                               >  2:4272509/1‑150 (MQ=255)
                                                                       aaCGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTAACCTACATTTGATTCTTTTTATCTCtt                                                                              >  2:4285105/1‑150 (MQ=255)
                                                                        aCGTGTCATTCTGATTAAGGCCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTa                                                                             >  2:4796392/1‑150 (MQ=255)
                                                                         cGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTaa                                                                            <  2:1200368/150‑1 (MQ=255)
                                                                                ttCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCa                                                                     <  1:4701809/150‑1 (MQ=255)
                                                                                ttCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCa                                                                     <  1:4701336/150‑1 (MQ=255)
                                                                                  atGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTACGACTTAATATTACACTACATTTGATGCTTTTTATCTCTTAAAGATTCATa                                                                   >  2:2203578/2‑150 (MQ=255)
                                                                                   tGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAg                                                                  >  1:348318/1‑150 (MQ=255)
                                                                                   tGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTAATAg                                                                  >  2:136042/1‑150 (MQ=255)
                                                                                    gATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAATATTCATAGa                                                                 >  1:660639/1‑150 (MQ=255)
                                                                                    gATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGa                                                                 <  1:4197955/150‑1 (MQ=255)
                                                                                    gATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGa                                                                 >  2:856975/1‑150 (MQ=255)
                                                                                     aatAAGACCCTTTAACTTAACCTTAAGATTTAACGAAATTTTGCACATTGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGAt                                                                <  2:4313961/148‑1 (MQ=255)
                                                                                     aTTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGAt                                                                <  2:3047950/150‑1 (MQ=255)
                                                                                               tttATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTAATTAAATTTTACATACCGCGATGTAAATGTTCCGAATTAATATTACCCTACAATTGATGCTTTTTATCACTTAAAGATTCATAGATCTGTTGACaa                                                      >  2:997816/1‑150 (MQ=255)
                                                                                                 tATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGt                                                    >  1:297268/1‑150 (MQ=255)
                                                                                                   tCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCa                                                  <  1:273063/150‑1 (MQ=255)
                                                                                                    cTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATTTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCAc                                                 >  1:1077234/1‑150 (MQ=255)
                                                                                                    cTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCAc                                                 >  1:1078068/1‑150 (MQ=255)
                                                                                                    cTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCAc                                                 >  2:1420691/1‑150 (MQ=255)
                                                                                                     ttACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACt                                                <  2:1884688/150‑1 (MQ=255)
                                                                                                        cccTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCt                                             <  1:1844814/150‑1 (MQ=255)
                                                                                                            tAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATATGATGCTTTTTATCTCTTAAAGATAAATAGATCTGTTGACAAGTCACACCTGCGa                                         >  1:2687637/1‑150 (MQ=255)
                                                                                                                     aGGAATTTTGGCTCATGGAAGAGTCCGTGTTATTTTAAATTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAACTCACACCCGCGATGTAGCGtt                                >  1:2863196/1‑150 (MQ=255)
                                                                                                                        aaTTTTGGCTCATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGt                             <  1:3831291/150‑1 (MQ=255)
                                                                                                                          ttttGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCa                           <  1:2135654/150‑1 (MQ=255)
                                                                                                                                  cATGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTc                   <  2:2631154/150‑1 (MQ=255)
                                                                                                                                  cATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTATCGTTCGTCAGCAATTTc                   >  1:3617287/1‑150 (MQ=255)
                                                                                                                                   aTGGAAGAGTCCTTTTTGTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCa                  >  2:3975649/1‑150 (MQ=255)
                                                                                                                                    tGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAg                 <  1:2449863/150‑1 (MQ=255)
                                                                                                                                          agTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATaaa           >  2:4217452/1‑150 (MQ=255)
                                                                                                                                           gtacTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAg          <  2:3346527/147‑1 (MQ=255)
                                                                                                                                           gTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAg          <  2:4696799/150‑1 (MQ=255)
                                                                                                                                            tCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGc         <  1:4350968/150‑1 (MQ=255)
                                                                                                                                             ccTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCt        >  2:724128/1‑150 (MQ=255)
                                                                                                                                                ttttATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAg     >  1:693458/1‑150 (MQ=255)
                                                                                                                                                ttttATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAg     <  2:2108928/150‑1 (MQ=255)
                                                                                                                                                ttttATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAg     <  2:3060279/150‑1 (MQ=255)
                                                                                                                                                  ttatttaAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCt   >  2:3856904/1‑150 (MQ=255)
                                                                                                                                                  ttatttaAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCt   <  1:856975/150‑1 (MQ=255)
                                                                                                                                                   tatttaAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCtt  <  2:756681/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATTTATTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTT  >  CP009273/1633091‑1633387

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