Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,873,087 G→T 92.7% 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,0870GT92.7% 1011.2 / 66.3 330V15V (GTC→GTAygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base G (9/15);  new base T (157/148);  total (167/163)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 330 total aligned reads displayed.

AGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGC  >  CP009273/2872940‑2873233
                                                                                                                                                   |                                                                                                                                                  
aGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCACCCGAGCCAACAGGAATATGAGTGCGGATCCCGGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAc                                                                                                                                                  <  2:2151432/150‑1 (MQ=255)
aGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAc                                                                                                                                                  <  2:2897984/150‑1 (MQ=255)
 gCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcg                                                                                                                                                 <  2:1818695/150‑1 (MQ=255)
   tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                               >  1:786483/1‑150 (MQ=255)
   tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                               >  1:380603/1‑150 (MQ=255)
       tGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCt                                                                                                                                           <  2:703467/150‑1 (MQ=255)
         cGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCttt                                                                                                                                         <  2:537246/150‑1 (MQ=255)
         cGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCttt                                                                                                                                         >  1:1513753/1‑150 (MQ=255)
            aaCCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAg                                                                                                                                      <  2:3365904/150‑1 (MQ=255)
              cccGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTg                                                                                                                                    <  1:4746042/150‑1 (MQ=255)
               ccGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTgg                                                                                                                                   <  1:652805/150‑1 (MQ=255)
                  tgtACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAAt                                                                                                                                >  1:1059061/1‑150 (MQ=255)
                    tACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATgg                                                                                                                              <  1:3689382/150‑1 (MQ=255)
                      ccAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGa                                                                                                                            >  2:3879074/1‑150 (MQ=255)
                        aGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGAtt                                                                                                                          <  1:2798378/150‑1 (MQ=255)
                             ccAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAggg                                                                                                                     <  2:251102/150‑1 (MQ=255)
                             ccAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAggg                                                                                                                     <  2:106173/150‑1 (MQ=255)
                              cAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgggg                                                                                                                    >  2:3706126/1‑150 (MQ=255)
                                 cATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAg                                                                                                                 >  2:4880678/1‑150 (MQ=255)
                                  atgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGc                                                                                                                <  1:3013681/150‑1 (MQ=255)
                                   tgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGcc                                                                                                               <  1:1732894/150‑1 (MQ=255)
                                   tgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGcc                                                                                                               <  2:1423797/150‑1 (MQ=255)
                                    gatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGACTGGGATTAAGGGGACGCCa                                                                                                              >  2:4310467/1‑150 (MQ=255)
                                    gatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCa                                                                                                              >  2:191906/1‑150 (MQ=255)
                                    gatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCa                                                                                                              >  1:2290157/1‑150 (MQ=255)
                                     atgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAg                                                                                                             >  2:2214521/1‑150 (MQ=255)
                                     atgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAg                                                                                                             <  1:1190494/150‑1 (MQ=255)
                                     atgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTAGTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAg                                                                                                             >  2:3097405/1‑150 (MQ=255)
                                      tgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAgg                                                                                                            <  1:4889376/150‑1 (MQ=255)
                                       gCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGt                                                                                                           <  2:1750054/150‑1 (MQ=255)
                                        cAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTc                                                                                                          >  2:1271016/1‑150 (MQ=255)
                                         aGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCa                                                                                                         <  2:5110630/150‑1 (MQ=255)
                                         aGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCa                                                                                                         >  1:3063640/1‑150 (MQ=255)
                                         aGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCa                                                                                                         >  2:1986732/1‑150 (MQ=255)
                                          ggCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAt                                                                                                        <  2:2782567/150‑1 (MQ=255)
                                            cAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAttt                                                                                                      <  2:3073749/150‑1 (MQ=255)
                                             aaCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAtttt                                                                                                     >  2:3174574/1‑150 (MQ=255)
                                             aaCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAtttt                                                                                                     >  2:3590288/1‑150 (MQ=255)
                                             aaCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAtttt                                                                                                     <  1:2568602/150‑1 (MQ=255)
                                             aaCCGAGACAACAGGAATATGAGTGCGGATCCCTGTCTTTTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGTATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAtttt                                                                                                     >  2:1798279/1‑150 (MQ=255)
                                                     cAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACAcc                                                                                             <  1:3484861/150‑1 (MQ=255)
                                                      aaCAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCt                                                                                            >  2:3465433/1‑150 (MQ=255)
                                                          ggAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAAt                                                                                        >  1:60224/1‑150 (MQ=255)
                                                             atatGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATcac                                                                                     >  1:4245533/1‑150 (MQ=255)
                                                                  aGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTgg                                                                                <  1:2006295/150‑1 (MQ=255)
                                                                  aGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTgg                                                                                >  1:1962259/1‑150 (MQ=255)
                                                                  aGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTgg                                                                                <  1:1022650/150‑1 (MQ=255)
                                                                   gTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGa                                                                               >  2:1350980/1‑150 (MQ=255)
                                                                   gTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGa                                                                               <  2:1159445/150‑1 (MQ=255)
                                                                    tGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAg                                                                              <  2:2972413/150‑1 (MQ=255)
                                                                    tGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAg                                                                              >  1:3381560/1‑150 (MQ=255)
                                                                        gATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCaa                                                                          >  1:2538784/1‑150 (MQ=255)
                                                                        gATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCaa                                                                          >  1:517421/1‑150 (MQ=255)
                                                                         aTCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCaaa                                                                         >  2:2028749/1‑150 (MQ=255)
                                                                             cTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATa                                                                     >  2:4511537/1‑150 (MQ=255)
                                                                              tGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAg                                                                    <  1:1852449/150‑1 (MQ=255)
                                                                                  ttGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAg                                                                <  1:1341209/150‑1 (MQ=255)
                                                                                   tGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGc                                                               >  1:2237446/1‑150 (MQ=255)
                                                                                      cGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCcac                                                            >  2:2677855/1‑150 (MQ=255)
                                                                                        aTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAt                                                          >  2:2368772/1‑150 (MQ=255)
                                                                                         tAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATc                                                         >  2:1200716/1‑150 (MQ=255)
                                                                                         tAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATc                                                         <  2:119908/150‑1 (MQ=255)
                                                                                           aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                       >  1:1695120/1‑150 (MQ=255)
                                                                                           aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                       >  2:4489958/1‑150 (MQ=255)
                                                                                           aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                       <  2:4472467/150‑1 (MQ=255)
                                                                                            gTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCa                                                      <  2:3908071/150‑1 (MQ=255)
                                                                                              aCAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATc                                                    <  1:475927/150‑1 (MQ=255)
                                                                                                aaaCGCGCCATCTATTGCATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGa                                                  <  1:3156686/150‑1 (MQ=255)
                                                                                                aaaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGa                                                  >  1:1864987/1‑150 (MQ=255)
                                                                                                aaaCGCGCCATCTATTACATCGATCTGCCCATATTACAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGa                                                  >  1:4653125/1‑150 (MQ=255)
                                                                                                 aaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAt                                                 <  2:4553287/150‑1 (MQ=255)
                                                                                                    gcgcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTa                                              <  2:4652124/150‑1 (MQ=255)
                                                                                                    gcgcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTa                                              >  2:2389634/1‑150 (MQ=255)
                                                                                                       ccATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCt                                           <  1:3269617/150‑1 (MQ=255)
                                                                                                            tATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTCTTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGccc                                      >  1:3588109/1‑150 (MQ=255)
                                                                                                             aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                     <  1:2667779/150‑1 (MQ=255)
                                                                                                              ttACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCaa                                    <  1:1775637/150‑1 (MQ=255)
                                                                                                               tACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAAt                                   <  2:4076904/150‑1 (MQ=255)
                                                                                                                aCATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATa                                  >  2:790797/1‑150 (MQ=255)
                                                                                                                aCATCGATCTGCCCATATTGCAGAAAGATCATGGAAACGCGAACTTTTAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATa                                  >  1:5039893/1‑150 (MQ=255)
                                                                                                                    cGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCtt                              <  1:236202/150‑1 (MQ=255)
                                                                                                                    cGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCtt                              <  1:895094/150‑1 (MQ=255)
                                                                                                                     gATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTg                             >  2:1642791/1‑150 (MQ=255)
                                                                                                                      aTCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGc                            >  2:4002801/1‑150 (MQ=255)
                                                                                                                       tCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCt                           <  2:1659861/150‑1 (MQ=255)
                                                                                                                       tCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCt                           >  1:2844413/1‑150 (MQ=255)
                                                                                                                        cTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTg                          <  1:4065069/150‑1 (MQ=255)
                                                                                                                           cccATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTAc                       <  1:4594774/150‑1 (MQ=255)
                                                                                                                           cccATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTAc                       >  2:4618517/1‑150 (MQ=255)
                                                                                                                             cATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACaa                     <  2:3293722/150‑1 (MQ=255)
                                                                                                                               tatTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGa                   <  1:2633339/150‑1 (MQ=255)
                                                                                                                                  tGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCt                >  1:4414129/1‑150 (MQ=255)
                                                                                                                                    cAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAt              <  2:1007266/150‑1 (MQ=255)
                                                                                                                                     aGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAtt             <  1:4457010/150‑1 (MQ=255)
                                                                                                                                        aaGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAAc          >  1:1126912/1‑150 (MQ=255)
                                                                                                                                         aGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACg         >  2:83512/1‑150 (MQ=255)
                                                                                                                                             cATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTgg     <  1:1424799/150‑1 (MQ=255)
                                                                                                                                             cATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTgg     >  2:4216026/1‑150 (MQ=255)
                                                                                                                                              aTGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc    <  1:449804/150‑1 (MQ=255)
                                                                                                                                                ggATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgcgc  >  1:2142849/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGC  >  CP009273/2872940‑2873233

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: