Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,880,202 G→T 5.7% S126* (TCA→TAA)  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,880,2020GT 5.7% 1126.9 / 28.2 297S126* (TCA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (139/141);  new base T (6/11);  total (145/152)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 302 total aligned reads displayed.

TTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTCCA  >  CP009273/2880053‑2880351
                                                                                                                                                     |                                                                                                                                                     
 ttCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGa                                                                                                                                                      >  2:916246/1‑150 (MQ=255)
  tCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAACCCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGaa                                                                                                                                                     <  2:4102130/150‑1 (MQ=255)
    ccAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATc                                                                                                                                                   <  2:136905/150‑1 (MQ=255)
    ccAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATc                                                                                                                                                   >  1:2160149/1‑150 (MQ=255)
      aaCGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCt                                                                                                                                                 >  2:2998052/1‑150 (MQ=255)
        cgcgACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTgg                                                                                                                                               >  2:3576723/1‑150 (MQ=255)
         gcgACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGt                                                                                                                                              <  1:3722523/150‑1 (MQ=255)
           gACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTa                                                                                                                                            >  1:443221/1‑150 (MQ=255)
           gACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTa                                                                                                                                            >  1:2085789/1‑150 (MQ=255)
              ttATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAAc                                                                                                                                         <  2:4107067/150‑1 (MQ=255)
                aTCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCa                                                                                                                                       <  1:3068049/150‑1 (MQ=255)
                 tCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCaa                                                                                                                                      >  1:1223109/1‑150 (MQ=255)
                     gATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACa                                                                                                                                  <  2:1148966/150‑1 (MQ=255)
                       tCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGa                                                                                                                                <  1:4112947/150‑1 (MQ=255)
                         cTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGAcc                                                                                                                              <  2:1963242/150‑1 (MQ=255)
                         cTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGAcc                                                                                                                              >  2:681170/1‑150 (MQ=255)
                          tGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCg                                                                                                                             >  1:2864906/1‑150 (MQ=255)
                          tGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCg                                                                                                                             >  1:965806/1‑150 (MQ=255)
                          tGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATACTGGTTAAACCAATACAGACCg                                                                                                                             <  1:1917541/150‑1 (MQ=255)
                          tGGGAATGTAATATAAAACCATGATGTCCTGTAACGACCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCg                                                                                                                             <  2:1923603/150‑1 (MQ=255)
                           gggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCgg                                                                                                                            >  1:7152/1‑150 (MQ=255)
                            ggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAATCCAATACAGACCGGc                                                                                                                           >  2:3944248/1‑150 (MQ=255)
                            ggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGc                                                                                                                           >  1:2040640/1‑150 (MQ=255)
                            ggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGc                                                                                                                           <  1:4156921/150‑1 (MQ=255)
                            ggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGc                                                                                                                           >  2:1688293/1‑150 (MQ=255)
                            ggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGc                                                                                                                           >  2:1702125/1‑150 (MQ=255)
                            ggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGc                                                                                                                           >  2:3961075/1‑150 (MQ=255)
                              aaTGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTg                                                                                                                         >  1:3946589/1‑150 (MQ=255)
                               aTGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGc                                                                                                                        <  1:105414/150‑1 (MQ=255)
                               aTGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGc                                                                                                                        >  1:3405116/1‑150 (MQ=255)
                                   aaTATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATTCTGGTTAAACAAATACAGACCGGCTGCACCa                                                                                                                    <  2:2079114/150‑1 (MQ=255)
                                   aaTATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCa                                                                                                                    >  1:1087143/1‑150 (MQ=255)
                                    atatAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCAt                                                                                                                   >  2:2276702/1‑150 (MQ=255)
                                       taAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGAt                                                                                                                <  1:1612601/150‑1 (MQ=255)
                                        aaaaCCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGAtt                                                                                                               <  1:2836805/150‑1 (MQ=255)
                                         aaaCCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTa                                                                                                              >  2:3422402/1‑150 (MQ=255)
                                              atgatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAttt                                                                                                         >  2:3784907/1‑150 (MQ=255)
                                              atgatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAttt                                                                                                         <  1:4006387/150‑1 (MQ=255)
                                               tgatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTa                                                                                                        >  2:2724149/1‑150 (MQ=255)
                                                    tCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGca                                                                                                   >  1:545741/1‑150 (MQ=255)
                                                     ccTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGcac                                                                                                  >  2:103874/1‑150 (MQ=255)
                                                       tGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAt                                                                                                <  2:3477633/150‑1 (MQ=255)
                                                        gTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAtt                                                                                               <  2:2809266/150‑1 (MQ=255)
                                                         tAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAttt                                                                                              <  2:2588447/150‑1 (MQ=255)
                                                         tAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAttt                                                                                              <  2:1625217/150‑1 (MQ=255)
                                                          aaCGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTtg                                                                                             <  2:2929392/150‑1 (MQ=255)
                                                               ccTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTg                                                                                        >  2:2148933/1‑150 (MQ=255)
                                                                    aCCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGa                                                                                   >  2:3132147/1‑150 (MQ=255)
                                                                    aCCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGa                                                                                   >  1:648838/1‑150 (MQ=255)
                                                                     cccATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGAc                                                                                  <  2:815303/150‑1 (MQ=255)
                                                                         tGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAACAACTGAAAAAATCCCCGTGAGTTTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATGTACGGCTCATTTGTGTTGATGGACCAtt                                                                              <  2:499943/150‑1 (MQ=255)
                                                                          ggAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTa                                                                             <  2:1419015/150‑1 (MQ=255)
                                                                          ggAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTa                                                                             <  1:1183715/150‑1 (MQ=255)
                                                                             aaCCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGt                                                                          >  2:875992/1‑150 (MQ=255)
                                                                                 aGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATg                                                                      <  2:3493365/150‑1 (MQ=255)
                                                                                 aGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATg                                                                      <  2:4437153/150‑1 (MQ=255)
                                                                                 aGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATg                                                                      <  2:1233345/150‑1 (MQ=255)
                                                                                 aGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATg                                                                      <  2:3493238/150‑1 (MQ=255)
                                                                                  ggACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATgg                                                                     <  1:257008/150‑1 (MQ=255)
                                                                                    aCTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCg                                                                   >  2:1599770/1‑150 (MQ=255)
                                                                                     cTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGt                                                                  >  1:1374884/1‑150 (MQ=255)
                                                                                      tCATAAGGATGAGGAGCGGCATCAAAAAAACTGACAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCggt                                                                 <  2:3697655/150‑3 (MQ=255)
                                                                                       cATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTg                                                                <  1:1859607/150‑1 (MQ=255)
                                                                                        aTAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGc                                                               <  2:2702636/150‑1 (MQ=255)
                                                                                          aaGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAg                                                             <  2:143757/150‑1 (MQ=255)
                                                                                           aGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAgg                                                            >  1:3285961/1‑150 (MQ=255)
                                                                                              aTGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGAtt                                                         >  2:4325651/1‑150 (MQ=255)
                                                                                                    agCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAAttt                                                   <  2:3988385/150‑1 (MQ=255)
                                                                                                      cGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCa                                                 <  1:2460491/150‑1 (MQ=255)
                                                                                                       ggCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAg                                                >  2:1764182/1‑150 (MQ=255)
                                                                                                            cAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAAc                                           >  1:1697234/1‑150 (MQ=255)
                                                                                                            cAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGGTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAAc                                           >  1:1309604/1‑150 (MQ=255)
                                                                                                             aaaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACt                                          <  1:2724153/150‑1 (MQ=255)
                                                                                                             aaaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACt                                          >  1:2215009/1‑150 (MQ=255)
                                                                                                              aaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAGTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACtg                                         >  1:2639525/1‑149 (MQ=255)
                                                                                                               aaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACttt                                        <  2:3820982/150‑1 (MQ=255)
                                                                                                               aaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACttt                                        <  2:568918/150‑1 (MQ=255)
                                                                                                                  aaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctg                                     <  1:2713538/150‑1 (MQ=255)
                                                                                                                  aaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctg                                     >  1:3185167/1‑150 (MQ=255)
                                                                                                                    cTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgct                                   <  1:1006911/150‑1 (MQ=255)
                                                                                                                     tGAACAAATCCCCGAGAGAGTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgctg                                  <  2:3309173/150‑1 (MQ=255)
                                                                                                                     tGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgctg                                  >  1:1095212/1‑150 (MQ=255)
                                                                                                                      gAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGa                                 >  1:461879/1‑150 (MQ=255)
                                                                                                                       aaaaaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGAt                                >  1:4309836/1‑150 (MQ=255)
                                                                                                                       aaaaaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGAt                                >  1:1717061/1‑150 (MQ=255)
                                                                                                                       aaaaaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAGTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGAt                                >  2:3254194/1‑150 (MQ=255)
                                                                                                                        aaaaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGAtt                               >  1:796359/1‑150 (MQ=255)
                                                                                                                           aaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTtat                            >  1:1437295/1‑150 (MQ=255)
                                                                                                                            aTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTtata                           <  2:3353043/150‑1 (MQ=255)
                                                                                                                             tCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCACTACAGACCGGCTGCACCATGATTAAATTTACGGCACAATTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTatat                          <  2:3731483/150‑1 (MQ=255)
                                                                                                                                  gagagaTTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGaa                     >  1:3250453/1‑150 (MQ=255)
                                                                                                                                   agagaTTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAAt                    <  2:2007174/150‑1 (MQ=255)
                                                                                                                                   agagaTTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAAt                    >  2:63062/1‑150 (MQ=255)
                                                                                                                                   agagaTTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAAt                    <  2:2448758/150‑1 (MQ=255)
                                                                                                                                    gagaTTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATGTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATc                   >  2:2333082/1‑150 (MQ=255)
                                                                                                                                        ttGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtat               <  1:875992/150‑1 (MQ=255)
                                                                                                                                        ttGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtat               <  2:731651/150‑1 (MQ=255)
                                                                                                                                          gCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTatat             <  2:465195/150‑1 (MQ=255)
                                                                                                                                           ctctGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATc            <  1:63062/150‑1 (MQ=255)
                                                                                                                                                   gTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTc    >  1:1352736/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTCCA  >  CP009273/2880053‑2880351

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: