Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,878,668 G→T 100% V637V (GTC→GTA)  ygcB ← R‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,878,6680GT99.3% 1020.0 / ‑1.8 309V637V (GTC→GTA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (138/168);  total (140/169)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 309 total aligned reads displayed.

AAAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAG  >  CP009273/2878522‑2878811
                                                                                                                                                  |                                                                                                                                               
aaaaTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTAc                                                                                                                                               <  2:1581478/149‑1 (MQ=255)
   aTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTg                                                                                                                                            >  2:1487705/1‑149 (MQ=255)
   aTGGGGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCCACGTCGAGTGATTGTTCCACTACCTg                                                                                                                                            <  2:3216805/149‑1 (MQ=255)
    tggtggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTgg                                                                                                                                           >  2:2973992/1‑149 (MQ=255)
       tggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGtt                                                                                                                                        <  2:326170/149‑1 (MQ=255)
         gCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGc                                                                                                                                      <  1:3027293/149‑1 (MQ=255)
         gCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGc                                                                                                                                      <  1:2910083/149‑1 (MQ=255)
           aaCAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGa                                                                                                                                    >  1:317984/1‑149 (MQ=255)
            aCAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGAc                                                                                                                                   <  2:1635507/149‑1 (MQ=255)
             cAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACa                                                                                                                                  <  1:4595886/149‑1 (MQ=255)
                gAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGt                                                                                                                               <  1:108797/149‑1 (MQ=255)
                 aaTCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTa                                                                                                                              <  1:334525/149‑1 (MQ=255)
                  aTCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTAt                                                                                                                             >  2:3916665/1‑149 (MQ=255)
                   tctcAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATc                                                                                                                            <  1:2117061/149‑1 (MQ=255)
                     tcAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCg                                                                                                                          <  2:1010476/149‑1 (MQ=255)
                          aCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCaa                                                                                                                     >  1:147772/1‑149 (MQ=255)
                          aCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCaa                                                                                                                     >  1:3456510/1‑149 (MQ=255)
                             aGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCTAGTGATTGTTCCACTACCAGGGTTGCGACAAGTATCCGTCCAACAt                                                                                                                  >  1:1403595/1‑149 (MQ=255)
                             aGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACAt                                                                                                                  <  1:1020049/149‑1 (MQ=255)
                                 ggACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCg                                                                                                              >  1:2130562/1‑149 (MQ=255)
                                  gACGATATGTGCGATTATGGCGATTTAAACGGCCCAATCGTTTGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGc                                                                                                             <  2:2411989/149‑1 (MQ=255)
                                   aCGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCt                                                                                                            >  2:3802679/1‑149 (MQ=255)
                                     gATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTc                                                                                                          >  1:355445/1‑149 (MQ=255)
                                     gATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTc                                                                                                          <  2:4086909/149‑1 (MQ=255)
                                       tttttGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTccc                                                                                                        <  1:2234000/147‑1 (MQ=255)
                                         tttGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAt                                                                                                      <  1:1787129/149‑1 (MQ=255)
                                            gCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAtttt                                                                                                   <  1:3775964/149‑1 (MQ=255)
                                            gCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAtttt                                                                                                   <  1:3477564/149‑1 (MQ=255)
                                                     gCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAAt                                                                                          <  2:3863631/149‑1 (MQ=255)
                                                      cGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAAtt                                                                                         <  2:3557149/149‑1 (MQ=255)
                                                      cGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAAtt                                                                                         <  1:3101097/149‑1 (MQ=255)
                                                        aTGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGCGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGc                                                                                       >  2:4322662/1‑149 (MQ=255)
                                                        aTGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGc                                                                                       >  2:4010150/1‑149 (MQ=255)
                                                         tGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCt                                                                                      >  1:3586676/1‑149 (MQ=255)
                                                         tGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCt                                                                                      >  1:579367/1‑149 (MQ=255)
                                                          gTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCta                                                                                     >  2:2044117/1‑149 (MQ=255)
                                                           tAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaa                                                                                    >  1:1676759/1‑149 (MQ=255)
                                                             aaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaata                                                                                  >  2:1721025/1‑149 (MQ=255)
                                                             aaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaata                                                                                  <  1:4449056/149‑1 (MQ=255)
                                                              aCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaataa                                                                                 <  2:606084/149‑1 (MQ=255)
                                                                   ccAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGa                                                                            <  1:4309141/149‑1 (MQ=255)
                                                                    cAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                           <  2:4356688/149‑1 (MQ=255)
                                                                    cAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                           <  2:1967239/149‑1 (MQ=255)
                                                                    cAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                           >  2:1284351/1‑149 (MQ=255)
                                                                      aTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtc                                                                         >  1:3938933/1‑149 (MQ=255)
                                                                      aTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtc                                                                         <  1:4430329/149‑1 (MQ=255)
                                                                      aTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtc                                                                         <  1:1759102/149‑1 (MQ=255)
                                                                        cGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtctc                                                                       <  1:1296591/149‑1 (MQ=255)
                                                                          ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCtt                                                                     >  2:2772209/1‑149 (MQ=255)
                                                                           tGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCttt                                                                    <  2:4781057/149‑1 (MQ=255)
                                                                             gAAAAGCAAATCTTCAGGACAATGATGAGTAATTAACAAATCAAAATCAACTTCGAGTGATTGTTACAATACCTGGTTTGCTACAAGTATCAGTCCAACATTTCTCTTCCCATTTTTGCCGAAATTGCTAATTACTCGATTCTCttttt                                                                  <  2:2245352/149‑1 (MQ=255)
                                                                               aaatcaaatcTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCa                                                                <  1:4282248/149‑1 (MQ=255)
                                                                                aaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCac                                                               <  2:3031889/149‑1 (MQ=255)
                                                                                      aTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGAt                                                         <  2:4779215/149‑1 (MQ=255)
                                                                                      aTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGAt                                                         <  1:186187/149‑1 (MQ=255)
                                                                                           cAGGACAATGCTGAGTAATTAACCAATCAAAATCACCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTc                                                    <  1:2339716/149‑1 (MQ=255)
                                                                                           cAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTc                                                    >  2:1922725/1‑149 (MQ=255)
                                                                                           cAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTc                                                    <  2:464683/149‑1 (MQ=255)
                                                                                           cAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTc                                                    <  1:4765148/149‑1 (MQ=255)
                                                                                             ggACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAg                                                  >  2:3576485/1‑149 (MQ=255)
                                                                                              gACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGc                                                 <  1:3802679/149‑1 (MQ=255)
                                                                                               aCAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCg                                                >  2:44986/1‑149 (MQ=255)
                                                                                                 aaTGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTa                                              >  2:3955509/1‑149 (MQ=255)
                                                                                                  aTGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTaa                                             >  2:1805608/1‑149 (MQ=255)
                                                                                                  aTGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTaa                                             >  1:4387246/1‑149 (MQ=255)
                                                                                                    gCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAg                                           <  2:2378372/149‑1 (MQ=255)
                                                                                                     cTGAGTAATTAACCAATCAAAATCAACGTCGATTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTACAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgc                                          <  2:1403595/149‑1 (MQ=255)
                                                                                                     cTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgc                                          <  2:237247/149‑1 (MQ=255)
                                                                                                       gAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgc                                        >  2:3178993/1‑149 (MQ=255)
                                                                                                        aGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcg                                       >  1:519930/1‑149 (MQ=255)
                                                                                                         gTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcgc                                      <  2:3085213/149‑1 (MQ=255)
                                                                                                         gTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcgc                                      >  2:3419313/1‑149 (MQ=255)
                                                                                                          tAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCa                                     <  2:4156182/149‑1 (MQ=255)
                                                                                                          tAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCa                                     <  2:1766265/149‑1 (MQ=255)
                                                                                                          tAATTAACCAATCAAAATCAACGACGAGTGATTGTTCAACTACCTGGGTTGCGACAAGTATCCGTCCAACAGTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCTTAAAGCGCGCa                                     >  1:612917/1‑149 (MQ=255)
                                                                                                            aTTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATg                                   >  2:7628/1‑149 (MQ=255)
                                                                                                             ttAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGa                                  >  1:86725/1‑149 (MQ=255)
                                                                                                                 ccAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACa                              <  2:3666329/149‑1 (MQ=255)
                                                                                                                 ccAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACa                              <  2:1285893/149‑1 (MQ=255)
                                                                                                                   aatcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACaaa                            <  2:595548/149‑1 (MQ=255)
                                                                                                                   aatcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACaaa                            >  2:4439240/1‑149 (MQ=255)
                                                                                                                   aatcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACaaa                            >  2:4024932/1‑149 (MQ=255)
                                                                                                                     tcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATc                          <  2:644472/149‑1 (MQ=255)
                                                                                                                      caaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCt                         >  2:1146278/1‑149 (MQ=255)
                                                                                                                       aaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCta                        >  1:2629861/1‑149 (MQ=255)
                                                                                                                       aaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCta                        <  2:1066510/149‑1 (MQ=255)
                                                                                                                        aaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtat                       <  1:1613377/149‑1 (MQ=255)
                                                                                                                           tcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATc                    <  1:2772209/149‑1 (MQ=255)
                                                                                                                           tcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATc                    >  1:1664761/1‑149 (MQ=255)
                                                                                                                            caaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCt                   <  1:789318/149‑1 (MQ=255)
                                                                                                                             aaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTa                  <  1:4145974/149‑1 (MQ=255)
                                                                                                                             aaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTa                  <  2:2603841/149‑1 (MQ=255)
                                                                                                                                   gAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCg            >  2:3709919/1‑149 (MQ=255)
                                                                                                                                     gTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCgtg          >  1:2559110/1‑149 (MQ=255)
                                                                                                                                      tGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCgtgt         <  1:3400273/149‑1 (MQ=255)
                                                                                                                                         ttGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCGGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAt      >  2:272415/1‑149 (MQ=255)
                                                                                                                                          tGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAtt     <  1:3897/149‑1 (MQ=255)
                                                                                                                                           gTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAttt    <  1:626306/149‑1 (MQ=255)
                                                                                                                                           gTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAttt    >  2:1904396/1‑149 (MQ=255)
                                                                                                                                            ttCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTTGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTa   <  1:3707076/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AAAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAG  >  CP009273/2878522‑2878811

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: ‑