Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,873,087 G→T 100% 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,0870GT95.2% 796.7 / 27.2 252V15V (GTC→GTAygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base G (9/3);  new base T (119/121);  total (128/124)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 252 total aligned reads displayed.

GCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCG  >  CP009273/2872941‑2873232
                                                                                                                                                  |                                                                                                                                                 
gCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGATAGATCATGGATAcg                                                                                                                                                <  1:2231685/150‑1 (MQ=255)
  tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                              >  1:3706209/1‑150 (MQ=255)
  tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                              <  2:643677/150‑1 (MQ=255)
   gCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGa                                                                                                                                             <  2:2628070/150‑1 (MQ=255)
     aTGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATc                                                                                                                                           <  2:2401953/150‑1 (MQ=255)
     aTGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATc                                                                                                                                           <  2:2402231/150‑1 (MQ=255)
      tGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCt                                                                                                                                          <  2:2050445/150‑1 (MQ=255)
         gAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTg                                                                                                                                       <  2:3232596/150‑1 (MQ=255)
          aaaCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGa                                                                                                                                      <  2:3706200/150‑1 (MQ=255)
          aaaCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGa                                                                                                                                      >  2:4514288/1‑150 (MQ=255)
          aaaCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGa                                                                                                                                      >  1:2591030/1‑150 (MQ=255)
              ccGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTgg                                                                                                                                  <  1:665319/150‑1 (MQ=255)
                 tgtACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAAt                                                                                                                               >  2:1948331/1‑150 (MQ=255)
                   tACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATgg                                                                                                                             >  1:1342293/1‑150 (MQ=255)
                    aCCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATggg                                                                                                                            <  1:3355877/150‑1 (MQ=255)
                      cAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGAt                                                                                                                          >  2:3636651/1‑150 (MQ=255)
                      cAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGAt                                                                                                                          <  2:22071/150‑1 (MQ=255)
                      cAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGAt                                                                                                                          <  2:23023/150‑1 (MQ=255)
                       aGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGAtt                                                                                                                         <  2:399236/150‑1 (MQ=255)
                        ggTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTa                                                                                                                        >  1:4341995/1‑150 (MQ=255)
                        ggTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTa                                                                                                                        >  1:2761126/1‑150 (MQ=255)
                           tCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTAGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgg                                                                                                                     >  2:3511115/1‑150 (MQ=255)
                           tCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTAGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgg                                                                                                                     >  2:3507325/1‑150 (MQ=255)
                                cATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCACTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACCTCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGACTGGGAGTAAGGGGAAg                                                                                                                >  2:2501202/1‑150 (MQ=255)
                                 atgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGc                                                                                                               >  2:1143709/1‑150 (MQ=255)
                                  tgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGcc                                                                                                              >  2:2176376/1‑150 (MQ=255)
                                      gCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGt                                                                                                          <  1:78961/150‑1 (MQ=255)
                                      gCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGt                                                                                                          <  1:460941/150‑1 (MQ=255)
                                       cAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTc                                                                                                         <  1:3551997/150‑1 (MQ=255)
                                           cAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAttt                                                                                                     >  2:3326321/1‑150 (MQ=255)
                                           cAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGAGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGAGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAttt                                                                                                     >  2:3326323/1‑150 (MQ=255)
                                            aaCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAtttt                                                                                                    >  2:3331363/1‑150 (MQ=255)
                                             aCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTtta                                                                                                   <  1:1415757/150‑1 (MQ=255)
                                               cGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTttatt                                                                                                 <  2:847930/150‑1 (MQ=255)
                                                 aGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTac                                                                                               <  2:1636562/150‑1 (MQ=255)
                                                   ccAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTacac                                                                                             <  1:1959923/150‑1 (MQ=255)
                                                    cAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACAcc                                                                                            >  1:1920600/1‑150 (MQ=255)
                                                      aCAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTc                                                                                          >  1:3214652/1‑150 (MQ=255)
                                                        aGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGGTTAAGGGGAAGCCAGGTCATTTTATTACACCTCaa                                                                                        >  1:3942741/1‑150 (MQ=255)
                                                        aGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCaa                                                                                        >  1:714501/1‑150 (MQ=255)
                                                         ggAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAAt                                                                                       >  1:1729021/1‑150 (MQ=255)
                                                         ggAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAAt                                                                                       >  2:423907/1‑150 (MQ=255)
                                                          gAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATc                                                                                      >  1:4366728/1‑150 (MQ=255)
                                                          gAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATc                                                                                      >  1:2089667/1‑150 (MQ=255)
                                                                   tGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAg                                                                             >  1:2973610/1‑150 (MQ=255)
                                                                     cGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGcc                                                                           <  1:1045110/150‑1 (MQ=255)
                                                                     cGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGcc                                                                           >  1:546912/1‑150 (MQ=255)
                                                                       gATCCCTGTCGTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGTGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCaa                                                                         <  2:2547706/150‑1 (MQ=255)
                                                                         tCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAg                                                                       >  2:2527957/1‑150 (MQ=255)
                                                                           ccTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGAt                                                                     <  2:693851/150‑1 (MQ=255)
                                                                            cTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATTGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATa                                                                    <  2:2202890/150‑1 (MQ=255)
                                                                            cTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATa                                                                    >  1:2058103/1‑150 (MQ=255)
                                                                             tGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAg                                                                   >  2:323872/1‑150 (MQ=255)
                                                                              gTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGc                                                                  <  2:4286213/150‑1 (MQ=255)
                                                                               tCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCa                                                                 <  1:4512952/150‑1 (MQ=255)
                                                                               tCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCa                                                                 >  2:512171/1‑150 (MQ=255)
                                                                                  tGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGc                                                              >  1:1670101/1‑150 (MQ=255)
                                                                                   gTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGcc                                                             >  1:3787164/1‑150 (MQ=255)
                                                                                    tCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCca                                                            <  1:3684648/150‑1 (MQ=255)
                                                                                       aTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAt                                                         >  1:3242567/1‑150 (MQ=255)
                                                                                       aTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAt                                                         >  1:3815580/1‑150 (MQ=255)
                                                                                           gTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCa                                                     >  1:1998616/1‑150 (MQ=255)
                                                                                           gTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCa                                                     >  1:1998633/1‑150 (MQ=255)
                                                                                            tACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCAt                                                    >  1:4134988/1‑150 (MQ=255)
                                                                                             aCAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATc                                                   <  1:3131048/150‑1 (MQ=255)
                                                                                                 aCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAtt                                               >  2:3395844/1‑150 (MQ=255)
                                                                                                 aCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAtt                                               <  2:666558/150‑1 (MQ=255)
                                                                                                 aCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAtt                                               <  2:2945604/150‑1 (MQ=255)
                                                                                                     gcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGc                                           <  1:2856580/150‑1 (MQ=255)
                                                                                                     gcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGc                                           <  1:2996543/150‑1 (MQ=255)
                                                                                                      ccATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCt                                          <  1:2838951/150‑1 (MQ=255)
                                                                                                      ccATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCt                                          <  1:1994922/150‑1 (MQ=255)
                                                                                                          cTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGcc                                      <  2:1662321/150‑1 (MQ=255)
                                                                                                              tACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAAt                                  >  2:4001322/1‑150 (MQ=255)
                                                                                                              tACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAAt                                  >  1:4445254/1‑150 (MQ=255)
                                                                                                              tACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAAt                                  >  2:4001346/1‑150 (MQ=255)
                                                                                                              tACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAAt                                  <  1:1168974/150‑1 (MQ=255)
                                                                                                                  tCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCt                              >  1:392848/1‑150 (MQ=255)
                                                                                                                     aTCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGc                           <  2:2761126/150‑1 (MQ=255)
                                                                                                                     aTCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGc                           >  2:1960718/1‑150 (MQ=255)
                                                                                                                     aTCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGc                           <  2:4341986/150‑1 (MQ=255)
                                                                                                                      tCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCt                          >  1:3387947/1‑150 (MQ=255)
                                                                                                                      tCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCt                          >  1:3388116/1‑150 (MQ=255)
                                                                                                                        tGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGt                        >  1:2709045/1‑150 (MQ=255)
                                                                                                                         gCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTa                       <  2:878421/150‑1 (MQ=255)
                                                                                                                         gCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTa                       <  2:2924945/150‑1 (MQ=255)
                                                                                                                         gCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTa                       <  2:2088061/150‑1 (MQ=255)
                                                                                                                             atatTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAg                   >  1:4265061/1‑150 (MQ=255)
                                                                                                                              tatTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGa                  <  1:1325263/150‑1 (MQ=255)
                                                                                                                               atTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGAt                 <  2:1269354/150‑1 (MQ=255)
                                                                                                                                 tGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCt               <  2:2988734/150‑1 (MQ=255)
                                                                                                                                    aGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAtt            >  1:1636924/1‑150 (MQ=255)
                                                                                                                                     gAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTa           >  2:4342148/1‑150 (MQ=255)
                                                                                                                                      aaaGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTaa          >  1:3470858/1‑150 (MQ=255)
                                                                                                                                          aTCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCt      <  2:3757844/150‑1 (MQ=255)
                                                                                                                                          aTCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCt      <  2:2225551/150‑1 (MQ=255)
                                                                                                                                          aTCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCt      >  2:1210812/1‑150 (MQ=255)
                                                                                                                                             aTGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc   <  1:995190/150‑1 (MQ=255)
                                                                                                                                             aTGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc   <  1:3273720/150‑1 (MQ=255)
                                                                                                                                             aTGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc   <  1:1840834/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCG  >  CP009273/2872941‑2873232

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: