Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,636,207 T→A 9.8% T96S (ACT→TCT)  rsmJ ← 16S rRNA m(2)G1516 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,636,2070TA 9.8% 1224.3 / 101.4 357T96S (ACT→TCT) rsmJ16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (175/146);  new base A (23/12);  total (198/159)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 359 total aligned reads displayed.

CCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGC  >  CP009273/3636060‑3636355
                                                                                                                                                   |                                                                                                                                                    
ccTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTg                                                                                                                                                     >  2:804699/1‑149 (MQ=255)
 cTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTgg                                                                                                                                                    >  2:610285/1‑149 (MQ=255)
  tGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGc                                                                                                                                                   <  1:2290290/149‑1 (MQ=255)
   gCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCa                                                                                                                                                  <  2:3724372/149‑1 (MQ=255)
            cgccgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACcaca                                                                                                                                         >  2:1872713/1‑149 (MQ=255)
            cgccgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACcaca                                                                                                                                         >  2:4704718/1‑149 (MQ=255)
               cgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATcc                                                                                                                                      <  2:3291010/149‑1 (MQ=255)
               cgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATcc                                                                                                                                      <  1:4732463/149‑1 (MQ=255)
               cgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATcc                                                                                                                                      >  2:1220027/1‑149 (MQ=255)
               cgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATcc                                                                                                                                      <  2:2374728/149‑1 (MQ=255)
                gATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTCCAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCg                                                                                                                                     >  1:1836119/1‑149 (MQ=255)
                gATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCg                                                                                                                                     >  1:5117747/1‑149 (MQ=255)
                gATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCg                                                                                                                                     >  1:1604347/1‑149 (MQ=255)
                  tttCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGc                                                                                                                                   <  1:1030648/149‑1 (MQ=255)
                  tttCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGc                                                                                                                                   >  2:2097850/1‑149 (MQ=255)
                   ttCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCa                                                                                                                                  <  1:69790/149‑1 (MQ=255)
                     ccgcatccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCaaa                                                                                                                                <  2:3463338/149‑1 (MQ=255)
                     ccgcatccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCaaa                                                                                                                                <  2:837639/149‑1 (MQ=255)
                     ccgcatccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCaaa                                                                                                                                >  2:5018481/1‑149 (MQ=255)
                       gcatccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAaata                                                                                                                              <  2:2866943/149‑1 (MQ=255)
                         atccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAaataat                                                                                                                            >  2:5149755/1‑149 (MQ=255)
                         atccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAaataat                                                                                                                            >  1:3446902/1‑149 (MQ=255)
                          tccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATc                                                                                                                           >  2:2172017/1‑149 (MQ=255)
                             gcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGcc                                                                                                                        >  2:5143322/1‑149 (MQ=255)
                              catAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCACATAATCGCCt                                                                                                                       >  2:3514333/1‑149 (MQ=255)
                              catAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCt                                                                                                                       <  2:146541/149‑1 (MQ=255)
                                tAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCttt                                                                                                                     >  2:2664214/1‑149 (MQ=255)
                                tAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCttt                                                                                                                     >  2:2003501/1‑149 (MQ=255)
                                tAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCttt                                                                                                                     >  1:2675635/1‑149 (MQ=255)
                                 aaCCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTa                                                                                                                    <  2:973034/149‑1 (MQ=255)
                                 aaCCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTa                                                                                                                    <  1:1810806/149‑1 (MQ=255)
                                 aaCCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTa                                                                                                                    <  1:1146249/149‑1 (MQ=255)
                                   ccACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAAt                                                                                                                  >  1:3761480/1‑149 (MQ=255)
                                     aCGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGc                                                                                                                >  2:177386/1‑149 (MQ=255)
                                     aCGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGc                                                                                                                >  1:3553660/1‑149 (MQ=255)
                                          ccAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCAcc                                                                                                           <  1:957620/149‑1 (MQ=255)
                                            aGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGc                                                                                                         >  1:850462/1‑149 (MQ=255)
                                             ggCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCt                                                                                                        >  2:4460048/1‑149 (MQ=255)
                                              gCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCtt                                                                                                       <  2:4919005/149‑1 (MQ=255)
                                               ccGTCGTCGAGTAGCGCGTCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCttt                                                                                                      <  2:4719050/149‑1 (MQ=255)
                                               ccGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCttt                                                                                                      <  1:2431397/149‑1 (MQ=255)
                                                cgtcgtCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTc                                                                                                     <  1:1639962/149‑1 (MQ=255)
                                                cgtcgtCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTc                                                                                                     <  1:1899926/149‑1 (MQ=255)
                                                cgtcgtCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTc                                                                                                     <  2:2484982/149‑1 (MQ=255)
                                                  tcgtcgAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGc                                                                                                   >  1:4262467/1‑149 (MQ=255)
                                                          tAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCt                                                                                           <  2:1208410/149‑1 (MQ=255)
                                                              gcgGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCAcc                                                                                       >  2:2051733/1‑149 (MQ=255)
                                                                 gCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACcgcg                                                                                    >  2:3196595/1‑149 (MQ=255)
                                                                       aCTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGaccacc                                                                              >  1:1498690/1‑149 (MQ=255)
                                                                         tGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCAccgc                                                                            >  2:3415322/1‑149 (MQ=255)
                                                                         tGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCAccgc                                                                            >  1:2630881/1‑149 (MQ=255)
                                                                           gATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCAccgccg                                                                          >  1:4341705/1‑149 (MQ=255)
                                                                           gATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCAccgccg                                                                          >  1:2993845/1‑149 (MQ=255)
                                                                             ttACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGaa                                                                        >  2:1165909/1‑149 (MQ=255)
                                                                             ttACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGaa                                                                        >  1:782878/1‑149 (MQ=255)
                                                                             ttACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGaa                                                                        >  1:2583939/1‑149 (MQ=255)
                                                                              tACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAAt                                                                       >  2:1495581/1‑149 (MQ=255)
                                                                                 gCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTg                                                                    >  2:5122677/1‑149 (MQ=255)
                                                                                   tCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCg                                                                  >  2:1138435/1‑149 (MQ=255)
                                                                                     cAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTc                                                                <  1:3645682/149‑1 (MQ=255)
                                                                                            cGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTgcgc                                                         <  1:382889/149‑1 (MQ=255)
                                                                                             gCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGcc                                                        >  1:4077279/1‑149 (MQ=255)
                                                                                               aCGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCAt                                                      <  1:4634452/149‑1 (MQ=255)
                                                                                                cgcgGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATc                                                     <  2:1785707/149‑1 (MQ=255)
                                                                                                 gcgGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCg                                                    >  1:3086622/1‑149 (MQ=255)
                                                                                                    gCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCCCCACATCCGGCACAGAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTc                                                 <  2:3922596/149‑1 (MQ=255)
                                                                                                    gCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTc                                                 <  1:3808384/149‑1 (MQ=255)
                                                                                                       gCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCCCCGc                                              <  2:1404316/149‑1 (MQ=255)
                                                                                                         cGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCca                                            <  2:4519524/149‑1 (MQ=255)
                                                                                                          gACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCcaa                                           <  2:2326031/149‑1 (MQ=255)
                                                                                                          gACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCcaa                                           >  1:1989239/1‑149 (MQ=255)
                                                                                                             tGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCcaacaa                                        <  2:4070512/149‑1 (MQ=255)
                                                                                                                aGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAAt                                     <  2:1911707/149‑1 (MQ=255)
                                                                                                                   cAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATcaa                                  <  1:1267057/149‑1 (MQ=255)
                                                                                                                       aCAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACaaa                              >  1:2195618/1‑149 (MQ=255)
                                                                                                                       aCAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACaaa                              >  1:900754/1‑149 (MQ=255)
                                                                                                                        cAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAg                             <  2:1687457/149‑1 (MQ=255)
                                                                                                                        cAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAg                             >  2:3096495/1‑149 (MQ=255)
                                                                                                                         aaaGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACaaaa                                      >  2:2245122/1‑139 (MQ=255)
                                                                                                                         aaaGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGa                            <  1:3269921/149‑1 (MQ=255)
                                                                                                                            ggCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgc                         >  1:1347714/1‑149 (MQ=255)
                                                                                                                            ggCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgc                         >  2:3653402/1‑149 (MQ=255)
                                                                                                                               aTCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgccgc                      >  1:2829921/1‑149 (MQ=255)
                                                                                                                               aTCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgccgc                      <  1:5149755/149‑1 (MQ=255)
                                                                                                                               aTCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgccgc                      >  1:1953795/1‑149 (MQ=255)
                                                                                                                                   cgcCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAg                  >  2:2932010/1‑149 (MQ=255)
                                                                                                                                   cgcCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAg                  <  1:3201205/149‑1 (MQ=255)
                                                                                                                                     cccAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGtt                >  2:1027124/1‑149 (MQ=255)
                                                                                                                                      ccAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGttt               >  2:2149479/1‑149 (MQ=255)
                                                                                                                                       cAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGtttt              <  2:4220332/149‑1 (MQ=255)
                                                                                                                                       cAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGtttt              <  1:4930137/149‑1 (MQ=255)
                                                                                                                                       cAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGtttt              >  2:3755502/1‑149 (MQ=255)
                                                                                                                                          gTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGc           <  2:1025786/149‑1 (MQ=255)
                                                                                                                                             cTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtca        <  2:4019884/149‑1 (MQ=255)
                                                                                                                                              tGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtcat       >  1:2387475/1‑149 (MQ=255)
                                                                                                                                               gCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtcatc      >  2:650885/1‑149 (MQ=255)
                                                                                                                                                 aGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCAc    <  2:1574397/149‑1 (MQ=255)
                                                                                                                                                  gTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACg   >  2:4749116/1‑149 (MQ=255)
                                                                                                                                                  gTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACg   >  1:4867766/1‑149 (MQ=255)
                                                                                                                                                   tGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGc  <  1:2797795/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGC  >  CP009273/3636060‑3636355

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: ‑