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,0870GT97.4% 898.1 / 12.8 269V15V (GTC→GTAygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (143/118);  total (147/122)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 269 total aligned reads displayed.

AGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCG  >  CP009273/2872940‑2873232
                                                                                                                                                   |                                                                                                                                                 
aGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAc                                                                                                                                                 <  1:2317045/150‑1 (MQ=255)
 gCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcg                                                                                                                                                <  1:2080973/150‑1 (MQ=255)
   tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                              <  1:1021636/150‑1 (MQ=255)
   tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                              <  1:122102/150‑1 (MQ=255)
   tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                              <  1:3495675/150‑1 (MQ=255)
   tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                              <  1:3400456/150‑1 (MQ=255)
       tGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCt                                                                                                                                          <  2:2386697/150‑1 (MQ=255)
        gCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCtt                                                                                                                                         <  2:1540978/150‑1 (MQ=255)
          gAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTg                                                                                                                                       <  1:1527518/150‑1 (MQ=255)
           naaCCCGTGTACCAGGNTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGANCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGa                                                                                                                                      >  2:2167258/2‑150 (MQ=255)
            aaCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAg                                                                                                                                     <  1:2899662/150‑1 (MQ=255)
            aaCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAg                                                                                                                                     >  2:615522/1‑150 (MQ=255)
                cGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGa                                                                                                                                 <  2:2039328/150‑1 (MQ=255)
                cGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGa                                                                                                                                 <  2:4638189/150‑1 (MQ=255)
                  tgtACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAAt                                                                                                                               >  2:3790598/1‑150 (MQ=255)
                     aCCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATggg                                                                                                                            <  2:1827761/150‑1 (MQ=255)
                       cAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGAt                                                                                                                          >  2:1127197/1‑150 (MQ=255)
                         ggTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTa                                                                                                                        <  2:871541/150‑1 (MQ=255)
                         ggTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTa                                                                                                                        <  2:3516442/150‑1 (MQ=255)
                          gTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTaa                                                                                                                       >  1:1944561/1‑150 (MQ=255)
                           ttCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAg                                                                                                                      >  2:1438579/1‑150 (MQ=255)
                           ttCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAg                                                                                                                      >  2:2270251/1‑150 (MQ=255)
                             ccAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAggg                                                                                                                    <  1:3677016/150‑1 (MQ=255)
                              cAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgggg                                                                                                                   <  1:843221/150‑1 (MQ=255)
                              cAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgggg                                                                                                                   <  2:2819960/150‑1 (MQ=255)
                              cAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgggg                                                                                                                   >  1:3942181/1‑150 (MQ=255)
                               aGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGa                                                                                                                  <  1:1767460/150‑1 (MQ=255)
                                gCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGaa                                                                                                                 >  1:30039/1‑150 (MQ=255)
                                gCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGaa                                                                                                                 <  2:309605/150‑1 (MQ=255)
                                 cATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAg                                                                                                                >  1:3211864/1‑150 (MQ=255)
                                       gCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGt                                                                                                          <  1:949876/150‑1 (MQ=255)
                                        cAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTc                                                                                                         >  2:1506443/1‑150 (MQ=255)
                                        cAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTc                                                                                                         >  2:3350304/1‑150 (MQ=255)
                                            cAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAttt                                                                                                     <  2:1312805/150‑1 (MQ=255)
                                            cAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAttt                                                                                                     >  1:3961291/1‑150 (MQ=255)
                                              aCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTtta                                                                                                   <  1:4127681/150‑1 (MQ=255)
                                              aCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTtta                                                                                                   <  1:4203677/150‑1 (MQ=255)
                                               ccGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTttat                                                                                                  >  1:553529/1‑150 (MQ=255)
                                                 gAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTttatta                                                                                                >  2:2728619/1‑150 (MQ=255)
                                                  aGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTac                                                                                               <  2:2393343/150‑1 (MQ=255)
                                                     cAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACAcc                                                                                            <  1:4034976/150‑1 (MQ=255)
                                                     cAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACAcc                                                                                            >  2:1778005/1‑150 (MQ=255)
                                                      aaCAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCt                                                                                           >  2:2138095/1‑150 (MQ=255)
                                                       aCAGGAATATGAGTGCGGCTCCCTGTCTTGTCGATAAGTTCAAACGCGCAATCTATTACATCGTTCTTCCCTTTTTGCAGTAATATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTGATTTTATTACACCTc                                                                                          <  2:1947565/150‑1 (MQ=255)
                                                       aCAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTc                                                                                          <  1:2122143/150‑1 (MQ=255)
                                                           gAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATc                                                                                      >  1:1396705/1‑150 (MQ=255)
                                                            aaTATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATca                                                                                     >  1:4467872/1‑150 (MQ=255)
                                                            aaTATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATca                                                                                     >  2:589627/1‑150 (MQ=255)
                                                                tGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGt                                                                                 <  2:1993134/150‑1 (MQ=255)
                                                                tGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGt                                                                                 <  1:548477/150‑1 (MQ=255)
                                                                       ggATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCa                                                                          >  2:2947476/1‑150 (MQ=255)
                                                                       ggATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCa                                                                          >  2:4396590/1‑150 (MQ=255)
                                                                       ggATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCa                                                                          <  2:74599/150‑1 (MQ=255)
                                                                         aTCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCaaa                                                                        >  1:3133274/1‑150 (MQ=255)
                                                                                 cTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCaa                                                                >  1:34421/1‑150 (MQ=255)
                                                                                     tCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCca                                                            <  2:3189144/150‑1 (MQ=255)
                                                                                       gATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCcaca                                                          >  2:3935338/1‑150 (MQ=255)
                                                                                        aTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAt                                                         >  2:2154855/1‑150 (MQ=255)
                                                                                          aaGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATcc                                                       >  1:1673236/1‑150 (MQ=255)
                                                                                             tACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCAt                                                    >  1:2638320/1‑150 (MQ=255)
                                                                                             tACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCAt                                                    <  2:2552693/150‑1 (MQ=255)
                                                                                              aCAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATc                                                   <  1:613771/150‑1 (MQ=255)
                                                                                              aCAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATc                                                   <  1:3770978/150‑1 (MQ=255)
                                                                                               cAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAAACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTAATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCg                                                  >  2:3690811/1‑150 (MQ=255)
                                                                                                 aaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAt                                                >  1:1547271/1‑150 (MQ=255)
                                                                                                 aaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAt                                                <  2:4053403/150‑1 (MQ=255)
                                                                                                  aCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAtt                                               >  1:725901/1‑150 (MQ=255)
                                                                                                   cgcgCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAttt                                              <  2:4326752/150‑1 (MQ=255)
                                                                                                    gcgcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTa                                             <  2:2816817/150‑1 (MQ=255)
                                                                                                     cgcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAg                                            >  2:4210200/1‑150 (MQ=255)
                                                                                                      gcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGc                                           >  2:693240/1‑150 (MQ=255)
                                                                                                       ccATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCt                                          <  1:615484/150‑1 (MQ=255)
                                                                                                         aTCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTgg                                        >  2:1234432/1‑150 (MQ=255)
                                                                                                           cTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGcc                                      <  1:1452004/150‑1 (MQ=255)
                                                                                                             aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                    <  1:4052464/150‑1 (MQ=255)
                                                                                                               tACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAAt                                  >  1:3563087/1‑150 (MQ=255)
                                                                                                               tACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAAt                                  <  2:965230/150‑1 (MQ=255)
                                                                                                                 cATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATAc                                >  2:2118248/1‑150 (MQ=255)
                                                                                                                 cATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATAc                                >  2:3348187/1‑150 (MQ=255)
                                                                                                                     gATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTg                            <  1:318180/150‑1 (MQ=255)
                                                                                                                      aTCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGc                           >  1:2591613/1‑150 (MQ=255)
                                                                                                                      aTCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGc                           <  1:2213860/150‑1 (MQ=255)
                                                                                                                       tCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCt                          >  1:1765024/1‑150 (MQ=255)
                                                                                                                       tCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCt                          <  1:2423611/150‑1 (MQ=255)
                                                                                                                         tGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGt                        >  1:1562526/1‑150 (MQ=255)
                                                                                                                          gCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTa                       <  2:3897184/150‑1 (MQ=255)
                                                                                                                            ccATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACa                     >  1:3879102/1‑150 (MQ=255)
                                                                                                                             cATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACaa                    >  1:932051/1‑150 (MQ=255)
                                                                                                                               tatTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGa                  >  2:3111539/1‑150 (MQ=255)
                                                                                                                               tatTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGa                  <  1:4520040/150‑1 (MQ=255)
                                                                                                                                atTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGAt                 <  1:1110858/150‑1 (MQ=255)
                                                                                                                                   gCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTa              >  2:245644/1‑150 (MQ=255)
                                                                                                                                    cAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAt             >  1:2752118/1‑150 (MQ=255)
                                                                                                                                          gATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGc       >  2:3633018/1‑150 (MQ=255)
                                                                                                                                           aTCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCt      >  1:881415/1‑150 (MQ=255)
                                                                                                                                             cATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTgg    >  2:4384332/1‑150 (MQ=255)
                                                                                                                                             cATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTgg    >  1:2425865/1‑150 (MQ=255)
                                                                                                                                              aTGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc   >  1:826791/1‑150 (MQ=255)
                                                                                                                                              aTGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc   >  1:2609125/1‑150 (MQ=255)
                                                                                                                                               tGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgcg  >  1:350017/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCG  >  CP009273/2872940‑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: