Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,666,401 C→A 5.5% G332G (GGC→GGA hcaD → phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,666,4010CA 5.5% 991.8 / 21.9 262G332G (GGC→GGAhcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base C (143/104);  new base A (7/8);  total (150/112)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

ATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATCCAGAGCGG  >  CP009273/2666253‑2666550
                                                                                                                                                    |                                                                                                                                                     
aTGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCg                                                                                                                                                      >  1:1850330/1‑150 (MQ=255)
   gTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGata                                                                                                                                                   <  2:4211058/150‑1 (MQ=255)
     gCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATg                                                                                                                                                 <  2:3265960/150‑1 (MQ=255)
      cACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGc                                                                                                                                                >  2:3376097/1‑150 (MQ=255)
           cacCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGc                                                                                                                                           <  2:4040598/150‑1 (MQ=255)
             ccGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGa                                                                                                                                         >  2:2437693/1‑150 (MQ=255)
               gCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGAGATATGCGTGGCGATg                                                                                                                                       <  1:1384196/150‑1 (MQ=255)
                 tGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTCCCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGAc                                                                                                                                     >  1:3662718/1‑150 (MQ=255)
                 tGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGAc                                                                                                                                     >  1:2405174/1‑150 (MQ=255)
                    gAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTgg                                                                                                                                  >  2:1293611/1‑150 (MQ=255)
                       aGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCtt                                                                                                                               <  1:2289940/150‑1 (MQ=255)
                        gCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCttt                                                                                                                              <  2:3324849/150‑1 (MQ=255)
                        gCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCttt                                                                                                                              >  1:2052036/1‑150 (MQ=255)
                        gCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCttt                                                                                                                              <  1:965170/150‑1 (MQ=255)
                          tGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGt                                                                                                                            <  2:3060389/150‑1 (MQ=255)
                          tGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGt                                                                                                                            >  1:267911/1‑150 (MQ=255)
                                aaCGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGc                                                                                                                      >  1:3059901/1‑150 (MQ=255)
                                aaCGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGc                                                                                                                      >  1:3059551/1‑150 (MQ=255)
                                aaCGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGc                                                                                                                      >  1:3059496/1‑150 (MQ=255)
                                  cGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCaa                                                                                                                    >  2:4471811/1‑150 (MQ=255)
                                     cAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATCACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAAccc                                                                                                                 <  1:1798772/150‑1 (MQ=255)
                                       aTAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCgg                                                                                                               >  2:1481569/1‑150 (MQ=255)
                                         aaCCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGaa                                                                                                             >  1:4380408/1‑150 (MQ=255)
                                         aaCCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGaa                                                                                                             >  2:1135196/1‑150 (MQ=255)
                                             aGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTc                                                                                                         >  1:2503969/1‑150 (MQ=255)
                                               gcgcAACTTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAg                                                                                                       <  2:4279978/150‑1 (MQ=255)
                                                cgcAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGa                                                                                                      <  2:181612/150‑1 (MQ=255)
                                                  cAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTTGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCAGGAAACTCAGAAg                                                                                                    <  2:2028385/150‑1 (MQ=255)
                                                  cAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATGGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAg                                                                                                    <  2:2653073/150‑1 (MQ=255)
                                                   aaaTTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAgg                                                                                                   <  2:368189/150‑1 (MQ=255)
                                                   aaaTTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAgg                                                                                                   <  2:1107963/150‑1 (MQ=255)
                                                    aaTTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGAGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGc                                                                                                  >  2:5501/1‑150 (MQ=255)
                                                     aTTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCg                                                                                                 >  1:1000769/1‑150 (MQ=255)
                                                         ccgctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGAttt                                                                                             <  1:1154701/150‑1 (MQ=255)
                                                          cgctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTg                                                                                            <  1:1970188/150‑1 (MQ=255)
                                                           gctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTgg                                                                                           >  2:1562174/1‑150 (MQ=255)
                                                            ctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                          >  1:2820710/1‑150 (MQ=255)
                                                            ctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                          >  1:1318473/1‑150 (MQ=255)
                                                            ctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                          >  1:535719/1‑150 (MQ=255)
                                                              gccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGttt                                                                                        <  2:2205599/150‑1 (MQ=255)
                                                                cgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCTGTACAGTGATAACTTTCAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTaa                                                                                      <  2:2627689/150‑1 (MQ=255)
                                                                cgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTaa                                                                                      <  1:3007783/150‑1 (MQ=255)
                                                                 gcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAAt                                                                                     >  1:1045921/1‑150 (MQ=255)
                                                                      gTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCa                                                                                <  2:203490/150‑1 (MQ=255)
                                                                      gTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCa                                                                                >  2:1925303/1‑150 (MQ=255)
                                                                        tAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCaaa                                                                              >  1:3158515/1‑150 (MQ=255)
                                                                         aGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCaaaa                                                                             >  1:712135/1‑150 (MQ=255)
                                                                           ggCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACg                                                                           <  1:1651612/150‑1 (MQ=255)
                                                                                ccgctaccgctaCTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCTACCCGGAAACTCAGAAGTCGTTTTGGTTTAATCTGCAAAACGGCGTg                                                                      <  2:1765269/150‑1 (MQ=255)
                                                                                    taccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTa                                                                  <  1:3640152/150‑1 (MQ=255)
                                                                                    taccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTa                                                                  >  2:3108166/1‑150 (MQ=255)
                                                                                    taccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTa                                                                  <  1:2323717/150‑1 (MQ=255)
                                                                                      ccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATc                                                                <  1:2881112/150‑1 (MQ=255)
                                                                                         ctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGt                                                             >  2:1406986/1‑150 (MQ=255)
                                                                                           acTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGc                                                           >  2:3437796/1‑150 (MQ=255)
                                                                                            cTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCg                                                          >  1:1229581/1‑150 (MQ=255)
                                                                                            cTGCCGCCGCCGTGGTTCTGGAGAGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCg                                                          >  1:4456305/1‑150 (MQ=255)
                                                                                              gccgccgccGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGt                                                        >  1:3574624/1‑150 (MQ=255)
                                                                                              gccgccgccGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGt                                                        >  1:3483706/1‑150 (MQ=255)
                                                                                                     ccgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTg                                                 >  1:92237/1‑150 (MQ=255)
                                                                                                     ccgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTg                                                 >  1:2030367/1‑150 (MQ=255)
                                                                                                      cgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCACAACGGCGTGCTTATCGGTGCGGTCACGCTGa                                                >  2:2318388/1‑150 (MQ=255)
                                                                                                      cgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGa                                                >  2:2316860/1‑150 (MQ=255)
                                                                                                       gTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGaa                                               >  1:2487263/1‑150 (MQ=255)
                                                                                                         ggTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATc                                             >  2:2709681/1‑150 (MQ=255)
                                                                                                         ggTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATc                                             >  2:2709591/1‑150 (MQ=255)
                                                                                                         ggTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATc                                             >  2:2284384/1‑150 (MQ=255)
                                                                                                          gTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCa                                            >  1:4058876/1‑150 (MQ=255)
                                                                                                          gTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCa                                            >  1:2452684/1‑150 (MQ=255)
                                                                                                             cTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAggg                                         <  1:126806/150‑1 (MQ=255)
                                                                                                             cTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAggg                                         <  2:279008/150‑1 (MQ=255)
                                                                                                              tGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAgggg                                        >  1:800039/1‑150 (MQ=255)
                                                                                                              tGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGAGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAgggg                                        >  1:3463472/1‑150 (MQ=255)
                                                                                                               ggAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGc                                       >  2:152531/1‑150 (MQ=255)
                                                                                                                gAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCg                                      >  1:3793464/1‑150 (MQ=255)
                                                                                                                  gCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTg                                    <  2:2204275/150‑1 (MQ=255)
                                                                                                                  gCGATCAGTACAGTGATAACTTACAGTTTATTGGAGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTg                                    <  1:2138664/150‑1 (MQ=255)
                                                                                                                   cGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTAGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTga                                   >  1:862939/1‑150 (MQ=255)
                                                                                                                       cAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGAtt                               >  1:2564064/1‑150 (MQ=255)
                                                                                                                        aGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTc                              <  1:176873/150‑1 (MQ=255)
                                                                                                                        aGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTc                              <  2:1850330/150‑1 (MQ=255)
                                                                                                                            cAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGccc                          <  1:2950882/150‑1 (MQ=255)
                                                                                                                             aGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCa                         >  2:4382011/1‑150 (MQ=255)
                                                                                                                             aGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCa                         >  1:4193186/1‑150 (MQ=255)
                                                                                                                             aGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCa                         >  1:585855/1‑150 (MQ=255)
                                                                                                                             aGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCa                         >  1:585718/1‑150 (MQ=255)
                                                                                                                                  tAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCg                    >  2:3779981/1‑150 (MQ=255)
                                                                                                                                   aaCTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGc                   >  1:2098333/1‑150 (MQ=255)
                                                                                                                                     cTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCaa                 >  2:3952963/1‑150 (MQ=255)
                                                                                                                                     cTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCaa                 >  1:2004585/1‑150 (MQ=255)
                                                                                                                                      ttACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCaaa                >  1:4049523/1‑150 (MQ=255)
                                                                                                                                           gTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCTGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGAt           >  2:4147325/1‑150 (MQ=255)
                                                                                                                                            tttATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATc          <  2:3167486/150‑1 (MQ=255)
                                                                                                                                              tATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCACATGGATCCa        >  1:3334047/1‑150 (MQ=255)
                                                                                                                                              tATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATCCa        >  1:2923165/1‑150 (MQ=255)
                                                                                                                                                ttGGAGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATCCaga      >  2:4288792/1‑150 (MQ=255)
                                                                                                                                                 tGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATCCagag     >  2:2077045/1‑150 (MQ=255)
                                                                                                                                                  ggCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATCCAGAGc    <  1:976993/150‑1 (MQ=255)
                                                                                                                                                   gCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATCCAGAGCg   <  2:4524304/150‑1 (MQ=255)
                                                                                                                                                    cGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATCCAGAGCgg  >  2:2951202/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATCCAGAGCGG  >  CP009273/2666253‑2666550

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: