Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 260,931 T→A 5.4% intergenic (‑14/+84) yafX ← / ← ykfF CP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273260,9310TA 5.4% 1145.2 / 4.9 381intergenic (‑14/+84)yafX/ykfFCP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (188/168);  new base A (2/23);  total (190/191)
Fisher's exact test for biased strand distribution p-value = 1.07e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 383 total aligned reads displayed.

GGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCG  >  CP009273/260782‑261078
                                                                                                                                                     |                                                                                                                                                   
ggAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCATATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTTTGTTGTCATGGGTATGTCTCCGt                                                                                                                                                     <  1:847628/150‑1 (MQ=255)
 gAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTc                                                                                                                                                    <  2:271242/150‑1 (MQ=255)
    tAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAAt                                                                                                                                                 <  2:2373183/150‑1 (MQ=255)
    tAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAAt                                                                                                                                                 <  1:2572489/150‑1 (MQ=255)
    tAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGACAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAAt                                                                                                                                                 >  1:4687130/1‑150 (MQ=255)
     aGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATa                                                                                                                                                <  1:1587771/150‑1 (MQ=255)
     aGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATa                                                                                                                                                >  2:3826474/1‑150 (MQ=255)
      gCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaa                                                                                                                                               >  1:3982751/1‑150 (MQ=255)
      gCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaa                                                                                                                                               >  1:4214482/1‑150 (MQ=255)
         aCTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaaaaa                                                                                                                                            <  1:2485509/150‑1 (MQ=255)
          cTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAAc                                                                                                                                           <  1:3049852/150‑1 (MQ=255)
          cTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAAc                                                                                                                                           <  1:2682803/150‑1 (MQ=255)
           tGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACg                                                                                                                                          >  1:4621234/1‑150 (MQ=255)
             cGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGcc                                                                                                                                        >  2:795980/1‑150 (MQ=255)
              gggATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCa                                                                                                                                       <  1:793488/150‑1 (MQ=255)
                 ttAGCGCAAATGTAGTGTGGCCAGAAATCGTTACGCTGTTCATGGGGTACGGGCGTGACGGTCTGTTCAATTTGGGTTTGTTTAGGGGGTGCGAGGTCGTGCTGGGTCTGTGTTGTCTTGGGTATGTCTCCGTCAATAAAAACGCCAgcg                                                                                                                                    <  2:2859378/149‑1 (MQ=255)
                  tAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAgcgg                                                                                                                                   >  2:2521927/1‑150 (MQ=255)
                      gcCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAt                                                                                                                               >  2:226820/1‑150 (MQ=255)
                       ccAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACTCCAGCGGCGAtg                                                                                                                              >  2:3219005/1‑150 (MQ=255)
                        cAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtgg                                                                                                                             >  2:2888330/1‑150 (MQ=255)
                        aaatgTACTGGGGAAACAAATCGATACGGTGTTCTTCGGGTACGGGCTTGACGTTATGTTAAATTTGGGGTTGTTTTGGGGGTGCTTGGTCGTGGTGGGTGTGTGTTTACTTGTTTTTGAACCCGAAAAAAAAAACGCCAGCGGCGAtgg                                                                                                                             <  2:2808027/146‑1 (MQ=255)
                          atgTACTGTGGCAATAATTCTATCCGCTTTTCTTGGGGGACGGGGGTGAGGGTCTGTTAAATTTGGGTTTGGTTTGCGGGTGGGTGGTGGTGCTGGGTTTGTTTTGTCTTGGTTTTTTCTCCTTCAATAAAAACGCCAGCGGCGAtggct                                                                                                                           <  2:64315/148‑1 (MQ=255)
                               cTGGGGCCAGAAATCGATCCGCTGTGCATGGGGGACGGGCGTGCCGGTCAGTTCAAATTCGGGTGGGTTAGCGGGTGCGAGGTCGTGCGGCGTCTGTGTTGACAGGGTTATGTCTCCGTCAATAAAAACGCCCGCGGCGATGGCTGGCGt                                                                                                                      <  2:3718131/150‑1 (MQ=255)
                                  tgGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTTGCGGGGGCGTGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCCCCTACAATAAAAACGCCAGCGGCGATGGCTGGCGTATg                                                                                                                   <  2:4508576/150‑1 (MQ=25)
                                  tgGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATg                                                                                                                   >  2:1899830/1‑150 (MQ=255)
                                    gCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATggg                                                                                                                 >  1:3171059/1‑150 (MQ=255)
                                         aaaTCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGatat                                                                                                            >  2:89627/1‑150 (MQ=255)
                                         aaaTCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGGTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGatat                                                                                                            >  2:3949776/1‑150 (MQ=255)
                                          aaTCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGAtata                                                                                                           >  1:4660775/1‑150 (MQ=255)
                                              gATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgt                                                                                                       >  1:147928/1‑150 (MQ=255)
                                              gATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgt                                                                                                       >  1:2306489/1‑150 (MQ=255)
                                                   ggtGTTCATGGGGGACCGGCGTGAGGGTCTGTACAAATGCGGGTTGGTTAGCGGGTGCGTGGTGGTGCTGCGTCTTTGTTGTCTTGGGTTTGTCTCCTAAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGtt                                                                                                  <  1:1066679/148‑1 (MQ=21)
                                                    cTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTc                                                                                                 >  2:2063150/1‑150 (MQ=255)
                                                     tGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCg                                                                                                >  1:1445907/1‑150 (MQ=255)
                                                     tGTTCAGGGGGGACGGGCGTGCCGGTTAGTTCAATTTGGGGTTGGTTTGGGGGTGCGTGGTGGTGCTGCGTCTGTGTTGTCTTGGGTAGGTCTCCGACAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCg                                                                                                <  2:442622/150‑1 (MQ=25)
                                                      gTTCATGGGGGACGGGCGTGACGGTCTTTTCAATTTGGGGTTGGTTAGCGGGTGCGTGGTGGTGCTGGTTCTGTTTTGTCAGGGGTTTGTCCCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCgg                                                                                               <  1:1088321/150‑1 (MQ=21)
                                                      gTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCgg                                                                                               >  1:2471637/1‑150 (MQ=255)
                                                       ttCATGGGGGACCGGGGTGACGGTCAGTTCAATTTGGGGTTGTTTTGGGGGTGCGTGGTCGTGCTGCTTCGGTGTTTTAATGGGTATGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCggg                                                                                              <  1:1938599/150‑1 (MQ=25)
                                                        tCAGCGGGGACCGGCGTGAGGGTCAGTTCAAATTCGGGTGGGTTAGCGGGGGCGTGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCGGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCgggg                                                                                             <  2:1763146/150‑1 (MQ=25)
                                                         cATGGGGGACGGGCGTGCCGGTCTTTACAATTTCGGGTTGGTTTGCGGGTGCGAGGTCGTGGGGCGTCTGTTTTGTCTTGGGTATGTCTCCTTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGa                                                                                            <  2:2000582/150‑1 (MQ=21)
                                                           tCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAgg                                                                                          >  1:3802763/1‑150 (MQ=255)
                                                            cGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGt                                                                                         >  1:2272120/1‑150 (MQ=255)
                                                               ggACCGGGGTGACGGTATGTTCAAATTGGGGTTGGTTAGGGGGTGCGTGGTGGTGCTGCGTCTGTTTTGTCTGGGTTATGTCTCCGTCAAAAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGaa                                                                                      <  1:1842145/150‑1 (MQ=21)
                                                                 aCCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATg                                                                                    <  2:1930269/150‑1 (MQ=255)
                                                                       gTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGCGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTa                                                                              <  2:120413/150‑1 (MQ=255)
                                                                        tGACGGTCAGTTCAAATTCGGGTTGGTTCGGGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTCTTGGGTATGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTaa                                                                             <  1:396439/150‑1 (MQ=25)
                                                                             gTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTAGGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCt                                                                        >  2:523418/1‑150 (MQ=255)
                                                                              tCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGCTATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCtt                                                                       >  2:717223/1‑150 (MQ=255)
                                                                              tCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCtt                                                                       <  2:1109419/150‑1 (MQ=255)
                                                                              tCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCtt                                                                       <  2:1836516/150‑1 (MQ=255)
                                                                               ctgTTCAAAGTGGGGTTGGTTCGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGTTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTc                                                                      <  1:3902256/148‑1 (MQ=255)
                                                                               cAGTTCAACTTCGGGGTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTGGTCAGGGGTAGGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTc                                                                      <  2:171410/150‑1 (MQ=25)
                                                                                    cAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGAt                                                                 >  1:933554/1‑150 (MQ=255)
                                                                                    cAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGAt                                                                 >  1:3876857/1‑150 (MQ=255)
                                                                                     aatttCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGACATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATc                                                                <  2:3359364/147‑1 (MQ=255)
                                                                                       aTTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAg                                                              >  2:459544/1‑150 (MQ=255)
                                                                                        ttCGGGTTGGTTTGGGGGTGCGTGGTGGTGCTGGGTCTGTGTTGTCATGGGTTTGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGc                                                             <  1:2224935/150‑1 (MQ=25)
                                                                                          cGGGTTGGTTAGCGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTCTGGGGTTTGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCgg                                                           <  1:362805/150‑1 (MQ=37)
                                                                                             gttggttAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTg                                                        >  1:3873004/1‑150 (MQ=255)
                                                                                               tggttAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGc                                                      >  1:3143503/1‑150 (MQ=255)
                                                                                                  ttcgggggtgggtGGTGGGGCGGCGTGTGTGTTGTCTGGGGTATGTCCCCGACAAAAAACACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCgt                                                   <  1:2031894/145‑1 (MQ=25)
                                                                                                   tAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCgtg                                                  >  2:1585876/1‑150 (MQ=255)
                                                                                                     gggggTGCGTGGTCGTGCTGCGTCTGTGTTGTCTTGGTTATGTCTCCTTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTc                                                <  1:2544043/148‑1 (MQ=37)
                                                                                                      cGGGTGCGTGGTGGTGCGGGGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCt                                               <  1:3426236/150‑1 (MQ=255)
                                                                                                          tGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTAc                                           <  2:3458154/150‑1 (MQ=255)
                                                                                                           gCGTGGTGGTGTTGGGTCTGTGTTGTCTTGGGTTTGGGTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACg                                          <  1:2888330/150‑1 (MQ=37)
                                                                                                            cGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACgg                                         >  2:257944/1‑150 (MQ=255)
                                                                                                              aGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGAt                                       <  2:2138087/150‑1 (MQ=255)
                                                                                                               gggggtggtgGGTGTGTGTGGTCCTGGGTTGGTCCCCGTAAATAAACACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATg                                      <  1:1246689/146‑1 (MQ=37)
                                                                                                                gTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGc                                     <  2:3000458/150‑1 (MQ=255)
                                                                                                                 tggTGCTGGTTCTGTGTGTTCATGGGTATGTCTCCGAAAAAAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGcc                                    <  1:656233/148‑1 (MQ=37)
                                                                                                                 tggTGCTGCGTCTGTGTTTGCTGGGGTTTGACTCCGACAATAAACACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGcc                                    <  1:2111732/148‑1 (MQ=38)
                                                                                                                   gTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTg                                  <  2:3018261/150‑1 (MQ=255)
                                                                                                                    tgctgcGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGa                                 >  2:4735441/1‑150 (MQ=255)
                                                                                                                    tgctgcGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGa                                 <  1:3474072/150‑1 (MQ=255)
                                                                                                                    tgctgcGTCTGTGTGGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGa                                 <  2:115662/150‑1 (MQ=255)
                                                                                                                     gctgcGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAg                                <  2:1483506/150‑1 (MQ=255)
                                                                                                                     gctgcGTCTGGGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAg                                <  2:137411/150‑1 (MQ=255)
                                                                                                                        gcGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTg                             <  2:4214482/150‑1 (MQ=255)
                                                                                                                        gcGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTg                             <  2:4087426/150‑1 (MQ=255)
                                                                                                                         cGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGc                            >  1:1783410/1‑150 (MQ=255)
                                                                                                                          gTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCg                           >  1:3327000/1‑150 (MQ=255)
                                                                                                                          gTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCg                           >  1:436361/1‑150 (MQ=255)
                                                                                                                            cTGTTTTGTCTTGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGa                         <  2:1051846/150‑1 (MQ=255)
                                                                                                                             tgtgTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGat                        >  1:3769081/1‑150 (MQ=255)
                                                                                                                                gttgtCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGAtata                     >  1:3571811/1‑150 (MQ=255)
                                                                                                                                gttgtCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGAtata                     >  1:1382737/1‑150 (MQ=255)
                                                                                                                                 ttgtCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAg                    <  2:303880/150‑1 (MQ=255)
                                                                                                                                 ttgtCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAg                    <  2:544545/150‑1 (MQ=255)
                                                                                                                                      aTGGGTAGGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATAGAAAGTGGGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGtt               <  2:3390811/150‑1 (MQ=255)
                                                                                                                                        gggTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTATGGATGCCTGAGGTGCGGATATAGCGGTTaa             >  2:3606954/1‑150 (MQ=255)
                                                                                                                                        gggTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTAGGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTaa             >  2:2210627/1‑150 (MQ=255)
                                                                                                                                           tATGTCTCCTAAAATAAACACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGAc          <  1:1266526/150‑1 (MQ=255)
                                                                                                                                            aTGTCTCCTAAAATAAACACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGAcc         <  1:76846/150‑1 (MQ=255)
                                                                                                                                            aTGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGAcc         >  1:2419990/1‑150 (MQ=255)
                                                                                                                                            aTGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGAcc         >  1:721864/1‑150 (MQ=255)
                                                                                                                                               tctcCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTc      <  2:3802763/150‑1 (MQ=255)
                                                                                                                                                ctcCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCa     >  2:3494560/1‑150 (MQ=255)
                                                                                                                                                 tcCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCAc    <  2:3873358/150‑1 (MQ=255)
                                                                                                                                                  ccGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCAcc   <  1:1652050/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCG  >  CP009273/260782‑261078

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: