Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,878,098 T→A 6.5% 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.5% 749.1 / 30.4 243A827A (GCA→GCTygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base T (113/114);  new base A (8/8);  total (121/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTA  >  CP009273/2877952‑2878245
                                                                                                                                                  |                                                                                                                                                   
   ccTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCaa                                                                                                                                                 <  1:3502868/148‑1 (MQ=255)
      gTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAgcg                                                                                                                                              >  1:3086286/1‑148 (MQ=255)
      gTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAgcg                                                                                                                                              <  2:2329157/148‑1 (MQ=255)
      gTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAgcg                                                                                                                                              >  2:2503854/1‑148 (MQ=255)
        ataGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGcc                                                                                                                                            <  2:1154410/148‑1 (MQ=255)
          aGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTc                                                                                                                                          <  1:1080614/148‑1 (MQ=255)
             taataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATa                                                                                                                                       <  1:2734833/148‑1 (MQ=255)
             taataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATa                                                                                                                                       >  1:3401187/1‑148 (MQ=255)
              aataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATAc                                                                                                                                      >  2:2126468/1‑148 (MQ=255)
              aataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATAc                                                                                                                                      >  1:532447/1‑148 (MQ=255)
                 aaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATActgc                                                                                                                                   <  1:751883/148‑1 (MQ=255)
                 aaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATActgc                                                                                                                                   >  1:3351548/1‑148 (MQ=255)
                    aaTACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGt                                                                                                                                >  1:497966/1‑148 (MQ=255)
                    aaTACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGt                                                                                                                                <  2:2438400/148‑1 (MQ=255)
                     aTACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGtt                                                                                                                               >  2:613969/1‑148 (MQ=255)
                      tACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTc                                                                                                                              >  2:2083589/1‑148 (MQ=255)
                      tACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTc                                                                                                                              <  1:2062972/148‑1 (MQ=255)
                      tACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTc                                                                                                                              >  1:566363/1‑148 (MQ=255)
                       aCTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTAAGAAAAACTACGTTTCCAGGTGAAGGTTAAATTGACGCGATTAAGAGAAAGCGCCTCATACTGCTGTTCa                                                                                                                             <  1:24051/148‑1 (MQ=255)
                          gTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGa                                                                                                                          >  1:3172921/1‑148 (MQ=255)
                          gTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGGACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGa                                                                                                                          <  2:630654/148‑1 (MQ=255)
                              cccTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTa                                                                                                                      <  2:245079/148‑1 (MQ=255)
                               ccTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAg                                                                                                                     >  1:800781/1‑148 (MQ=255)
                                  gCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTc                                                                                                                  <  2:1837551/148‑1 (MQ=255)
                                   ccAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTcc                                                                                                                 <  2:362782/148‑1 (MQ=255)
                                           atccatccAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGAc                                                                                                         <  1:2675634/148‑1 (MQ=255)
                                               atccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTgg                                                                                                     >  1:2559446/1‑148 (MQ=255)
                                               atccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTgg                                                                                                     <  2:456555/148‑1 (MQ=255)
                                                    gATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATc                                                                                                >  1:1697827/1‑148 (MQ=255)
                                                     aTTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCg                                                                                               >  2:1314908/1‑148 (MQ=255)
                                                       tCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAg                                                                                             <  2:1480541/148‑1 (MQ=255)
                                                       tCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAg                                                                                             >  2:2892306/1‑148 (MQ=255)
                                                       tCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTGTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAg                                                                                             <  1:2016287/148‑1 (MQ=255)
                                                        cTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGc                                                                                            <  1:910450/148‑1 (MQ=255)
                                                         tGTTTCCCTTCCAGCCAAAGTAACCCCTCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCa                                                                                           <  1:2155177/148‑1 (MQ=255)
                                                         tGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCa                                                                                           <  1:2155136/148‑1 (MQ=255)
                                                                 ttCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTa                                                                                   >  1:1084405/1‑148 (MQ=255)
                                                                 ttCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTa                                                                                   >  1:1785364/1‑148 (MQ=255)
                                                                 ttCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTa                                                                                   <  1:1906845/148‑1 (MQ=255)
                                                                 ttCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTa                                                                                   <  2:1894566/148‑1 (MQ=255)
                                                                   ccAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTAcc                                                                                 <  2:1771294/148‑1 (MQ=255)
                                                                      gCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGa                                                                              <  1:3522814/148‑1 (MQ=255)
                                                                        cAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAg                                                                            <  1:3446499/148‑1 (MQ=255)
                                                                         aaaGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGa                                                                           >  2:2044309/1‑148 (MQ=255)
                                                                          aaGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGAc                                                                          >  2:1401528/1‑148 (MQ=255)
                                                                          aaGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGAc                                                                          <  2:3308672/148‑1 (MQ=255)
                                                                          aaGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGAc                                                                          <  2:837069/148‑1 (MQ=255)
                                                                             tAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGtt                                                                       >  2:880382/1‑148 (MQ=255)
                                                                             tAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGtt                                                                       >  1:2810889/1‑148 (MQ=255)
                                                                                 ccATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTGTGTa                                                                   <  1:298360/148‑1 (MQ=255)
                                                                                  cATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTAc                                                                  <  1:1589354/148‑1 (MQ=255)
                                                                                   aTCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGTCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACa                                                                 <  1:249512/148‑1 (MQ=255)
                                                                                   aTCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACa                                                                 <  1:1282522/148‑1 (MQ=255)
                                                                                    tCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACAt                                                                >  2:1215875/1‑148 (MQ=255)
                                                                                      ttCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCTTCGAGCAGTTGTTTACCTGAAGACGTTTGTACATaa                                                              <  2:255937/148‑1 (MQ=255)
                                                                                      ttCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATaa                                                              <  1:1719618/148‑1 (MQ=255)
                                                                                           ctactaCTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaa                                                         >  1:2905003/1‑148 (MQ=255)
                                                                                             actactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaata                                                       <  1:46701/148‑1 (MQ=255)
                                                                                             actactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaata                                                       >  2:518946/1‑148 (MQ=255)
                                                                                              ctactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataa                                                      <  2:3162175/148‑1 (MQ=255)
                                                                                               tactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataat                                                     <  2:681273/148‑1 (MQ=255)
                                                                                                 ctTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATgg                                                   >  1:1013745/1‑148 (MQ=255)
                                                                                                    cAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAg                                                >  1:560044/1‑148 (MQ=255)
                                                                                                     aGAAAAACTACGTTTCCAGGTGAAGGGTATATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAgg                                               >  2:1957991/1‑148 (MQ=255)
                                                                                                     aGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAgg                                               >  1:1233040/1‑148 (MQ=255)
                                                                                                      gAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGc                                              >  1:323568/1‑148 (MQ=255)
                                                                                                         aaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCa                                           >  1:1161422/1‑148 (MQ=255)
                                                                                                         aaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCa                                           <  2:497925/148‑1 (MQ=255)
                                                                                                         aaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCa                                           >  2:3407130/1‑148 (MQ=255)
                                                                                                           aCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCAtt                                         <  2:1156366/148‑1 (MQ=255)
                                                                                                            cTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCAttt                                        >  1:3115042/1‑148 (MQ=255)
                                                                                                            cTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCAttt                                        >  2:1016546/1‑148 (MQ=255)
                                                                                                             tACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCTTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTc                                       >  1:2855899/1‑148 (MQ=255)
                                                                                                                    ccAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATcc                                <  1:1501860/148‑1 (MQ=255)
                                                                                                                       ggTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTc                             >  2:109912/1‑148 (MQ=255)
                                                                                                                        gTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCg                            >  1:3469734/1‑148 (MQ=255)
                                                                                                                        gTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCg                            >  2:476522/1‑148 (MQ=255)
                                                                                                                        gTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCg                            >  2:2263342/1‑148 (MQ=255)
                                                                                                                        gTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCg                            >  2:3369918/1‑148 (MQ=255)
                                                                                                                        gTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCg                            >  2:1711854/1‑148 (MQ=255)
                                                                                                                        gTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCg                            >  1:2783047/1‑148 (MQ=255)
                                                                                                                         tGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGACCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGt                           >  2:580441/1‑148 (MQ=255)
                                                                                                                         tGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGACCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGt                           >  2:580327/1‑148 (MQ=255)
                                                                                                                          gAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGtt                          >  2:2611712/1‑148 (MQ=255)
                                                                                                                            aGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAc                        <  2:2504110/148‑1 (MQ=255)
                                                                                                                              ggTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCg                      >  2:1026371/1‑148 (MQ=255)
                                                                                                                               gTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGc                     >  2:3079514/1‑148 (MQ=255)
                                                                                                                                tACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCa                    >  1:2627092/1‑148 (MQ=255)
                                                                                                                                tACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCa                    <  1:656353/148‑1 (MQ=255)
                                                                                                                                 aCATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCaa                   >  2:2880821/1‑148 (MQ=255)
                                                                                                                                    ttGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGaa                >  1:2888600/1‑148 (MQ=255)
                                                                                                                                    ttGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATATCCCCATCCCTCGTTACCGCAAGaa                >  1:883484/1‑148 (MQ=255)
                                                                                                                                     tGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAAt               <  1:3545938/148‑1 (MQ=255)
                                                                                                                                       aCGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATgg             <  1:2126468/148‑1 (MQ=255)
                                                                                                                                          cgATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGttt          <  1:614080/148‑1 (MQ=255)
                                                                                                                                                aGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGt    >  2:1699027/1‑148 (MQ=255)
                                                                                                                                                 gTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGtt   <  1:2083589/148‑1 (MQ=255)
                                                                                                                                                  tGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTa  <  2:1697827/148‑1 (MQ=255)
                                                                                                                                                  tGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTa  >  2:1651998/1‑148 (MQ=255)
                                                                                                                                                  tGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCAACGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTa  >  2:3286308/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTA  >  CP009273/2877952‑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: