Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,878,668 G→T 85.2% 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,6680GT85.2% 646.3 / 130.2 279V637V (GTC→GTAygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (14/27);  new base T (113/125);  total (127/152)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

Only 100 of 279 total aligned reads displayed.

AAAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTCCT  >  CP009273/2878522‑2878816
                                                                                                                                                  |                                                                                                                                                    
   aTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTgg                                                                                                                                                >  2:3605740/1‑150 (MQ=255)
   aTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTgg                                                                                                                                                >  1:1919128/1‑150 (MQ=255)
   aTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTgg                                                                                                                                                >  2:3439406/1‑150 (MQ=255)
   aTGGTGGAAATAGGAATCTCAAACCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTgg                                                                                                                                                <  1:1135710/150‑1 (MQ=255)
     ggtggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATTTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGt                                                                                                                                              <  2:3028733/150‑1 (MQ=255)
      gtggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGtt                                                                                                                                             >  1:4000520/1‑150 (MQ=255)
        ggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGc                                                                                                                                           >  2:2654469/1‑150 (MQ=255)
         gCAACAGGAATCTCAAAACCAGCGGGACTATTTTTGCGTTGATGGCTATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCg                                                                                                                                          <  1:1984592/150‑1 (MQ=255)
           aaCAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGAc                                                                                                                                        >  1:384741/1‑150 (MQ=255)
           aaCAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGAc                                                                                                                                        >  1:1641945/1‑150 (MQ=255)
              aGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAg                                                                                                                                     >  1:4052248/1‑150 (MQ=255)
                  aTCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATc                                                                                                                                 >  2:159874/1‑150 (MQ=255)
                        aaaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCa                                                                                                                           >  1:364341/1‑150 (MQ=255)
                         aaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCaa                                                                                                                          >  1:1656654/1‑150 (MQ=255)
                            cAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACAt                                                                                                                       <  2:2436405/150‑1 (MQ=255)
                              gCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACAttt                                                                                                                     <  2:851986/150‑1 (MQ=255)
                               cGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTGCCTGGGTTGCGACAAGTATCCGTCCAACATTTc                                                                                                                    >  1:98034/1‑150 (MQ=255)
                                 ggACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGc                                                                                                                  >  1:4384308/1‑150 (MQ=255)
                                  gACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCt                                                                                                                 <  2:3009701/150‑1 (MQ=255)
                                      atatTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTccc                                                                                                             >  1:1382166/1‑150 (MQ=255)
                                            gCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAttttt                                                                                                       <  2:141117/150‑1 (MQ=255)
                                                     gCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAAtt                                                                                              <  1:1027705/150‑1 (MQ=255)
                                                     gCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAAtt                                                                                              <  1:3836116/150‑1 (MQ=255)
                                                     gCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAAtt                                                                                              <  1:1342898/150‑1 (MQ=255)
                                                      cGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                             <  2:1714528/150‑1 (MQ=255)
                                                      cGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                             <  2:1449829/150‑1 (MQ=255)
                                                       gATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGc                                                                                            <  1:3681816/150‑1 (MQ=255)
                                                       gATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGc                                                                                            >  1:520141/1‑150 (MQ=255)
                                                             aaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaataa                                                                                      <  1:4526729/150‑1 (MQ=255)
                                                               cGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACt                                                                                    >  1:2140518/1‑150 (MQ=255)
                                                                   ccAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                                >  2:3548310/1‑150 (MQ=255)
                                                                   ccAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                                >  1:196448/1‑150 (MQ=255)
                                                                   ccAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                                <  1:3395909/150‑1 (MQ=255)
                                                                       tCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtctc                                                                            <  1:1011625/150‑1 (MQ=255)
                                                                       tCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtctc                                                                            <  1:717014/150‑1 (MQ=255)
                                                                          ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCttt                                                                         <  2:599101/150‑1 (MQ=255)
                                                                            ggAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCttttt                                                                       >  2:2411304/1‑150 (MQ=255)
                                                                             gAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTc                                                                      >  1:3942100/1‑150 (MQ=255)
                                                                               aaaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCac                                                                    >  1:4141611/1‑150 (MQ=255)
                                                                                aaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacg                                                                   >  1:2661233/1‑150 (MQ=255)
                                                                                aaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacg                                                                   <  2:1921404/150‑1 (MQ=255)
                                                                                 aGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTCATAACTCGATTCTCTTTTTCacga                                                                  >  2:1157694/1‑150 (MQ=255)
                                                                                   cAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacgacg                                                                >  2:2720268/1‑150 (MQ=255)
                                                                                     aaTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGAt                                                              >  1:1533163/1‑150 (MQ=255)
                                                                                      aTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATc                                                             <  2:247326/150‑1 (MQ=255)
                                                                                      aTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATc                                                             <  2:2354628/150‑1 (MQ=255)
                                                                                      aTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATc                                                             >  2:3590915/1‑150 (MQ=255)
                                                                                         tGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGtt                                                          >  2:2029203/1‑150 (MQ=255)
                                                                                         tGCAGGACAATACTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATGCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGtt                                                          <  2:3127976/150‑1 (MQ=255)
                                                                                          gCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTc                                                         <  2:1292639/150‑1 (MQ=255)
                                                                                           cAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCa                                                        >  2:2516409/1‑150 (MQ=255)
                                                                                            aGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAg                                                       >  1:3510555/1‑150 (MQ=255)
                                                                                               aCAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGt                                                    <  1:443984/150‑1 (MQ=255)
                                                                                                cAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTa                                                   <  2:1751837/150‑1 (MQ=255)
                                                                                                   tGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAg                                                <  2:4417670/150‑1 (MQ=255)
                                                                                                    gCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgc                                               <  2:3755907/150‑1 (MQ=255)
                                                                                                     cTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcg                                              <  2:520141/150‑1 (MQ=255)
                                                                                                     cTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcg                                              >  1:369251/1‑150 (MQ=255)
                                                                                                     cTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcg                                              <  2:4018724/150‑1 (MQ=255)
                                                                                                        aGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcgc                                           >  1:17735/1‑150 (MQ=255)
                                                                                                        aGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcgc                                           >  1:4148184/1‑150 (MQ=255)
                                                                                                        aGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcgc                                           >  1:1873846/1‑150 (MQ=255)
                                                                                                           aaTTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATg                                        <  2:818332/150‑1 (MQ=255)
                                                                                                           aaTTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATg                                        <  2:4344708/150‑1 (MQ=255)
                                                                                                               aaacactcataagcaaCGTCGAGGGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATACCTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAAc                                    <  2:3738124/144‑1 (MQ=255)
                                                                                                               aaCCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAAc                                    <  2:1068198/150‑1 (MQ=255)
                                                                                                               aaCCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAAc                                    >  2:1992094/1‑150 (MQ=255)
                                                                                                               aaCCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAAc                                    >  2:4231992/1‑150 (MQ=255)
                                                                                                                aCCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACa                                   <  1:269485/150‑1 (MQ=255)
                                                                                                                   aatcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAAt                                <  1:3603957/150‑1 (MQ=255)
                                                                                                                     tcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCt                              <  1:2366332/150‑1 (MQ=255)
                                                                                                                      caaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCta                             <  1:4251315/150‑1 (MQ=255)
                                                                                                                       aaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtat                            <  2:4138023/150‑1 (MQ=255)
                                                                                                                        aaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTACAGCGCGCATGAAACAAATCtata                           >  2:1916054/1‑150 (MQ=255)
                                                                                                                        aaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtata                           <  1:200645/150‑1 (MQ=255)
                                                                                                                         aatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTatat                          >  1:4384648/1‑150 (MQ=255)
                                                                                                                          atcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATc                         <  1:3076467/150‑1 (MQ=255)
                                                                                                                          atcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATc                         >  1:2949597/1‑150 (MQ=255)
                                                                                                                          atcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATc                         <  2:1641945/150‑1 (MQ=255)
                                                                                                                          atcaaCGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATc                         <  1:4245408/150‑1 (MQ=255)
                                                                                                                             aaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTAc                      <  2:2074013/150‑1 (MQ=255)
                                                                                                                             aaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTAc                      >  2:3686840/1‑150 (MQ=255)
                                                                                                                             aaCGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTAc                      >  2:3795256/1‑150 (MQ=255)
                                                                                                                              aCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACt                     <  1:93809/150‑1 (MQ=255)
                                                                                                                              aCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACt                     <  1:3337012/150‑1 (MQ=255)
                                                                                                                              aCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACt                     <  2:4000520/150‑1 (MQ=255)
                                                                                                                               cGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACtt                    <  2:3305351/150‑1 (MQ=255)
                                                                                                                                gTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTg                   >  1:688928/1‑150 (MQ=255)
                                                                                                                                gTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTg                   <  1:589225/150‑1 (MQ=255)
                                                                                                                                gTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTg                   >  1:3749715/1‑150 (MQ=255)
                                                                                                                                gTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTg                   <  2:4448029/150‑1 (MQ=255)
                                                                                                                                       gATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTa            >  2:4133784/1‑150 (MQ=255)
                                                                                                                                        aTTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAt           >  2:2363022/1‑150 (MQ=255)
                                                                                                                                         ttGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAtt          <  1:1807209/150‑1 (MQ=255)
                                                                                                                                          tGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAttt         <  1:1171406/150‑1 (MQ=255)
                                                                                                                                          tGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAttt         <  1:4482822/150‑1 (MQ=255)
                                                                                                                                          tGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAttt         <  1:2223527/150‑1 (MQ=255)
                                                                                                                                              ccACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCt     <  1:2536546/150‑1 (MQ=255)
                                                                                                                                              ccACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCt     <  2:135309/150‑1 (MQ=255)
                                                                                                                                                 cGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTCCt  <  1:2848093/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTCCT  >  CP009273/2878522‑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: