Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,873,087 G→T 60.9% V15V (GTC→GTA ygbT ← multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,873,0870GT60.9% 207.4 / inf 310V15V (GTC→GTAygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base G (65/56);  new base T (88/101);  total (153/157)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 316 total aligned reads displayed.

CAGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGT  >  CP009273/2872939‑2873235
                                                                                                                                                    |                                                                                                                                                    
cAGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGgaga                                                                                                                                                     <  2:4753114/150‑1 (MQ=255)
  gCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTATAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGAcg                                                                                                                                                   >  2:2645112/1‑150 (MQ=255)
  gCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcg                                                                                                                                                   <  1:2962055/150‑1 (MQ=255)
    tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                                 >  1:2625956/1‑150 (MQ=255)
         gCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCtt                                                                                                                                            <  1:4800535/150‑1 (MQ=255)
          cGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCttt                                                                                                                                           <  1:2550089/150‑1 (MQ=255)
            aaaCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGa                                                                                                                                         >  1:5177248/1‑150 (MQ=255)
              aCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGt                                                                                                                                       <  2:697397/150‑1 (MQ=255)
              aCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGt                                                                                                                                       <  1:2340089/150‑1 (MQ=255)
              aCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGt                                                                                                                                       >  1:5107102/1‑150 (MQ=255)
               cccGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTg                                                                                                                                      <  1:801163/150‑1 (MQ=255)
                 cGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGa                                                                                                                                    <  2:1067784/150‑1 (MQ=255)
                 cGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGa                                                                                                                                    <  2:3694054/150‑1 (MQ=255)
                   tgtACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCCTATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAAt                                                                                                                                  >  2:3631069/1‑150 (MQ=255)
                     tACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATgg                                                                                                                                >  2:3555341/1‑150 (MQ=255)
                      aCCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATggg                                                                                                                               <  2:1057122/150‑1 (MQ=255)
                         aGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGAtt                                                                                                                            <  1:4406690/150‑1 (MQ=255)
                          ggTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTa                                                                                                                           >  2:4548868/1‑150 (MQ=255)
                           gTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTaa                                                                                                                          >  2:3833474/1‑150 (MQ=255)
                              ccAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAggg                                                                                                                       <  2:5124999/150‑1 (MQ=255)
                               cAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgggg                                                                                                                      <  2:444566/150‑1 (MQ=255)
                               cAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgggg                                                                                                                      >  1:3946662/1‑150 (MQ=255)
                                aGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGa                                                                                                                     >  1:1961552/1‑150 (MQ=255)
                                aGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGa                                                                                                                     >  1:1598707/1‑150 (MQ=255)
                                aGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGa                                                                                                                     >  1:4322525/1‑150 (MQ=255)
                                 gCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGaa                                                                                                                    >  1:4497314/1‑150 (MQ=255)
                                  cATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAg                                                                                                                   <  1:1274750/150‑1 (MQ=255)
                                    tgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACACACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGcc                                                                                                                 <  1:2229429/150‑1 (MQ=255)
                                     gatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCa                                                                                                                <  2:2197445/150‑1 (MQ=255)
                                         cAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTc                                                                                                            <  1:1164504/150‑1 (MQ=255)
                                          aGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCa                                                                                                           <  1:4104855/150‑1 (MQ=255)
                                          aGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCa                                                                                                           >  1:2598452/1‑150 (MQ=255)
                                          aGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCGTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCa                                                                                                           >  2:2712178/1‑150 (MQ=255)
                                          aGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCa                                                                                                           >  1:4459396/1‑150 (MQ=255)
                                          aGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCa                                                                                                           >  1:4549132/1‑150 (MQ=255)
                                           ggCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAt                                                                                                          <  2:2881063/150‑1 (MQ=255)
                                             cAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAttt                                                                                                        <  2:1898515/150‑1 (MQ=255)
                                             cAACCGAGCCAACAGGAATATGAGTGCGGATCCCGGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCATTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAttt                                                                                                        <  2:4307995/150‑1 (MQ=255)
                                              aaCCGAGCCAACATGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAtttt                                                                                                       <  2:4947713/150‑1 (MQ=255)
                                               aCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTtta                                                                                                      <  2:3575268/150‑1 (MQ=255)
                                                ccGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTttat                                                                                                     >  1:3885714/1‑150 (MQ=255)
                                                   aGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTac                                                                                                  <  1:544017/150‑1 (MQ=255)
                                                       aaCAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCt                                                                                              >  1:4461439/1‑150 (MQ=255)
                                                        aCAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTc                                                                                             <  1:2610995/150‑1 (MQ=255)
                                                          aGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCaa                                                                                           >  2:603900/1‑150 (MQ=255)
                                                              atatGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATcac                                                                                       >  1:2518798/1‑150 (MQ=255)
                                                                atGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAg                                                                                     >  2:1099796/1‑150 (MQ=255)
                                                                 tGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGt                                                                                    >  2:1814966/1‑150 (MQ=255)
                                                                 tGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGACATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGt                                                                                    <  2:1996444/150‑1 (MQ=255)
                                                                     tGCGGATCCCTGTCTTGTCGTTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAg                                                                                <  1:4944702/150‑1 (MQ=255)
                                                                       cGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGcc                                                                              >  1:4660775/1‑150 (MQ=255)
                                                                                  cTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGTATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCaa                                                                   <  1:5230013/150‑1 (MQ=255)
                                                                                   ttGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAg                                                                  <  2:3581819/150‑1 (MQ=255)
                                                                                      tCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCca                                                               <  2:2975893/150‑1 (MQ=255)
                                                                                         aTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAt                                                            <  2:2753152/150‑1 (MQ=255)
                                                                                         aTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAt                                                            <  2:2296607/150‑1 (MQ=255)
                                                                                           aaGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATcc                                                          >  1:4485071/1‑150 (MQ=255)
                                                                                           aaGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAGcc                                                          >  1:2700672/1‑150 (MQ=255)
                                                                                            aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                         <  1:3459978/150‑1 (MQ=255)
                                                                                            aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                         >  2:1528442/1‑150 (MQ=255)
                                                                                              tACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCAt                                                       >  1:4491620/1‑150 (MQ=255)
                                                                                                cAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCg                                                     <  2:3118502/150‑1 (MQ=255)
                                                                                                   aCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAtt                                                  <  1:3051642/150‑1 (MQ=255)
                                                                                                    cgcgCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAttt                                                 <  2:1952853/150‑1 (MQ=255)
                                                                                                         cATCTATTACATCGATCTGCCCATATTGCAGAAAGATCCTGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTg                                            <  2:4097644/150‑1 (MQ=255)
                                                                                                         cATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTg                                            >  2:1632227/1‑150 (MQ=255)
                                                                                                         cATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTg                                            <  2:1774574/150‑1 (MQ=255)
                                                                                                          aTCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTgg                                           >  1:2322239/1‑150 (MQ=255)
                                                                                                          aTCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTgg                                           <  1:3555341/150‑1 (MQ=255)
                                                                                                              aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                       <  1:1253931/150‑1 (MQ=255)
                                                                                                              aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                       <  1:1242419/150‑1 (MQ=255)
                                                                                                              aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                       <  2:2458782/150‑1 (MQ=255)
                                                                                                              aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                       <  2:4621123/150‑1 (MQ=255)
                                                                                                              aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                       <  2:4814645/150‑1 (MQ=255)
                                                                                                              aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                       <  1:3548478/150‑1 (MQ=255)
                                                                                                                 aCATCGATCTGCCACTATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATa                                    <  2:3994299/150‑1 (MQ=255)
                                                                                                                  cATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATAc                                   >  2:1625198/1‑150 (MQ=255)
                                                                                                                  cATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATAc                                   >  2:763657/1‑150 (MQ=255)
                                                                                                                   aTCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATAcc                                  >  1:1084642/1‑150 (MQ=255)
                                                                                                                    tCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATACCATCGATTTAGCTGGCCCAATACCt                                 >  2:3665922/1‑150 (MQ=255)
                                                                                                                     cGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATTTGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATAACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCtt                                >  1:1545980/1‑150 (MQ=255)
                                                                                                                      gATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTg                               >  1:4061916/1‑150 (MQ=255)
                                                                                                                      gATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTg                               >  1:284978/1‑150 (MQ=255)
                                                                                                                      gATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTg                               <  2:4129740/150‑1 (MQ=255)
                                                                                                                      gATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTg                               >  1:3362533/1‑150 (MQ=255)
                                                                                                                           gCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCCCAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTATCTGGCCCAATACCTTGCTGTa                          <  2:5106806/150‑1 (MQ=255)
                                                                                                                           gCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTa                          >  1:2891827/1‑150 (MQ=255)
                                                                                                                              cATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACaa                       >  1:1688753/1‑150 (MQ=255)
                                                                                                                              cATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACaa                       >  1:979885/1‑150 (MQ=255)
                                                                                                                                tatTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGa                     >  2:144241/1‑150 (MQ=255)
                                                                                                                                   tGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCt                  <  1:1892469/150‑1 (MQ=255)
                                                                                                                                     cAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAt                >  2:4193156/1‑150 (MQ=255)
                                                                                                                                     cAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAt                >  2:2346673/1‑150 (MQ=255)
                                                                                                                                      aGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAtt               <  2:2330702/150‑1 (MQ=255)
                                                                                                                                      aGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAtt               <  1:3824432/150‑1 (MQ=255)
                                                                                                                                          aGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACg           <  1:3697153/150‑1 (MQ=255)
                                                                                                                                           gATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACATTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGc          >  2:3483287/1‑150 (MQ=255)
                                                                                                                                           gATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACATTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGc          >  2:2452041/1‑150 (MQ=255)
                                                                                                                                                  gagaCGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGcgcg   >  2:4754587/1‑150 (MQ=255)
                                                                                                                                                   agaCGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGcgt  >  2:5196410/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGT  >  CP009273/2872939‑2873235

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: