Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,990,354 C→A 10.6% H236Q (CAC→CAA)  xerC → site‑specific tyrosine recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,990,3540CA10.6% 1028.3 / 63.3 310H236Q (CAC→CAA) xerCsite‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base C (151/126);  new base A (17/16);  total (168/142)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 317 total aligned reads displayed.

GCGTGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGT  >  CP009273/3990205‑3990503
                                                                                                                                                     |                                                                                                                                                     
gCGTGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAc                                                                                                                                                       <  2:1472071/150‑1 (MQ=255)
gCGTGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAc                                                                                                                                                       <  2:3252696/150‑1 (MQ=255)
 cGTGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACg                                                                                                                                                      >  1:4109556/1‑150 (MQ=255)
 cGTGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACg                                                                                                                                                      >  1:3690393/1‑150 (MQ=255)
 cGTGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACg                                                                                                                                                      >  1:1791389/1‑150 (MQ=255)
 cGTGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACg                                                                                                                                                      >  1:4381929/1‑150 (MQ=255)
   tGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGtt                                                                                                                                                    <  2:3970557/150‑1 (MQ=255)
    ggATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTc                                                                                                                                                   <  2:212166/150‑1 (MQ=255)
    ggATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTc                                                                                                                                                   <  2:725247/150‑1 (MQ=255)
     gATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCa                                                                                                                                                  >  1:3958834/1‑150 (MQ=255)
     gATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCa                                                                                                                                                  >  1:2095700/1‑150 (MQ=255)
     gATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCa                                                                                                                                                  <  2:2652960/150‑1 (MQ=255)
         gAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCg                                                                                                                                              <  2:3999172/150‑1 (MQ=255)
          aGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGc                                                                                                                                             >  1:4316624/1‑150 (MQ=255)
          aGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGc                                                                                                                                             >  1:3008756/1‑150 (MQ=255)
              cTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATaa                                                                                                                                         >  1:3543724/1‑150 (MQ=255)
                 gCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAAtt                                                                                                                                      <  1:471979/150‑1 (MQ=255)
                 gCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAAtt                                                                                                                                      >  2:872078/1‑150 (MQ=255)
                 gCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAAtt                                                                                                                                      <  1:3707040/150‑1 (MQ=255)
                   ttGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTAc                                                                                                                                    <  1:3876532/150‑1 (MQ=255)
                     gATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGt                                                                                                                                  <  1:4244515/150‑1 (MQ=255)
                      aTTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTc                                                                                                                                 <  1:2687430/150‑1 (MQ=255)
                       tttGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCa                                                                                                                                <  2:1215322/150‑1 (MQ=255)
                       tttGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCa                                                                                                                                <  1:4550209/150‑1 (MQ=255)
                         tGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACt                                                                                                                              <  1:4881859/150‑1 (MQ=255)
                         tGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCACt                                                                                                                              <  2:5147252/150‑1 (MQ=255)
                           cgcgACCTGTTGGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCCTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCg                                                                                                                            <  1:3760768/150‑1 (MQ=255)
                             cgACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGtt                                                                                                                          <  1:2707918/150‑1 (MQ=255)
                              gACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTc                                                                                                                         >  2:1989135/1‑150 (MQ=255)
                              gACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTc                                                                                                                         <  2:4620706/150‑1 (MQ=255)
                              gACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTc                                                                                                                         <  2:4867529/150‑1 (MQ=255)
                                   gTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCAc                                                                                                                    >  1:3140921/1‑150 (MQ=255)
                                        gTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCata                                                                                                               >  2:3886620/1‑150 (MQ=255)
                                          aGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATg                                                                                                             <  1:2024837/150‑1 (MQ=255)
                                             gAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTg                                                                                                          <  1:369951/150‑1 (MQ=255)
                                             gAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTg                                                                                                          >  1:4419774/1‑150 (MQ=255)
                                               aGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGa                                                                                                        <  1:1660893/150‑1 (MQ=255)
                                               aGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGa                                                                                                        <  1:2393084/150‑1 (MQ=255)
                                                gacgacGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAg                                                                                                       <  1:1936313/150‑1 (MQ=255)
                                                gacgacGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGCGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGCCACTCGTTCGCCACGCATATGCTGGAg                                                                                                       <  1:3065550/150‑1 (MQ=255)
                                                  cgacgCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTc                                                                                                     >  2:4720550/1‑150 (MQ=255)
                                                  cgacgCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTc                                                                                                     <  2:3712605/150‑1 (MQ=255)
                                                  cgacgCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTc                                                                                                     >  1:4989274/1‑150 (MQ=255)
                                                   gacgCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCg                                                                                                    >  2:654384/1‑150 (MQ=255)
                                                    acgCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGa                                                                                                   <  1:3462837/150‑1 (MQ=255)
                                                    acgCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGa                                                                                                   <  2:2252920/150‑1 (MQ=255)
                                                       cgcTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAgcg                                                                                                >  2:2832356/1‑150 (MQ=255)
                                                              tCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCt                                                                                         >  2:202989/1‑150 (MQ=255)
                                                                tGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTc                                                                                       >  2:5053460/1‑150 (MQ=255)
                                                                  tCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCgt                                                                                     >  1:4487872/1‑150 (MQ=255)
                                                                       aCTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGgtgt                                                                                >  2:2737878/1‑150 (MQ=255)
                                                                          gggCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCa                                                                             >  1:628898/1‑150 (MQ=255)
                                                                          gggCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCa                                                                             >  2:1179183/1‑150 (MQ=255)
                                                                           ggCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAg                                                                            >  2:1673543/1‑150 (MQ=255)
                                                                            gCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAgg                                                                           >  1:1500746/1‑150 (MQ=255)
                                                                             cAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGa                                                                          >  1:1426453/1‑150 (MQ=255)
                                                                              aaGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAg                                                                         >  2:3350557/1‑150 (MQ=255)
                                                                                gCGTATCTCCGCGCGTACTGTGCAGAAACGCTTTGCCGAATGGGGCATAAACCCAGGGCTGAAGAATCACGTTCATCCGCAGAAAGTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAgct                                                                       <  2:240675/150‑1 (MQ=255)
                                                                                    aTCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGctgctg                                                                   <  1:4937087/150‑1 (MQ=255)
                                                                                      ctcCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTggg                                                                 <  2:5020146/150‑1 (MQ=255)
                                                                                      ctcCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTggg                                                                 >  1:1063460/1‑150 (MQ=255)
                                                                                        ccGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTc                                                               <  2:4925837/150‑1 (MQ=255)
                                                                                               tAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCaa                                                        >  2:1803936/1‑150 (MQ=255)
                                                                                                  tgtgCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACct                                                     >  2:100620/1‑150 (MQ=255)
                                                                                                    tgCAGAAACGCTTTGCCGAATGGGGCATAAAACAATGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTTGTGTGCAGGAGCTGCTGGGTCATGCCAACctct                                                   <  1:717304/150‑1 (MQ=255)
                                                                                                    tgCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACctct                                                   <  1:3264963/150‑1 (MQ=255)
                                                                                                     gCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCtctc                                                  >  1:1122120/1‑150 (MQ=255)
                                                                                                       aGAAACGCTTTGCCGAATGGGGCATAAAACAATGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCTCCACGCATATGCTGGAGTCGAGCGGAGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTcca                                                >  1:2783801/1‑150 (MQ=255)
                                                                                                        gAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTccac                                               <  2:3435493/150‑1 (MQ=255)
                                                                                                         aaaCGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTccacc                                              >  1:4738107/1‑150 (MQ=255)
                                                                                                             gCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGc                                          >  1:1638292/1‑150 (MQ=255)
                                                                                                              cTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCa                                         <  2:2358490/150‑1 (MQ=255)
                                                                                                              cTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCa                                         >  1:87858/1‑150 (MQ=255)
                                                                                                              cTTGGCCGAATGGGGCATCACACAAGGGCTGAATAATCACGTTCAGCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCa                                         <  2:2957791/150‑1 (MQ=255)
                                                                                                                ttGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCaaa                                       <  1:2893439/150‑1 (MQ=255)
                                                                                                                ttGCCGAATGGGGCATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCaaa                                       >  1:2230234/1‑150 (MQ=255)
                                                                                                                  gCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATc                                     >  2:3778967/1‑150 (MQ=255)
                                                                                                                  gCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATc                                     <  2:1516049/150‑1 (MQ=255)
                                                                                                                   ccGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCt                                    <  2:795143/150‑1 (MQ=255)
                                                                                                                        tGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACt                               >  2:4338407/1‑150 (MQ=255)
                                                                                                                        tGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACt                               <  1:531810/150‑1 (MQ=255)
                                                                                                                            gCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATc                           >  1:2963170/1‑150 (MQ=255)
                                                                                                                             cATAAAACAAGGGCTGAATAATCACGTTCATCCGCCTAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGAGCTTCGTGGTGTGCAGGAGCTGCTGGGTGATGCCAACCTCTCCACCCCGCTCATCTAGACTCagcg                          >  2:5076780/1‑147 (MQ=255)
                                                                                                                             cATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCt                          >  1:2162953/1‑150 (MQ=255)
                                                                                                                             cATAAAACAAGGGCTGAATAATCAAGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCt                          <  2:4134190/150‑1 (MQ=255)
                                                                                                                                   aCAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGAttt                    >  1:4160616/1‑150 (MQ=255)
                                                                                                                                     aaGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTc                  >  2:645933/1‑150 (MQ=255)
                                                                                                                                     aaGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTc                  <  2:1080781/150‑1 (MQ=255)
                                                                                                                                      aGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCa                 <  1:1359861/150‑1 (MQ=255)
                                                                                                                                       gggCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCaa                >  1:4657030/1‑150 (MQ=255)
                                                                                                                                           tGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACAcc            >  2:3289421/1‑150 (MQ=255)
                                                                                                                                            gAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCt           >  1:2551113/1‑150 (MQ=255)
                                                                                                                                            gAATAATCAAGTTCATCCGCATAAATTACGTCTCTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTGCTCTTGATTTTCTACACCt           >  2:4835655/1‑150 (MQ=255)
                                                                                                                                             aataatCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCtt          <  2:4165338/150‑1 (MQ=255)
                                                                                                                                                aatCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGcc       >  1:5036298/1‑150 (MQ=255)
                                                                                                                                                 atCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCt      <  2:589167/150‑1 (MQ=255)
                                                                                                                                                 atCAAGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCGGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCt      <  2:2709047/150‑1 (MQ=255)
                                                                                                                                                  tCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTc     <  2:3690393/150‑1 (MQ=255)
                                                                                                                                                   cACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCg    >  1:4007559/1‑150 (MQ=255)
                                                                                                                                                    aCGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCgg   >  1:1036914/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTGGATTGAGCACTGGCTTGATTTGCGCGACCTGTTTGGTAGCGAAGACGACGCGCTTTTTCTGTCGAAACTGGGCAAGCGTATCTCCGCGCGTAATGTGCAGAAACGCTTTGCCGAATGGGGCATAAAACAAGGGCTGAATAATCACGTTCATCCGCATAAATTACGTCACTCGTTCGCCACGCATATGCTGGAGTCGAGCGGCGATCTTCGTGGTGTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGT  >  CP009273/3990205‑3990503

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: ‑