Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,878,098 T→A 6.9% A827A (GCA→GCT 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,0980TA 6.9% 812.8 / 37.9 260A827A (GCA→GCTygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base T (116/126);  new base A (11/7);  total (127/133)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 260 total aligned reads displayed.

ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTA  >  CP009273/2877951‑2878245
                                                                                                                                                   |                                                                                                                                                   
 tCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGc                                                                                                                                                   >  2:3814480/1‑149 (MQ=255)
   cccTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCaa                                                                                                                                                 <  2:2304005/149‑1 (MQ=255)
       gTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAgcgc                                                                                                                                             >  2:1344372/1‑149 (MQ=255)
            ggTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCAt                                                                                                                                        <  2:3982527/149‑1 (MQ=255)
            ggTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGAGTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCAt                                                                                                                                        <  2:694048/149‑1 (MQ=255)
              taataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATAc                                                                                                                                      >  1:249082/1‑149 (MQ=255)
               aataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATAct                                                                                                                                     >  2:424306/1‑149 (MQ=255)
                ataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATActg                                                                                                                                    <  2:1273824/149‑1 (MQ=255)
                 taaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATActgc                                                                                                                                   <  2:2941416/149‑1 (MQ=255)
                  aaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCATAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATActgct                                                                                                                                  <  2:3604131/149‑1 (MQ=255)
                      aTACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTc                                                                                                                              >  1:2763945/1‑149 (MQ=255)
                      aTACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTc                                                                                                                              >  2:1843386/1‑149 (MQ=255)
                      aTACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTc                                                                                                                              >  2:4109196/1‑149 (MQ=255)
                       tACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCa                                                                                                                             <  2:246069/149‑1 (MQ=255)
                       tACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTAAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCa                                                                                                                             <  1:489839/149‑1 (MQ=255)
                         cTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATg                                                                                                                           >  1:2265154/1‑149 (MQ=255)
                          tGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAGGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGa                                                                                                                          >  1:2823947/1‑149 (MQ=255)
                          tGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGa                                                                                                                          >  1:1840540/1‑149 (MQ=255)
                           gTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGAc                                                                                                                         >  1:1607161/1‑149 (MQ=255)
                           gTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGAc                                                                                                                         <  2:3327219/149‑1 (MQ=255)
                              aCCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTa                                                                                                                      <  2:641173/149‑1 (MQ=255)
                                ccTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAgg                                                                                                                    <  1:2026495/149‑1 (MQ=255)
                                    ccAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCt                                                                                                                <  2:1331589/149‑1 (MQ=255)
                                    ccAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCt                                                                                                                <  2:1671873/149‑1 (MQ=255)
                                       gACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGt                                                                                                             >  1:197072/1‑149 (MQ=255)
                                            atccatccAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGAcc                                                                                                        <  1:2175372/149‑1 (MQ=255)
                                             tccatccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCt                                                                                                       >  2:228956/1‑149 (MQ=255)
                                              ccatccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTg                                                                                                      <  2:450274/149‑1 (MQ=255)
                                                atccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGc                                                                                                    >  2:1763120/1‑149 (MQ=255)
                                                 tccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGcc                                                                                                   <  2:569593/149‑1 (MQ=255)
                                                 tccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGcc                                                                                                   >  1:3776650/1‑149 (MQ=255)
                                                 tccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGcc                                                                                                   >  1:1733912/1‑149 (MQ=255)
                                                 tccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCAACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAAATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGcc                                                                                                   >  1:3796002/1‑149 (MQ=255)
                                                    aGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATc                                                                                                <  1:1987186/149‑1 (MQ=255)
                                                    aGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATc                                                                                                >  2:2193508/1‑149 (MQ=255)
                                                     gATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCg                                                                                               >  2:2204908/1‑149 (MQ=255)
                                                     gATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCg                                                                                               >  1:4071791/1‑149 (MQ=255)
                                                       ttCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAg                                                                                             >  1:553952/1‑149 (MQ=255)
                                                          tGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAg                                                                                          >  1:3294673/1‑149 (MQ=255)
                                                          tGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAg                                                                                          >  1:1241097/1‑149 (MQ=255)
                                                            tttCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAgtt                                                                                        <  2:4095817/149‑1 (MQ=255)
                                                            tttCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAgtt                                                                                        >  1:2096687/1‑149 (MQ=255)
                                                            tttCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTGTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAgtt                                                                                        <  2:769542/149‑1 (MQ=255)
                                                              tCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAgttgt                                                                                      >  2:2887924/1‑149 (MQ=255)
                                                               cccTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAgttgtt                                                                                     <  2:3197628/149‑1 (MQ=255)
                                                               cccTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAgttgtt                                                                                     >  1:1967192/1‑149 (MQ=255)
                                                                ccTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGttt                                                                                    <  2:1258679/149‑1 (MQ=255)
                                                                   tCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATAAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTAcc                                                                                 <  1:1776546/149‑1 (MQ=255)
                                                                      aGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGa                                                                              >  1:3015112/1‑149 (MQ=255)
                                                                         cAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGa                                                                           <  1:3648541/149‑1 (MQ=255)
                                                                         cAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGa                                                                           <  1:2666212/149‑1 (MQ=255)
                                                                             gTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGtt                                                                       <  1:4027765/149‑1 (MQ=255)
                                                                              tAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGttt                                                                      <  1:3989432/149‑1 (MQ=255)
                                                                               aaCCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTg                                                                     <  2:2994563/149‑1 (MQ=255)
                                                                                aCCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGt                                                                    <  1:2618954/149‑1 (MQ=255)
                                                                                 cccATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTa                                                                   <  1:1130955/149‑1 (MQ=255)
                                                                                     tCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATa                                                               >  1:832261/1‑149 (MQ=255)
                                                                                     tCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATa                                                               >  1:2550257/1‑149 (MQ=255)
                                                                                      cTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATaa                                                              <  2:1012829/149‑1 (MQ=255)
                                                                                       ttCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAg                                                             <  2:2823850/149‑1 (MQ=255)
                                                                                       ttCATCTACTACTTCAGAAAAACTACGTGTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAg                                                             <  1:3548841/149‑1 (MQ=255)
                                                                                         cATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGc                                                           <  1:2474030/149‑1 (MQ=255)
                                                                                          aTCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCa                                                          >  2:1722770/1‑149 (MQ=255)
                                                                                           tCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaa                                                         <  1:99621/149‑1 (MQ=255)
                                                                                           tCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaa                                                         <  1:3405826/149‑1 (MQ=255)
                                                                                            ctactaCTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaat                                                        <  1:4109196/149‑1 (MQ=255)
                                                                                              actactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataa                                                      >  1:2137986/1‑149 (MQ=255)
                                                                                               ctactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataat                                                     <  2:1534986/149‑1 (MQ=255)
                                                                                               ctactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataat                                                     <  1:2286768/149‑1 (MQ=255)
                                                                                               ctactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataat                                                     <  1:1258303/149‑1 (MQ=255)
                                                                                               ctactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataat                                                     <  1:3221700/149‑1 (MQ=255)
                                                                                                 actTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATgg                                                   >  1:2339855/1‑149 (MQ=255)
                                                                                                  ctTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGc                                                  >  1:2724756/1‑149 (MQ=255)
                                                                                                   ttCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCa                                                 <  2:3941872/149‑1 (MQ=255)
                                                                                                    tCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAg                                                <  1:3465848/149‑1 (MQ=255)
                                                                                                       gAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCt                                             <  1:418098/149‑1 (MQ=255)
                                                                                                        aaaaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTc                                            <  1:1359696/149‑1 (MQ=255)
                                                                                                        aaaaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTc                                            <  2:1579739/149‑1 (MQ=255)
                                                                                                         aaaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCa                                           >  2:551713/1‑149 (MQ=255)
                                                                                                         aaaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCCCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCa                                           <  2:1506689/149‑1 (MQ=255)
                                                                                                          aaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCAt                                          >  2:242456/1‑149 (MQ=255)
                                                                                                            aCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCAttt                                        >  1:1999796/1‑149 (MQ=255)
                                                                                                               aCGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTccc                                     >  2:3271303/1‑149 (MQ=255)
                                                                                                                 gTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCa                                   <  1:3179069/149‑1 (MQ=255)
                                                                                                                   ttCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATc                                 <  2:1990379/149‑1 (MQ=255)
                                                                                                                      cAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCt                              <  1:1079459/149‑1 (MQ=255)
                                                                                                                       aGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTc                             >  1:2225402/1‑149 (MQ=255)
                                                                                                                            aaGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAc                        <  2:1150864/149‑1 (MQ=255)
                                                                                                                             aGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAcc                       >  2:316987/1‑149 (MQ=255)
                                                                                                                             aGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAcc                       >  2:1903418/1‑149 (MQ=255)
                                                                                                                              gggTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCg                      >  2:2480399/1‑149 (MQ=255)
                                                                                                                                 tACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCaa                   >  1:3721168/1‑149 (MQ=255)
                                                                                                                                     ttGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAAt               >  2:1419690/1‑149 (MQ=255)
                                                                                                                                      tGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATg              <  1:1569424/149‑1 (MQ=255)
                                                                                                                                          gcgATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGttt          <  2:1614427/149‑1 (MQ=255)
                                                                                                                                          gcgATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGttt          >  2:826945/1‑149 (MQ=255)
                                                                                                                                           cgATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATTATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTc         <  1:31384/149‑1 (MQ=255)
                                                                                                                                                aaGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGt    <  1:951451/149‑1 (MQ=255)
                                                                                                                                                  gTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTa  <  1:624390/149‑1 (MQ=255)
                                                                                                                                                  gTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTa  >  2:888452/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTA  >  CP009273/2877951‑2878245

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: