Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 280,433 G→T 9.6% G508C (GGC→TGC)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273280,4330GT 9.6% 1616.2 / 100.2 449G508C (GGC→TGC) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (195/210);  new base T (25/19);  total (220/229)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 475 total aligned reads displayed.

CGTCGGTGGTGGGCGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGCCG  >  CP009273/280284‑280582
                                                                                                                                                     |                                                                                                                                                     
 gTCGGTGGTGGGCGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTccggcaccgg                                                                                                                                                      >  2:3395266/1‑150 (MQ=255)
   cGGTGGTGGGCGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCCTGGTGGTGATCGGCGGCGGGCCGTccggcaccggca                                                                                                                                                    >  2:3323653/1‑150 (MQ=255)
   cGGTGGTGGGCGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTccggcaccggca                                                                                                                                                    >  2:1368839/1‑150 (MQ=255)
    ggtggtGGGCGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCAt                                                                                                                                                   >  1:3627684/1‑150 (MQ=255)
       ggtggGCGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGa                                                                                                                                                >  2:3170443/1‑150 (MQ=255)
         tggGCGAAGATGGCGTATACCACCACCCCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAag                                                                                                                                              <  2:273331/150‑1 (MQ=255)
             cGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGAc                                                                                                                                          >  2:2788555/1‑150 (MQ=255)
             cGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGAc                                                                                                                                          >  2:3889695/1‑150 (MQ=255)
                aGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTa                                                                                                                                       <  1:3456920/150‑1 (MQ=255)
                aGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTa                                                                                                                                       <  2:3739183/150‑1 (MQ=255)
                      cGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCt                                                                                                                                 <  1:1733564/150‑1 (MQ=255)
                      cGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCt                                                                                                                                 <  1:2855063/150‑1 (MQ=255)
                       gTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTc                                                                                                                                <  1:2565689/150‑1 (MQ=255)
                        tataCCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCa                                                                                                                               >  1:1572163/1‑150 (MQ=255)
                        tataCCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCa                                                                                                                               <  2:364730/150‑1 (MQ=255)
                          taCCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCAcc                                                                                                                             >  1:1203288/1‑150 (MQ=255)
                          taCCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCAcc                                                                                                                             <  1:3322263/150‑1 (MQ=255)
                           accaccACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCt                                                                                                                            >  2:2469823/1‑150 (MQ=255)
                             caccacCCCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTcc                                                                                                                          <  2:615238/150‑1 (MQ=255)
                              accacACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCcg                                                                                                                         >  2:1540198/1‑150 (MQ=255)
                              accacACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCcg                                                                                                                         >  1:887950/1‑150 (MQ=255)
                                 acacCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCgcgc                                                                                                                      >  2:2301915/1‑150 (MQ=255)
                                      ggccggGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAg                                                                                                                 <  2:4308397/150‑1 (MQ=255)
                                      ggccggGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAg                                                                                                                 <  2:1924734/150‑1 (MQ=255)
                                       gccggGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGc                                                                                                                <  2:709521/150‑1 (MQ=255)
                                             tgcgggtgTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATAtct                                                                                                          <  2:1247227/150‑1 (MQ=255)
                                               cgggtgTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCg                                                                                                        >  2:2375962/1‑150 (MQ=255)
                                                 ggtgTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTg                                                                                                      >  1:410887/1‑150 (MQ=255)
                                                   tgtTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTggg                                                                                                    <  1:2050219/150‑1 (MQ=255)
                                                   tgtTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTggg                                                                                                    <  1:3245222/150‑1 (MQ=255)
                                                   tgtTTGTATCGGAAGAGCAGGCGATTAAGGAGATGTAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGTCCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTggg                                                                                                    <  1:874129/150‑1 (MQ=255)
                                                     tttGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGc                                                                                                  <  2:4438548/150‑1 (MQ=255)
                                                        gTCTCGGAAGCGCAGTCGATGATGTCGATCATGCGGGAAGAGATTGTGCAGGGCGATATCTTGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGATGAGACCTACCAGCTCACCTCCGCGATAAAGCTTATCTCGTGTGGCAAg                                                                                               <  2:76878/150‑1 (MQ=255)
                                                        gTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAg                                                                                               <  2:2671723/150‑1 (MQ=255)
                                                        gTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAg                                                                                               >  1:4526248/1‑150 (MQ=255)
                                                         tctcGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGa                                                                                              <  2:4144635/150‑1 (MQ=255)
                                                         tctcGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGa                                                                                              <  1:711318/150‑1 (MQ=255)
                                                             ggAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGgt                                                                                          <  2:303551/150‑1 (MQ=255)
                                                             ggAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGgt                                                                                          >  2:1494605/1‑150 (MQ=255)
                                                                  cgcAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGc                                                                                     <  2:1802539/150‑1 (MQ=255)
                                                                   gcAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCt                                                                                    <  1:2236472/150‑1 (MQ=255)
                                                                    cAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCtc                                                                                   <  1:4378584/150‑1 (MQ=255)
                                                                      ggcgatcaaggcgatcaaGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCtcat                                                                                 <  2:2214034/150‑1 (MQ=255)
                                                                          atcaaggcgatcaagCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCAcc                                                                             <  1:2648202/150‑1 (MQ=255)
                                                                          atcaaggcgatcaagCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCAcc                                                                             <  1:1277992/150‑1 (MQ=255)
                                                                             aaggcgatcaagCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGAt                                                                          >  1:530162/1‑150 (MQ=255)
                                                                             aaggcgatcaagCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGCCCCTCGAGAGGCTGTATTGAAGACACCTAGCAGAACACCACACCGCTCAACCCAAGAAAGTGGTGCAAGACGCTCTCGCTCATCACCGAt                                                                          >  1:1268225/1‑150 (MQ=255)
                                                                              aggcgatcaagCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATg                                                                         <  1:4703929/150‑1 (MQ=255)
                                                                              aggcgatcaagCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATg                                                                         >  1:975068/1‑150 (MQ=255)
                                                                                  gatcaagCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATgcgcg                                                                     <  1:3207817/150‑1 (MQ=255)
                                                                                        gCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTtctc                                                               <  1:1389321/150‑1 (MQ=255)
                                                                                         cGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCg                                                              >  2:3496650/1‑150 (MQ=255)
                                                                                             aaGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCg                                                          <  2:1685635/150‑1 (MQ=255)
                                                                                              agagATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCgt                                                         <  2:1214680/150‑1 (MQ=255)
                                                                                              agagATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCgt                                                         >  1:2679116/1‑150 (MQ=255)
                                                                                                   ttGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGa                                                    >  1:4347473/1‑150 (MQ=255)
                                                                                                    tgtgCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGAc                                                   <  1:407495/150‑1 (MQ=255)
                                                                                                    tgtgCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGAc                                                   <  1:2813858/150‑1 (MQ=255)
                                                                                                     gtgCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACg                                                  <  1:3776821/150‑1 (MQ=255)
                                                                                                     gtgCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACg                                                  <  1:3818750/150‑1 (MQ=255)
                                                                                                     gtgCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACg                                                  <  1:1949796/150‑1 (MQ=255)
                                                                                                     gtgCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACg                                                  <  1:1343809/150‑1 (MQ=255)
                                                                                                        cAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGgc                                               <  1:2331711/150‑1 (MQ=255)
                                                                                                        cAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGgc                                               >  1:1296505/1‑150 (MQ=255)
                                                                                                         aGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGgcg                                              <  2:3750749/150‑1 (MQ=255)
                                                                                                            gCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCt                                           <  2:2839297/150‑1 (MQ=255)
                                                                                                            gCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCt                                           <  2:1967884/150‑1 (MQ=255)
                                                                                                             cGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTg                                          <  2:466432/150‑1 (MQ=255)
                                                                                                              gATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGc                                         >  1:3360113/1‑150 (MQ=255)
                                                                                                                tatCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCtt                                       >  1:4319121/1‑150 (MQ=255)
                                                                                                                tatCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCtt                                       <  2:405100/150‑1 (MQ=255)
                                                                                                                tatCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCtt                                       >  1:2422335/1‑150 (MQ=255)
                                                                                                                  tCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCg                                     <  1:651959/150‑1 (MQ=255)
                                                                                                                     tggtggTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGcc                                  <  2:934976/150‑1 (MQ=255)
                                                                                                                       gtggtgATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCAc                                >  2:408112/1‑150 (MQ=255)
                                                                                                                       gtggtgATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCAc                                <  2:4609454/150‑1 (MQ=255)
                                                                                                                        tggtgATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACg                               >  2:2876693/1‑150 (MQ=255)
                                                                                                                        tggtgATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACg                               >  2:730538/1‑150 (MQ=255)
                                                                                                                           tgATCGGCGGCGGGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACgtgt                            <  2:49790/150‑1 (MQ=255)
                                                                                                                           tgATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCTCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACgtgt                            >  1:1000292/1‑150 (MQ=255)
                                                                                                                            gATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTc                           <  2:1336506/150‑1 (MQ=255)
                                                                                                                             aTCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCg                          <  1:2769330/150‑1 (MQ=255)
                                                                                                                               cggcggcggGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGcc                        <  2:368794/150‑1 (MQ=255)
                                                                                                                                ggcggcggGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCg                       >  2:1091234/1‑150 (MQ=255)
                                                                                                                                   ggcggGCCGTCCGGCACCTGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAg                    <  2:291908/150‑1 (MQ=255)
                                                                                                                                   ggcggGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAg                    >  1:1143274/1‑150 (MQ=255)
                                                                                                                                   ggcggGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAg                    >  1:2182596/1‑150 (MQ=255)
                                                                                                                                     cggGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGgc                  <  1:4017642/150‑1 (MQ=255)
                                                                                                                                     cggGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGgc                  <  2:2124553/150‑1 (MQ=255)
                                                                                                                                     cggGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGgc                  <  2:3234307/150‑1 (MQ=255)
                                                                                                                                       ggCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGgcgc                >  1:4215680/1‑150 (MQ=255)
                                                                                                                                         ccGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTg              <  1:4635592/150‑1 (MQ=255)
                                                                                                                                          cGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTgg             <  1:2991614/150‑1 (MQ=255)
                                                                                                                                          cGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTgg             >  1:1495126/1‑150 (MQ=255)
                                                                                                                                          cGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTgg             <  1:1788839/150‑1 (MQ=255)
                                                                                                                                          cGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTgg             >  1:3886326/1‑150 (MQ=255)
                                                                                                                                                gcaccggcaTGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTggcgggcg       <  2:2416716/150‑1 (MQ=255)
                                                                                                                                                  accggcaTGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGcgggcggg     <  1:3079690/150‑1 (MQ=255)
                                                                                                                                                    cggcaTGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGcc   >  2:2610153/1‑150 (MQ=255)
                                                                                                                                                    cggcaTGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGcc   >  2:3608689/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTCGGTGGTGGGCGAAGATGGCGTATACCACCACACCGGCCGGGTGCGGGTGTTTGTCTCGGAAGCGCAGGCGATCAAGGCGATCAAGCGGGAAGAGATTGTGCAGGGCGATATCATGGTGGTGATCGGCGGCGGGCCGTCCGGCACCGGCATGGAAGAGACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGCCG  >  CP009273/280284‑280582

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: