Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,636,207 T→A 9.3% 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.3% 969.6 / 70.2 281T96S (ACT→TCT) rsmJ16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (141/114);  new base A (14/12);  total (155/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 289 total aligned reads displayed.

TCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGC  >  CP009273/3636059‑3636355
                                                                                                                                                    |                                                                                                                                                    
tCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGt                                                                                                                                                      <  1:4053421/149‑1 (MQ=255)
        ccAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATcca                                                                                                                                              >  1:690947/1‑149 (MQ=255)
        ccAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATcca                                                                                                                                              <  1:3496907/149‑1 (MQ=255)
           gccgccGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATccacca                                                                                                                                           >  2:2822574/1‑149 (MQ=255)
            ccgccgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCcaccac                                                                                                                                          <  1:609775/149‑1 (MQ=255)
             cgccgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACcaca                                                                                                                                         >  1:2634242/1‑149 (MQ=255)
               ccgATTTCCGCATCCGCATAACCACGCGCTAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATc                                                                                                                                       <  2:3860812/149‑1 (MQ=255)
               ccgATTTCCGCATCCGCATAACCACGCGCTAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATc                                                                                                                                       <  2:2987656/149‑1 (MQ=255)
               ccgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATc                                                                                                                                       >  2:3960244/1‑149 (MQ=255)
                cgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATcc                                                                                                                                      <  2:3577525/149‑1 (MQ=255)
                cgATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATcc                                                                                                                                      >  2:3403870/1‑149 (MQ=255)
                          atccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAaataat                                                                                                                            <  2:152559/149‑1 (MQ=255)
                           tccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATc                                                                                                                           <  1:1025434/149‑1 (MQ=255)
                           tccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATc                                                                                                                           >  1:224645/1‑149 (MQ=255)
                            ccgcatAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCg                                                                                                                          <  1:1634224/149‑1 (MQ=255)
                                  aaCCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTa                                                                                                                    >  1:2017329/1‑149 (MQ=255)
                                   aCCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTaa                                                                                                                   <  2:718858/149‑1 (MQ=255)
                                   aCCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTaa                                                                                                                   >  1:436930/1‑149 (MQ=255)
                                    ccACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAAt                                                                                                                  <  1:3077207/149‑1 (MQ=255)
                                    ccACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAAt                                                                                                                  <  2:4124385/149‑1 (MQ=255)
                                      aCGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGc                                                                                                                >  2:307783/1‑149 (MQ=255)
                                         cgcCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCa                                                                                                             <  1:996392/149‑1 (MQ=255)
                                         cgcCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCa                                                                                                             <  1:1602220/149‑1 (MQ=255)
                                         cgcCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCa                                                                                                             <  1:3333301/149‑1 (MQ=255)
                                           ccAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCAcc                                                                                                           >  1:321026/1‑149 (MQ=255)
                                              ggCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCt                                                                                                        >  2:3947976/1‑149 (MQ=255)
                                              ggCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCt                                                                                                        >  2:3052379/1‑149 (MQ=255)
                                                ccGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCttt                                                                                                      <  2:2742333/149‑1 (MQ=255)
                                                ccGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCttt                                                                                                      >  2:187429/1‑149 (MQ=255)
                                                 cgtcgtCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTc                                                                                                     <  2:2613371/149‑1 (MQ=255)
                                                  gtcgtcGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCg                                                                                                    >  1:968481/1‑149 (MQ=255)
                                                  gtcgtcGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCg                                                                                                    >  1:4422066/1‑149 (MQ=255)
                                                     gtcgAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCa                                                                                                 >  1:4659663/1‑149 (MQ=255)
                                                      tcgAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCAc                                                                                                >  2:1891062/1‑149 (MQ=255)
                                                      tcgAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCAGCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCCCATAATCGCCTTTAATGCCCACCGCTTTCGCCAc                                                                                                >  1:1215077/1‑149 (MQ=255)
                                                             gcgcGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCa                                                                                         >  2:2946045/1‑149 (MQ=255)
                                                                 ggCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACcgc                                                                                     >  1:3967723/1‑149 (MQ=255)
                                                                  gCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACcgcg                                                                                    >  2:4391005/1‑149 (MQ=255)
                                                                  gCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACcgcg                                                                                    >  1:1295071/1‑149 (MQ=255)
                                                                  gCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACcgcg                                                                                    >  2:2182118/1‑149 (MQ=255)
                                                                    aaCCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGac                                                                                  >  2:587445/1‑149 (MQ=255)
                                                                       cACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGaccac                                                                               <  1:2185712/149‑1 (MQ=255)
                                                                          tGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCAccgc                                                                            >  2:1499460/1‑149 (MQ=255)
                                                                              ttACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGaa                                                                        <  2:4028344/149‑1 (MQ=255)
                                                                                aCGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAAtt                                                                      <  2:3552407/149‑1 (MQ=255)
                                                                                 cGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAAttt                                                                     <  2:4354513/149‑1 (MQ=255)
                                                                                   cTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGc                                                                   >  2:3335958/1‑149 (MQ=255)
                                                                                     ccAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGt                                                                 <  1:4555653/149‑1 (MQ=255)
                                                                                     ccAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGt                                                                 >  2:777879/1‑149 (MQ=255)
                                                                                     ccAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGt                                                                 <  1:3442451/149‑1 (MQ=255)
                                                                                     ccAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGt                                                                 >  2:160708/1‑149 (MQ=255)
                                                                                      cAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTc                                                                >  1:4264709/1‑149 (MQ=255)
                                                                                      cAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTc                                                                >  2:4569482/1‑149 (MQ=255)
                                                                                       aGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCg                                                               >  1:1185620/1‑149 (MQ=255)
                                                                                        gCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCgg                                                              <  1:906355/149‑1 (MQ=255)
                                                                                          aTCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTg                                                            >  1:2677719/1‑149 (MQ=255)
                                                                                              gCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGcc                                                        >  1:147272/1‑149 (MQ=255)
                                                                                              gCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGcc                                                        >  1:3576236/1‑149 (MQ=255)
                                                                                               cACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCa                                                       >  2:4053578/1‑149 (MQ=255)
                                                                                                 cgcgGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGAGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATc                                                     >  1:3158823/1‑149 (MQ=255)
                                                                                                     gCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTc                                                 <  2:4105012/149‑1 (MQ=255)
                                                                                                     gCAGCCGACTGAAGCCAGTACAAAGGAATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTc                                                 >  2:1274253/1‑149 (MQ=255)
                                                                                                       aGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCg                                               >  1:1599822/1‑149 (MQ=255)
                                                                                                        gCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGc                                              >  2:3695860/1‑149 (MQ=255)
                                                                                                        gCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGc                                              >  2:574022/1‑149 (MQ=255)
                                                                                                        gCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTGTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGc                                              >  2:574120/1‑149 (MQ=255)
                                                                                                          cGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCca                                            <  2:3671406/149‑1 (MQ=255)
                                                                                                          cGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCca                                            >  1:2947442/1‑149 (MQ=255)
                                                                                                          cGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCca                                            <  2:1335291/149‑1 (MQ=255)
                                                                                                           gACTGAAGCCAGTACAAAGGCCTCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGACTTTGCGTCGGTGCGCCATCGCTCCGCcac                                           >  1:1776294/1‑148 (MQ=255)
                                                                                                            aCTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCcaac                                          >  1:4295541/1‑149 (MQ=255)
                                                                                                               gAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACaaa                                       >  1:2133183/1‑149 (MQ=255)
                                                                                                               gAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACaaa                                       >  1:4482688/1‑149 (MQ=255)
                                                                                                                 aGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAAt                                     >  1:2878388/1‑149 (MQ=255)
                                                                                                                 aGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAAt                                     >  1:3883267/1‑149 (MQ=255)
                                                                                                                 aGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAAt                                     >  1:1523901/1‑149 (MQ=255)
                                                                                                                      gTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATcaaca                                <  2:237569/149‑1 (MQ=255)
                                                                                                                       tACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATcaacaa                               >  1:1112949/1‑149 (MQ=255)
                                                                                                                         cAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAg                             >  1:3633241/1‑149 (MQ=255)
                                                                                                                          aaaGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGa                            >  1:1587945/1‑149 (MQ=255)
                                                                                                                           aaGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAACTCAACAAAGAt                           >  2:3820855/1‑149 (MQ=255)
                                                                                                                           aaGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGAt                           >  2:3172207/1‑149 (MQ=255)
                                                                                                                           aaGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGAt                           >  2:214162/1‑149 (MQ=255)
                                                                                                                             ggCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgc                         <  1:663600/149‑1 (MQ=255)
                                                                                                                                 tCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgccgcc                     <  2:4242572/149‑1 (MQ=255)
                                                                                                                                 tCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgccgcc                     >  2:4546792/1‑149 (MQ=255)
                                                                                                                                 tCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgccgcc                     >  1:2925216/1‑149 (MQ=255)
                                                                                                                                   gcgcCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCaa                   <  1:636300/149‑1 (MQ=255)
                                                                                                                                   gcgcCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCaa                   <  1:2320352/149‑1 (MQ=255)
                                                                                                                                     gcCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGt                 <  2:690947/149‑1 (MQ=255)
                                                                                                                                         aaGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTg             <  2:3273958/149‑1 (MQ=255)
                                                                                                                                          aGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTgg            >  2:2485836/1‑149 (MQ=255)
                                                                                                                                                gCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAAGGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtcatc      <  2:887553/149‑1 (MQ=255)
                                                                                                                                                 cAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtcatca     <  2:501106/149‑1 (MQ=255)
                                                                                                                                                 cAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtcatca     <  2:1106435/149‑1 (MQ=255)
                                                                                                                                                  aGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCAc    <  1:1948709/149‑1 (MQ=255)
                                                                                                                                                   gTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACg   >  1:2848244/1‑149 (MQ=255)
                                                                                                                                                   gTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACg   >  1:4190706/1‑149 (MQ=255)
                                                                                                                                                    tGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGc  <  1:440300/149‑1 (MQ=255)
                                                                                                                                                    tGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGc  <  1:819844/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGC  >  CP009273/3636059‑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: