Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,878,668 G→T 92.6% 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,6680GT92.6% 573.4 / 35.9 190V637V (GTC→GTAygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (8/6);  new base T (99/77);  total (107/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 190 total aligned reads displayed.

CAAAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTCCT  >  CP009273/2878521‑2878816
                                                                                                                                                   |                                                                                                                                                    
 aaaaTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTAcc                                                                                                                                                   <  2:3271213/150‑1 (MQ=255)
  aaaTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCt                                                                                                                                                  >  2:1147662/1‑150 (MQ=255)
       gtggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGtt                                                                                                                                             <  2:2281504/150‑1 (MQ=255)
          gCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCg                                                                                                                                          <  2:2916549/150‑1 (MQ=255)
          gCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCg                                                                                                                                          <  1:301984/150‑1 (MQ=255)
                ggAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGt                                                                                                                                    <  1:2490908/150‑1 (MQ=255)
                   aTCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATc                                                                                                                                 <  1:2627821/150‑1 (MQ=255)
                     ctcAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCg                                                                                                                               >  2:2772476/1‑150 (MQ=255)
                        aaaaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTcc                                                                                                                            <  2:353491/150‑1 (MQ=255)
                        aaaaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTcc                                                                                                                            <  2:2527085/150‑1 (MQ=255)
                          aaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCaa                                                                                                                          >  1:484527/1‑150 (MQ=255)
                             cAGCGGGACTATATTTGCGATGATGGCTAGGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGAAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACAt                                                                                                                       <  1:2369188/150‑1 (MQ=255)
                                cGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTc                                                                                                                    >  1:3807136/1‑150 (MQ=255)
                                cGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTc                                                                                                                    <  2:2205611/150‑1 (MQ=255)
                                cGGGACGATATTGGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTc                                                                                                                    <  1:83657/150‑1 (MQ=255)
                                 gggACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCg                                                                                                                   >  2:1122469/1‑150 (MQ=255)
                                   gACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCt                                                                                                                 <  1:2688310/150‑1 (MQ=255)
                                   gACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCt                                                                                                                 <  1:4216148/150‑1 (MQ=255)
                                   gACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCt                                                                                                                 <  2:2480060/150‑1 (MQ=255)
                                     ggATTTTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTc                                                                                                               <  1:737130/149‑1 (MQ=255)
                                     cGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTc                                                                                                               <  2:547766/150‑1 (MQ=255)
                                      gATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTcc                                                                                                              >  1:2492124/1‑150 (MQ=255)
                                      gATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTcc                                                                                                              >  2:3123466/1‑150 (MQ=255)
                                      gATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTcc                                                                                                              >  2:3422665/1‑150 (MQ=255)
                                      gATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTcc                                                                                                              >  2:1405902/1‑150 (MQ=255)
                                         atTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAt                                                                                                           >  2:3803140/1‑150 (MQ=255)
                                           ttGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAttt                                                                                                         <  2:2999712/150‑1 (MQ=255)
                                            tGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAtttt                                                                                                        <  1:28757/150‑1 (MQ=255)
                                             gCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAttttt                                                                                                       <  2:1748725/150‑1 (MQ=255)
                                                  gatgGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGa                                                                                                  <  2:1087534/150‑1 (MQ=255)
                                                   atgGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCTCTTCCCATTTTTGACGaa                                                                                                 >  1:3161640/1‑150 (MQ=255)
                                                    tgGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGaaa                                                                                                >  1:3520595/1‑150 (MQ=255)
                                                    tgGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGaaa                                                                                                >  2:1585431/1‑150 (MQ=255)
                                                       cGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                             >  1:2517918/1‑150 (MQ=255)
                                                       cGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                             >  1:477467/1‑150 (MQ=255)
                                                          tttAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCta                                                                                          <  1:1161697/148‑1 (MQ=255)
                                                           gTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGACATTGCtaa                                                                                         >  1:3051349/1‑150 (MQ=255)
                                                            tAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaat                                                                                        <  1:1435310/150‑1 (MQ=255)
                                                                    ccAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                                >  2:3439262/1‑150 (MQ=255)
                                                                    ccAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                                >  2:1579195/1‑150 (MQ=255)
                                                                     cAATCGTTTGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAtt                                                                               <  2:1272162/150‑1 (MQ=255)
                                                                     cAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAtt                                                                               >  1:274294/1‑150 (MQ=255)
                                                                      aaTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtc                                                                              <  2:1486528/150‑1 (MQ=255)
                                                                      aaTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtc                                                                              >  2:1810341/1‑150 (MQ=255)
                                                                      aaTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtc                                                                              <  1:501561/150‑1 (MQ=255)
                                                                      aaTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTAGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtc                                                                              >  2:3503543/1‑150 (MQ=255)
                                                                       aTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtct                                                                             >  2:3863026/1‑150 (MQ=255)
                                                                         cGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGACATTGCTAATAACTCGATTctct                                                                           >  2:2163373/1‑150 (MQ=255)
                                                                           ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCttt                                                                         <  1:2712596/150‑1 (MQ=255)
                                                                           ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCttt                                                                         <  1:541525/150‑1 (MQ=255)
                                                                             ggAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCttttt                                                                       >  1:2531351/1‑150 (MQ=255)
                                                                               aaaaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCa                                                                     >  1:1220926/1‑150 (MQ=255)
                                                                                 aaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacg                                                                   <  1:139174/150‑1 (MQ=255)
                                                                                    cAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacgacg                                                                <  2:3659480/150‑1 (MQ=255)
                                                                                    cAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacgacg                                                                <  2:2931496/150‑1 (MQ=255)
                                                                                      aaTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGAt                                                              <  2:3266791/150‑1 (MQ=255)
                                                                                       aTCTGCAGGACAATGCTGAGTAATTACCCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATc                                                             <  2:2140396/150‑1 (MQ=255)
                                                                                       aTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATc                                                             <  2:1009818/150‑1 (MQ=255)
                                                                                        tCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCg                                                            >  2:1908235/1‑150 (MQ=255)
                                                                                         cTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGt                                                           >  1:2390989/1‑150 (MQ=255)
                                                                                         cTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGt                                                           >  2:455517/1‑150 (MQ=255)
                                                                                         cTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGt                                                           >  2:455583/1‑150 (MQ=255)
                                                                                             aGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAg                                                       >  2:39848/1‑150 (MQ=255)
                                                                                             aGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAg                                                       >  2:2897424/1‑150 (MQ=255)
                                                                                               gACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCg                                                     >  1:4294464/1‑150 (MQ=255)
                                                                                                aCAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATAGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGt                                                    >  1:2171256/1‑150 (MQ=255)
                                                                                                     gCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgc                                               >  2:1234587/1‑150 (MQ=255)
                                                                                                       tGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgc                                             >  2:2445451/1‑150 (MQ=255)
                                                                                                        gAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcg                                            >  1:1850991/1‑150 (MQ=255)
                                                                                                        gAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcg                                            >  1:951268/1‑150 (MQ=255)
                                                                                                         aGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcgc                                           >  1:2271506/1‑150 (MQ=255)
                                                                                                         aGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcgc                                           >  1:3994387/1‑150 (MQ=255)
                                                                                                            aaTTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATg                                        >  2:4261164/1‑150 (MQ=255)
                                                                                                             aTTAACCAATCAAAATCAACGTCGAGTGATTGTCCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGa                                       >  2:3214218/1‑150 (MQ=255)
                                                                                                               tAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGaaa                                     <  1:2228478/150‑1 (MQ=255)
                                                                                                                aaCCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAAc                                    <  2:274294/150‑1 (MQ=255)
                                                                                                                    aatcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAAt                                <  1:1147662/150‑1 (MQ=255)
                                                                                                                    aatcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAAt                                >  2:3993396/1‑150 (MQ=255)
                                                                                                                       caaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCta                             <  2:4219782/150‑1 (MQ=255)
                                                                                                                         aaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtata                           <  2:484527/150‑1 (MQ=255)
                                                                                                                           atcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATc                         >  2:847283/1‑150 (MQ=255)
                                                                                                                              aaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTAc                      >  1:172503/1‑150 (MQ=255)
                                                                                                                                cGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACtt                    <  1:4039061/150‑1 (MQ=255)
                                                                                                                                 gTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTg                   <  2:3931801/150‑1 (MQ=255)
                                                                                                                                   cGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCg                 >  2:3620155/1‑150 (MQ=255)
                                                                                                                                   cGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCg                 >  1:882212/1‑150 (MQ=255)
                                                                                                                                    gAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCgt                >  1:4117517/1‑150 (MQ=255)
                                                                                                                                      gTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCgtgt              >  1:19549/1‑150 (MQ=255)
                                                                                                                                       tGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCTGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGtt             <  1:3803140/150‑1 (MQ=255)
                                                                                                                                       tGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGtt             >  2:3087418/1‑150 (MQ=255)
                                                                                                                                       tGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGtt             >  2:2720544/1‑150 (MQ=255)
                                                                                                                                        gATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTa            <  1:3113685/150‑1 (MQ=255)
                                                                                                                                        gATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTa            >  2:2309552/1‑150 (MQ=255)
                                                                                                                                         aTTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAt           <  1:290735/150‑1 (MQ=255)
                                                                                                                                         aTTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAt           <  1:3886779/150‑1 (MQ=255)
                                                                                                                                         aTTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAt           <  1:628982/150‑1 (MQ=255)
                                                                                                                                         aTTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAt           <  1:2014262/150‑1 (MQ=255)
                                                                                                                                           tGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAttt         <  1:4136528/150‑1 (MQ=255)
                                                                                                                                              tCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGc      <  2:1673950/150‑1 (MQ=255)
                                                                                                                                                  cGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTCCt  >  1:3248443/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAAAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTCCT  >  CP009273/2878521‑2878816

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: