Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,429,000 C→A 15.7% intergenic (‑519/+215) ynaE ← / ← ttcC cold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,429,0000CA15.7% 784.0 / 94.3 253intergenic (‑519/+215)ynaE/ttcCcold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (85/128);  new base A (23/17);  total (108/145)
Fisher's exact test for biased strand distribution p-value = 5.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

ACACGGATTGGGTTTAAAAGTATACATAGATAACAGTTTTTATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATG  >  CP009273/1428851‑1429149
                                                                                                                                                     |                                                                                                                                                     
tttaagtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCAcgtc                                                                                                                                                       <  1:3745296‑M1/107‑1 (MQ=255)
tttaagtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCAcgtc                                                                                                                                                       <  1:3475400‑M1/107‑1 (MQ=255)
tttaagtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCAcgtc                                                                                                                                                       <  1:2929730‑M1/107‑1 (MQ=255)
tttaagtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCAcgtc                                                                                                                                                       <  2:1920059‑M1/107‑1 (MQ=255)
 ttaagtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCAcgtcg                                                                                                                                                      >  2:1949266‑M1/43‑150 (MQ=255)
  taagtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCAcgtcgt                                                                                                                                                     <  1:1068329‑M1/109‑1 (MQ=255)
  taagtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCAcgtcgt                                                                                                                                                     <  2:3733295‑M1/109‑1 (MQ=255)
    agtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTAt                                                                                                                                                   <  1:1630650‑M1/111‑1 (MQ=255)
     gtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATg                                                                                                                                                  >  2:1233881‑M1/39‑150 (MQ=255)
     gtgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATg                                                                                                                                                  <  1:716645‑M1/112‑1 (MQ=255)
           ttggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGaa                                                                                                                                            <  1:3589719‑M1/118‑1 (MQ=255)
              gttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGt                                                                                                                                         <  1:2581304‑M1/121‑1 (MQ=255)
               ttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGtg                                                                                                                                        >  2:2864807‑M1/29‑150 (MQ=255)
               ttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGtg                                                                                                                                        >  2:1877770‑M1/29‑150 (MQ=255)
                   ctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGtgctgc                                                                                                                                    <  1:4826758‑M1/126‑1 (MQ=255)
                   ctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGtgctgc                                                                                                                                    <  1:4267284‑M1/126‑1 (MQ=255)
                    tggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCg                                                                                                                                   >  2:4783626‑M1/24‑150 (MQ=255)
                     ggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCgg                                                                                                                                  <  1:2205063‑M1/128‑1 (MQ=255)
                     ggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCgg                                                                                                                                  <  1:530801‑M1/128‑1 (MQ=255)
                       agattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCggtt                                                                                                                                <  2:3369537‑M1/130‑1 (MQ=255)
                           tcagggtgccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCggttggtt                                                                                                                            <  2:534417‑M1/134‑1 (MQ=255)
                             agggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTg                                                                                                                          <  1:4521567‑M1/136‑1 (MQ=255)
                             agggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTg                                                                                                                          >  2:2242367‑M1/15‑150 (MQ=255)
                             agggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTg                                                                                                                          <  1:3031368‑M1/136‑1 (MQ=255)
                                gggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACt                                                                                                                       >  1:1012043‑M1/12‑150 (MQ=255)
                                   ccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttg                                                                                                                    >  2:418035‑M1/9‑150 (MQ=255)
                                   ccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttg                                                                                                                    <  2:1784672‑M1/142‑1 (MQ=255)
                                     agtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttgat                                                                                                                  >  2:2597471‑M1/7‑150 (MQ=255)
                                       ttATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttgattg                                                                                                                <  2:2542842/150‑1 (MQ=255)
                                       tctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGGATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttgattg                                                                                                                >  2:1993368‑M1/5‑150 (MQ=255)
                                       tctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttgattg                                                                                                                >  1:4119054‑M1/5‑150 (MQ=255)
                                       tctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttgattg                                                                                                                >  1:4119172‑M1/5‑150 (MQ=255)
                                          aCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTgaga                                                                                                             >  2:4694920‑M1/2‑150 (MQ=255)
                                          aCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTgaga                                                                                                             >  2:1458884‑M1/2‑150 (MQ=255)
                                           cTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGAc                                                                                                            <  1:4868318/150‑1 (MQ=255)
                                            tGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACg                                                                                                           >  2:1586558/1‑150 (MQ=255)
                                            tGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACg                                                                                                           <  2:3612608/150‑1 (MQ=255)
                                              aagaagaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGtt                                                                                                         >  1:432442/1‑150 (MQ=255)
                                              aagaagaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGtt                                                                                                         <  1:1952697/150‑1 (MQ=255)
                                              aagaagaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGtt                                                                                                         >  2:3404493/1‑150 (MQ=255)
                                                gaagaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGtttt                                                                                                       >  1:568497/1‑150 (MQ=255)
                                                       aaTATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAttt                                                                                                <  2:3254869/150‑1 (MQ=255)
                                                         tatCAAGGTTATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAttttt                                                                                              <  2:678201/150‑1 (MQ=255)
                                                         tatCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCTGTTGGTTTGACTTTGATTGAGACGTTTTGGAAttttt                                                                                              <  1:467743/150‑1 (MQ=255)
                                                          atCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAtttttt                                                                                             <  1:4495363/150‑1 (MQ=255)
                                                             aaGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCg                                                                                          <  1:2596836/150‑1 (MQ=255)
                                                               ggTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGt                                                                                        <  1:1393063/150‑1 (MQ=255)
                                                                gTGATAGAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTg                                                                                       <  2:1348855/150‑1 (MQ=255)
                                                                  gATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGc                                                                                     <  2:3033227/150‑1 (MQ=255)
                                                                    tataGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaa                                                                                   >  2:1234702/1‑150 (MQ=255)
                                                                    tataGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaa                                                                                   >  1:3261066/1‑150 (MQ=255)
                                                                    tataGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTGTGGAATTTTTTTCGGTGGCaa                                                                                   >  2:4499977/1‑150 (MQ=255)
                                                                    tataGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaa                                                                                   >  2:1370450/1‑150 (MQ=255)
                                                                      taGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaaaa                                                                                 <  2:3333908/150‑1 (MQ=255)
                                                                       aGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaaaaa                                                                                <  2:4393545/150‑1 (MQ=255)
                                                                       aGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaaaaa                                                                                >  1:335285/1‑150 (MQ=255)
                                                                        gCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAAt                                                                               <  1:198332/150‑1 (MQ=255)
                                                                          cTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGACTTTTTTTCGGTGGCAAAAATgg                                                                             >  2:4969516/1‑150 (MQ=255)
                                                                           tataTGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATggg                                                                            <  2:5057702/150‑1 (MQ=255)
                                                                            atatGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATgggg                                                                           >  1:628309/1‑150 (MQ=255)
                                                                              atGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCa                                                                         >  2:1023956/1‑150 (MQ=255)
                                                                              atGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCa                                                                         >  2:2424078/1‑150 (MQ=255)
                                                                                  gtttGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAAc                                                                     >  1:3490752/2‑150 (MQ=255)
                                                                                  cTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAAc                                                                     <  1:5060093/150‑1 (MQ=255)
                                                                                      gATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTg                                                                 <  2:2894899/150‑1 (MQ=255)
                                                                                        tGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCa                                                               >  1:1777209/1‑150 (MQ=255)
                                                                                          ctgaggaATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaa                                                             <  2:5097672/148‑1 (MQ=255)
                                                                                             aggaATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaagg                                                          <  1:2652949/150‑1 (MQ=255)
                                                                                              ggaATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggg                                                         >  1:3153337/1‑150 (MQ=255)
                                                                                              ggaATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggg                                                         >  1:3903752/1‑150 (MQ=255)
                                                                                               gaatgaatGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaagggg                                                        >  1:1811051/1‑150 (MQ=255)
                                                                                               gaatgaatGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaagggg                                                        >  1:2976731/1‑150 (MQ=255)
                                                                                                 atgaatgTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggca                                                      <  2:977101/150‑1 (MQ=255)
                                                                                                     atgTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaa                                                  <  2:4072895/150‑1 (MQ=255)
                                                                                                      tgtgATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaa                                                 <  1:3544610/150‑1 (MQ=255)
                                                                                                          aTGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaagggg                                             <  2:1462889/150‑1 (MQ=255)
                                                                                                           tGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaaggggc                                            >  1:1518362/1‑150 (MQ=255)
                                                                                                             ggAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaaggggcaa                                          >  2:2712731/1‑150 (MQ=255)
                                                                                                                   aTGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAag                                    <  2:3768114/150‑1 (MQ=255)
                                                                                                                     gTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAagag                                  >  1:1577239/1‑150 (MQ=255)
                                                                                                                       aTCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTg                                <  1:799098/150‑1 (MQ=255)
                                                                                                                       aTCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTg                                <  1:1198429/150‑1 (MQ=255)
                                                                                                                         cTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGa                              <  2:606384/150‑1 (MQ=255)
                                                                                                                         cTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGa                              <  2:523482/150‑1 (MQ=255)
                                                                                                                           gAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGAtt                            <  1:2401733/150‑1 (MQ=255)
                                                                                                                            aaTAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTa                           <  2:1548935/150‑1 (MQ=255)
                                                                                                                              tAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAATGGGGCACAAAAAGAGTGGATTATc                         >  1:1768413/1‑150 (MQ=255)
                                                                                                                              tAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATc                         <  1:1537878/150‑1 (MQ=255)
                                                                                                                              tAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATc                         <  2:755750/150‑1 (MQ=255)
                                                                                                                              tAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCGGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATc                         >  1:704374/1‑150 (MQ=255)
                                                                                                                               aGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCg                        <  1:1876052/150‑1 (MQ=255)
                                                                                                                                 ttGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTa                      >  2:4300647/1‑150 (MQ=255)
                                                                                                                                        aCCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTAtt               <  1:4585982/150‑1 (MQ=255)
                                                                                                                                         ccGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTAttg              <  1:4021083/150‑1 (MQ=255)
                                                                                                                                            cAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTAttgttg           >  2:1794536/1‑150 (MQ=255)
                                                                                                                                             aGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATtgttgt          >  1:859902/1‑150 (MQ=255)
                                                                                                                                               tCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCg        <  2:90543/150‑1 (MQ=255)
                                                                                                                                                   gtcgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtga    >  2:3905509/1‑150 (MQ=255)
                                                                                                                                                    tcgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCACAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtgat   <  2:3475400/150‑1 (MQ=255)
                                                                                                                                                     cgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtgatg  >  1:931159/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACACGGATTGGGTTTAAAAGTATACATAGATAACAGTTTTTATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATG  >  CP009273/1428851‑1429149

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: