Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,757,927 A→T 17.7% S196T (TCC→ACC)  yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,757,9270AT17.7% 878.6 / 163.1 320S196T (TCC→ACC) yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base A (122/141);  new base T (27/30);  total (149/171)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

TCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGT  >  CP009273/2757779‑2758076
                                                                                                                                                    |                                                                                                                                                     
tCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGaa                                                                                                                                                      <  1:3486815/150‑1 (MQ=255)
  aTCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACt                                                                                                                                                    <  1:637453/150‑1 (MQ=255)
     ttCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTgag                                                                                                                                                 <  2:815716/150‑1 (MQ=255)
      tCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTgaga                                                                                                                                                <  1:4527311/150‑1 (MQ=255)
       cGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGaa                                                                                                                                               >  2:3690298/1‑150 (MQ=255)
        gACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAAt                                                                                                                                              >  2:1594467/1‑150 (MQ=255)
        gACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAAt                                                                                                                                              >  1:962811/1‑150 (MQ=255)
         acacGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                                                             <  1:1286969/150‑1 (MQ=255)
            cGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGAcc                                                                                                                                          <  2:4458892/150‑1 (MQ=255)
                gctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCAt                                                                                                                                      >  1:1077014/1‑150 (MQ=255)
                gctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCAt                                                                                                                                      <  1:5058670/150‑1 (MQ=255)
                gctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCAt                                                                                                                                      <  1:4559961/150‑1 (MQ=255)
                 ctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                                                     <  2:2482388/150‑1 (MQ=255)
                    gAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAg                                                                                                                                  <  2:3972742/150‑1 (MQ=255)
                      aCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGAt                                                                                                                                >  2:1739788/1‑150 (MQ=255)
                         aCTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCtt                                                                                                                             <  2:1820357/150‑1 (MQ=255)
                           tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTgg                                                                                                                           >  1:1361765/1‑150 (MQ=255)
                           tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTgg                                                                                                                           <  2:4824322/150‑1 (MQ=255)
                           tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTgg                                                                                                                           <  2:2249053/150‑1 (MQ=255)
                            cAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTggg                                                                                                                          <  1:2337975/150‑1 (MQ=255)
                             aaGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGa                                                                                                                         <  2:416913/150‑1 (MQ=255)
                               gTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAgg                                                                                                                       <  2:3261644/150‑1 (MQ=255)
                                tAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGc                                                                                                                      <  1:28955/150‑1 (MQ=255)
                                     cTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGcc                                                                                                                 >  2:4377671/1‑150 (MQ=255)
                                      tGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCg                                                                                                                >  1:245474/1‑150 (MQ=255)
                                      tGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCg                                                                                                                <  2:4490783/150‑1 (MQ=255)
                                      tGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCg                                                                                                                <  1:4294723/150‑1 (MQ=255)
                                       ggTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGt                                                                                                               >  1:1624281/1‑150 (MQ=255)
                                        gTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTa                                                                                                              >  1:2715464/1‑150 (MQ=255)
                                        gTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTa                                                                                                              <  2:1439108/150‑1 (MQ=255)
                                         tCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTAt                                                                                                             <  2:898125/150‑1 (MQ=255)
                                          cGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATg                                                                                                            <  1:1758559/150‑1 (MQ=255)
                                           ggaaGGAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGt                                                                                                           <  2:4034245/148‑1 (MQ=255)
                                           gTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGt                                                                                                           >  1:1147718/1‑150 (MQ=255)
                                            tAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGtt                                                                                                          >  1:3182110/1‑150 (MQ=255)
                                             aaGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTc                                                                                                         >  1:1008932/1‑150 (MQ=255)
                                               gTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCtt                                                                                                       <  1:3377985/150‑1 (MQ=255)
                                                tAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTg                                                                                                      >  2:2644518/1‑150 (MQ=255)
                                                        gTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAg                                                                                              >  1:2663583/1‑150 (MQ=255)
                                                        gTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAg                                                                                              <  2:4726328/150‑1 (MQ=255)
                                                            tcatcaCTGTGAATGCGTCGAGTATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGATAATGAACCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATACCTAGAAGt                                                                                          <  2:5036088/150‑1 (MQ=255)
                                                            tcaGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATATTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGt                                                                                          <  2:3456103/150‑1 (MQ=255)
                                                             caGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGtt                                                                                         >  1:2609301/1‑150 (MQ=255)
                                                                  cTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGt                                                                                    >  1:1483677/1‑150 (MQ=255)
                                                                     tgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTcg                                                                                 <  1:4437543/150‑1 (MQ=255)
                                                                     tgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTcg                                                                                 >  2:2150511/1‑150 (MQ=255)
                                                                           cGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAgg                                                                           >  2:4552884/1‑150 (MQ=255)
                                                                            gTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAggg                                                                          <  2:3486841/150‑1 (MQ=255)
                                                                             tCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGa                                                                         >  1:4240267/1‑150 (MQ=255)
                                                                              cGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGAt                                                                        >  1:3307978/1‑150 (MQ=255)
                                                                               gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGACCTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATg                                                                       >  2:2055823/1‑150 (MQ=255)
                                                                               gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATg                                                                       <  1:4595844/150‑1 (MQ=255)
                                                                                agaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGa                                                                      <  2:4896158/150‑1 (MQ=255)
                                                                                agaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGa                                                                      >  1:2244854/1‑150 (MQ=255)
                                                                                agaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGa                                                                      <  2:678799/150‑1 (MQ=255)
                                                                                  aaTGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaaa                                                                    >  2:1245781/1‑150 (MQ=255)
                                                                                   aTGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaaaa                                                                   >  1:4068799/1‑150 (MQ=255)
                                                                                     gCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAaca                                                                 >  1:1179553/1‑150 (MQ=255)
                                                                                      cTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATCTAGATCATGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAacac                                                                >  2:2989758/1‑150 (MQ=255)
                                                                                       tCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACt                                                               >  2:5141776/1‑150 (MQ=255)
                                                                                        cGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACtt                                                              <  1:450177/150‑1 (MQ=255)
                                                                                               gcgcGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGt                                                       <  2:1149920/150‑1 (MQ=255)
                                                                                                cgcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTa                                                      <  1:2544558/150‑1 (MQ=255)
                                                                                                cgcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTa                                                      >  2:2789739/1‑150 (MQ=255)
                                                                                                cgcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTa                                                      <  1:3062863/150‑1 (MQ=255)
                                                                                                cgcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTa                                                      >  2:3398273/1‑150 (MQ=255)
                                                                                                cgcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTa                                                      <  1:2152758/150‑1 (MQ=255)
                                                                                                  cgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGc                                                    >  1:1553598/1‑150 (MQ=255)
                                                                                                   gTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGcc                                                   >  1:4679101/1‑150 (MQ=255)
                                                                                                      acaGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                <  2:1797491/150‑1 (MQ=255)
                                                                                                       caGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCAt                                               <  1:3896090/150‑1 (MQ=255)
                                                                                                       caGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCAt                                               <  2:2586173/150‑1 (MQ=255)
                                                                                                        aGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATc                                              >  1:4401485/1‑150 (MQ=255)
                                                                                                        aGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATc                                              >  1:1847728/1‑150 (MQ=255)
                                                                                                          tATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGa                                            >  2:4691970/1‑150 (MQ=255)
                                                                                                           aTCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGAt                                           >  1:3689811/1‑150 (MQ=255)
                                                                                                                ggATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAAt                                      <  2:4342175/150‑1 (MQ=255)
                                                                                                                   tGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTg                                   <  2:4068799/150‑1 (MQ=255)
                                                                                                                      ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaaa                                <  2:105935/150‑1 (MQ=255)
                                                                                                                           cAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCtt                           <  2:693434/150‑1 (MQ=255)
                                                                                                                           cAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTATTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACTCTTCCCATGTAGCCGAATCGATCCAATCTGAAACGCtt                           >  2:3837639/1‑150 (MQ=255)
                                                                                                                            aaaTGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTta                          <  2:5061367/150‑1 (MQ=255)
                                                                                                                             aaTGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTtat                         >  1:4804755/1‑150 (MQ=255)
                                                                                                                              aTGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTtata                        >  1:397729/1‑150 (MQ=255)
                                                                                                                               tGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATAc                       <  1:891164/150‑1 (MQ=255)
                                                                                                                                gTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATAcc                      <  1:2485812/150‑1 (MQ=255)
                                                                                                                                gTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATAcc                      >  1:692523/1‑150 (MQ=255)
                                                                                                                                 tCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCt                     >  2:1812721/1‑150 (MQ=255)
                                                                                                                                  cTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTg                    >  1:1407181/1‑150 (MQ=255)
                                                                                                                                  cTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTg                    <  2:3095569/150‑1 (MQ=255)
                                                                                                                                      ggCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTc                <  2:2815476/150‑1 (MQ=255)
                                                                                                                                         gACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGt             <  2:3662053/150‑1 (MQ=255)
                                                                                                                                          aCGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGtt            <  2:4270803/150‑1 (MQ=255)
                                                                                                                                            ggTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaa          >  1:2070179/1‑150 (MQ=255)
                                                                                                                                            ggTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaa          >  2:1880118/1‑150 (MQ=255)
                                                                                                                                             gTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaac         >  2:2870152/1‑149 (MQ=255)
                                                                                                                                              tttCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaaa        <  2:1407181/150‑1 (MQ=255)
                                                                                                                                              tttCGGAACTGAGTATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaaa        <  1:4846312/150‑1 (MQ=255)
                                                                                                                                                 cGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAAt     >  1:3603298/1‑150 (MQ=255)
                                                                                                                                                  ggAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAAtt    <  1:284302/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGT  >  CP009273/2757779‑2758076

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: