Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,879,222 C→A 6.6% E453* (GAA→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,879,2220CA 6.6% 1113.8 / 32.2 301E453* (GAA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base C (154/126);  new base A (13/8);  total (167/134)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

TATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAA  >  CP009273/2879075‑2879370
                                                                                                                                                   |                                                                                                                                                    
     aGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAAAAACACTCCCTCCCACATCAGCCTGAGCCGTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATcaac                                                                                                                                                >  1:2792572/1‑149 (MQ=255)
      gACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATcaaca                                                                                                                                               >  1:3554838/1‑149 (MQ=255)
       aCCATAAGTATCCAGAAGTTTCTGTTTTTTTTTCATTGGTAGGGTTGCGGGAAGAAGAATAACACTCCCTCCCACAACAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGGTCATATAGGTGTCGTAAGCATGAACATAATcaacaa                                                                                                                                              >  1:4123407/1‑149 (MQ=255)
       aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTAATcaacaa                                                                                                                                              >  2:4696210/1‑149 (MQ=255)
        ccATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTAATCAACAAt                                                                                                                                             <  2:312743/149‑1 (MQ=255)
         cATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAAtt                                                                                                                                            <  2:464477/149‑1 (MQ=255)
          aTAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTTAGCACTTCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTa                                                                                                                                           >  2:438135/1‑149 (MQ=255)
          aTAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTa                                                                                                                                           >  2:4312028/1‑149 (MQ=255)
             aGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTaaaa                                                                                                                                        <  2:906767/149‑1 (MQ=255)
             aGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTaaaa                                                                                                                                        <  1:2855370/149‑1 (MQ=255)
              gTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAac                                                                                                                                       <  1:3287557/149‑1 (MQ=255)
                   cAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTc                                                                                                                                  >  2:3891290/1‑149 (MQ=255)
                    aGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTACCTCCCACATCAGCCTGAGCCTTGAGCACAGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCAAGAAATTCATCAACAATTAAAACACTTCg                                                                                                                                 >  2:3038379/1‑149 (MQ=255)
                          ttCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAAt                                                                                                                           >  2:1885146/1‑149 (MQ=255)
                          ttCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAAt                                                                                                                           <  1:2041096/149‑1 (MQ=255)
                           tCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAAtt                                                                                                                          <  1:1254108/149‑1 (MQ=255)
                              gTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTccc                                                                                                                       >  2:3385444/1‑149 (MQ=255)
                               tttttGTTTCATTGGTAGGGTTGCGGAAAGAAGAGTAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCa                                                                                                                      >  1:1361900/1‑149 (MQ=255)
                               tttttGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCa                                                                                                                      >  1:3422910/1‑149 (MQ=255)
                                ttttGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCaa                                                                                                                     >  1:3700375/1‑149 (MQ=255)
                                 tttGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCaaa                                                                                                                    <  2:119019/149‑1 (MQ=255)
                                    gTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCa                                                                                                                 <  1:1883289/149‑1 (MQ=255)
                                         aTTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGAt                                                                                                            <  1:1532930/149‑1 (MQ=255)
                                          ttGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATTAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATa                                                                                                           >  2:501661/1‑149 (MQ=255)
                                          ttGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATa                                                                                                           >  2:4279224/1‑149 (MQ=255)
                                           tGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATaa                                                                                                          <  2:2243388/149‑1 (MQ=255)
                                            ggTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATaaa                                                                                                         <  2:3971442/149‑1 (MQ=255)
                                             gTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTAGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAg                                                                                                        >  1:4215818/1‑149 (MQ=255)
                                              tAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGc                                                                                                       <  2:278134/149‑1 (MQ=255)
                                               aGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCg                                                                                                      >  2:142167/1‑149 (MQ=255)
                                                 ggTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGgt                                                                                                    <  1:929749/149‑1 (MQ=255)
                                                    tGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGtt                                                                                                 <  2:4371714/149‑1 (MQ=255)
                                                      cGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTa                                                                                               <  2:102794/149‑1 (MQ=255)
                                                         aaaGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACt                                                                                            <  1:2458335/149‑1 (MQ=255)
                                                         aaaGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACt                                                                                            <  1:3614876/149‑1 (MQ=255)
                                                         aaaGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCAGCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACt                                                                                            >  2:1308888/1‑149 (MQ=255)
                                                           agaagaATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTgg                                                                                          >  2:72106/1‑149 (MQ=255)
                                                              agaaTAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCaa                                                                                       >  2:4604537/1‑149 (MQ=255)
                                                               gaaTAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAAt                                                                                      >  2:1255693/1‑149 (MQ=255)
                                                                aaTAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATa                                                                                     <  1:967844/149‑1 (MQ=255)
                                                                aaTAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATa                                                                                     <  1:1571650/149‑1 (MQ=255)
                                                                      acTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGata                                                                               <  1:1140313/149‑1 (MQ=255)
                                                                      acTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGata                                                                               <  1:4547105/149‑1 (MQ=255)
                                                                         ccctcccACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTAATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCa                                                                            >  1:3499135/1‑149 (MQ=255)
                                                                            tcccACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTAATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAaca                                                                         >  1:4649843/1‑149 (MQ=255)
                                                                             cccACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAacac                                                                        >  2:2930605/1‑149 (MQ=255)
                                                                              ccACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTAATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACAcc                                                                       <  1:4468599/149‑1 (MQ=255)
                                                                                acaTCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTg                                                                     <  1:3861385/149‑1 (MQ=255)
                                                                                 caTCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGa                                                                    <  1:2957308/149‑1 (MQ=255)
                                                                                   tCAGCCGGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACACTTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATc                                                                  <  2:1142211/149‑1 (MQ=255)
                                                                                    cAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCa                                                                 <  2:4478307/149‑1 (MQ=255)
                                                                                     aGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCACACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCaa                                                                <  2:2016664/149‑1 (MQ=255)
                                                                                     aGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCaa                                                                >  1:1971706/1‑149 (MQ=255)
                                                                                      gCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAAt                                                               >  1:613945/1‑149 (MQ=255)
                                                                                      gCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAAt                                                               <  2:656100/149‑1 (MQ=255)
                                                                                      gCCTGAGCCATGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAAt                                                               <  2:2168052/149‑1 (MQ=255)
                                                                                         tGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGt                                                            >  1:4329761/1‑149 (MQ=255)
                                                                                          gAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTg                                                           >  1:2738908/1‑149 (MQ=255)
                                                                                            gacTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAAAACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATAACGATATCAACACCTGATCAATAGTGCa                                                         >  1:1842584/3‑149 (MQ=255)
                                                                                              cTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCaaa                                                       >  2:3264669/1‑149 (MQ=255)
                                                                                              cTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCaaa                                                       <  1:4021831/149‑1 (MQ=255)
                                                                                              cTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCACAAACCACGGATAAAGCGGTGTTTAACTGGCAATAACGATATCAACACCTGATCAATCGTGCaaa                                                       >  2:980122/1‑149 (MQ=255)
                                                                                               ttGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAAc                                                      <  1:1769855/149‑1 (MQ=255)
                                                                                                 gAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGc                                                    <  1:845336/149‑1 (MQ=255)
                                                                                                       tGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGAttt                                              <  1:472953/149‑1 (MQ=255)
                                                                                                       tGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGAttt                                              >  1:68750/1‑149 (MQ=255)
                                                                                                        gCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTg                                             >  2:4275236/1‑149 (MQ=255)
                                                                                                          cTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGcc                                           >  2:2034336/1‑149 (MQ=255)
                                                                                                           tCCATCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTACCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGAGCAAACGCAGATATGacc                                          >  1:164172/1‑149 (MQ=255)
                                                                                                           tCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGccc                                          >  1:2078593/1‑149 (MQ=255)
                                                                                                             cAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCaa                                        >  1:3414344/1‑149 (MQ=255)
                                                                                                             cAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTAATCAACAATTAAAACACTTCGACCACTTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCaa                                        <  2:1583478/149‑1 (MQ=255)
                                                                                                              aGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAg                                       >  2:3354233/1‑149 (MQ=255)
                                                                                                               gCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGa                                      >  1:1826152/1‑149 (MQ=255)
                                                                                                                cAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGaa                                     >  2:2731780/1‑149 (MQ=255)
                                                                                                                   gCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAaca                                  >  2:1114635/1‑149 (MQ=255)
                                                                                                                    ccGTTCATATAGGTGTCGTAAGCACGAACTACATCAACAATTAAAACACTACGACCAATTCCCAAACCACGGATAAAGCGTTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAacac                                 <  2:1930303/149‑1 (MQ=255)
                                                                                                                      gTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACtt                               >  2:2344554/1‑149 (MQ=255)
                                                                                                                       ttCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACttt                              >  2:207764/1‑149 (MQ=255)
                                                                                                                         cATATAGGTGTCGTAAGCATGAACTTCACCAACAACTAAAACACTTCGACCACTTCCCAACCCACGCATAAAGCGGTGTTTAACGGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCt                            <  2:1728906/149‑1 (MQ=255)
                                                                                                                           tataGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTa                          >  1:3355330/1‑149 (MQ=255)
                                                                                                                            ataGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTAt                         <  2:1763828/149‑1 (MQ=255)
                                                                                                                             taGGTGTCGTAATCATGAACTTCATAAACAATTAAAACACTTCGACCAATTCACAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTAtt                        >  1:4781684/1‑149 (MQ=255)
                                                                                                                             taGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTAtt                        >  1:4045555/1‑149 (MQ=255)
                                                                                                                             taGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTAtt                        >  2:1521052/1‑149 (MQ=255)
                                                                                                                             taGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATACCCACACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTAtt                        <  2:3641277/149‑1 (MQ=255)
                                                                                                                             taGGTGTCGTAAGCATGAACTTAATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGTATACAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTAtt                        <  1:2207245/149‑1 (MQ=255)
                                                                                                                                gtgtCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAAGACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCt                     <  2:3397458/149‑1 (MQ=255)
                                                                                                                                   tCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTtg                  <  2:3562007/149‑1 (MQ=255)
                                                                                                                                   tCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTtg                  >  1:3259110/1‑149 (MQ=255)
                                                                                                                                     gTAAGCATGAACTTAATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTtgtg                >  1:2429378/1‑149 (MQ=255)
                                                                                                                                      tAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGa               >  2:731394/1‑149 (MQ=255)
                                                                                                                                       aaGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGAc              >  1:1500524/1‑149 (MQ=255)
                                                                                                                                         gCATGAACTTAATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACaa            >  1:3790133/1‑149 (MQ=255)
                                                                                                                                           aTGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAAcc          >  1:4154154/1‑149 (MQ=255)
                                                                                                                                            tGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCa         >  2:420739/1‑149 (MQ=255)
                                                                                                                                             gAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCAc        <  1:2338471/149‑1 (MQ=255)
                                                                                                                                               aCTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTg      >  1:2481370/1‑149 (MQ=255)
                                                                                                                                                 ttCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGac    >  2:3834971/1‑149 (MQ=255)
                                                                                                                                                   catcaACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGacaa  <  2:1826500/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAA  >  CP009273/2879075‑2879370

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: