Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,056,703 C→A 6.5% R120S (CGC→AGC)  argK → membrane ATPase/protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,056,7030CA 6.5% 1217.9 / 39.8 331R120S (CGC→AGC) argKmembrane ATPase/protein kinase
Reads supporting (aligned to +/- strand):  ref base C (146/160);  new base A (18/7);  total (164/167)
Fisher's exact test for biased strand distribution p-value = 2.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 335 total aligned reads displayed.

TTGAGGCCTTTGGCATGTTGTTGATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGA  >  CP009273/3056557‑3056851
                                                                                                                                                  |                                                                                                                                                    
ttGAGGCCTTTGGCATGTTGTTGATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGc                                                                                                                                                    >  2:1405543/1‑149 (MQ=255)
 tGAGGCCTTTGGCATGTTGTTGATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGGTTATTCGcc                                                                                                                                                   >  1:4156163/1‑149 (MQ=255)
    ggCCTTTGGCATGTTGTTGATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCgg                                                                                                                                                >  2:4035320/1‑149 (MQ=255)
      ccTTTGGCATGTTGTTGATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTa                                                                                                                                              >  1:1308541/1‑149 (MQ=255)
       cTTTGGCATGTTGTTGATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTAc                                                                                                                                             <  1:4864533/149‑1 (MQ=255)
            gCATGTTGTTGATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCAtcc                                                                                                                                        <  1:3857282/149‑1 (MQ=255)
                    ttgATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCa                                                                                                                                >  1:4166258/1‑149 (MQ=255)
                    ttgATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCACCGGTCa                                                                                                                                >  1:3585677/1‑149 (MQ=255)
                       aTTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCACCGGTCATCt                                                                                                                             >  1:1052364/1‑149 (MQ=255)
                        ttCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCTCCGGTCATCTg                                                                                                                            >  2:5009745/1‑149 (MQ=255)
                         tCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTgg                                                                                                                           <  1:2389828/149‑1 (MQ=255)
                            agagGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCg                                                                                                                        >  1:4487625/1‑149 (MQ=255)
                            agagGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCg                                                                                                                        >  1:2782082/1‑149 (MQ=255)
                                ggATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGc                                                                                                                    <  2:5092860/149‑1 (MQ=255)
                                 gATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGcc                                                                                                                   >  2:2778117/1‑149 (MQ=255)
                                 gATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCAACCTCCGGTCATCTGGGCGGTGcc                                                                                                                   >  1:4257150/1‑149 (MQ=255)
                                    tAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGt                                                                                                                <  2:4901037/149‑1 (MQ=255)
                                      aaGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCa                                                                                                              <  1:2428276/149‑1 (MQ=255)
                                      aaGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCa                                                                                                              >  2:453258/1‑149 (MQ=255)
                                       aGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCag                                                                                                             >  2:1956963/1‑149 (MQ=255)
                                        ggTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCagc                                                                                                            <  1:373363/149‑1 (MQ=255)
                                         gTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCagcg                                                                                                           >  2:3485167/1‑149 (MQ=255)
                                         gTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCCGTACCAACCACCGGTCATCTGGGCGGTGCCAGTCagcg                                                                                                           >  2:887723/1‑149 (MQ=255)
                                          tCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCagcga                                                                                                          >  2:745659/1‑149 (MQ=255)
                                               gTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGcgcg                                                                                                     >  2:3240922/1‑149 (MQ=255)
                                                 tATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCggg                                                                                                   <  2:5042024/149‑1 (MQ=255)
                                                   ttGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGaa                                                                                                 <  2:2132983/149‑1 (MQ=255)
                                                    tGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCTTGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAAt                                                                                                <  1:4654129/149‑1 (MQ=255)
                                                       ggTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTaa                                                                                             >  1:1664946/1‑149 (MQ=255)
                                                        gTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCTCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAaataat                                                                                            >  1:2529874/1‑149 (MQ=255)
                                                        gTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAAt                                                                                            >  2:2389499/1‑149 (MQ=255)
                                                        gTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAAt                                                                                            <  2:2286721/149‑1 (MQ=255)
                                                        gTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAAt                                                                                            >  2:3820456/1‑149 (MQ=255)
                                                        gTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAAt                                                                                            >  2:3820425/1‑149 (MQ=255)
                                                         tCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATg                                                                                           >  2:4216908/1‑149 (MQ=255)
                                                          cGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGc                                                                                          <  2:551447/149‑1 (MQ=255)
                                                           gATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCt                                                                                         <  2:4768075/149‑1 (MQ=255)
                                                            aTCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTg                                                                                        >  1:4579198/1‑149 (MQ=255)
                                                            aTCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTg                                                                                        <  2:3356508/149‑1 (MQ=255)
                                                            aTCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTg                                                                                        >  1:1546012/1‑149 (MQ=255)
                                                            aTCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCTCCGGTCAACTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTg                                                                                        >  1:4819964/1‑149 (MQ=255)
                                                             tCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGt                                                                                       <  1:2613005/149‑1 (MQ=255)
                                                              cccAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGtt                                                                                      <  1:2955777/149‑1 (MQ=255)
                                                               ccAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTa                                                                                     >  1:3139099/1‑149 (MQ=255)
                                                                     gcCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGCGCGCGGGAATTAATGCTGTTATGCGaa                                                                               >  1:776705/1‑149 (MQ=255)
                                                                     gcCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGaa                                                                               <  2:265823/149‑1 (MQ=255)
                                                                        cGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAagca                                                                            <  2:1441824/149‑1 (MQ=255)
                                                                         ggTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAagcag                                                                           <  1:5195653/149‑1 (MQ=255)
                                                                         ggTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAagcag                                                                           <  2:4647029/149‑1 (MQ=255)
                                                                          gTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAagcagc                                                                          <  1:806827/149‑1 (MQ=255)
                                                                          gTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAagcagc                                                                          <  1:4773795/149‑1 (MQ=255)
                                                                           tCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCg                                                                         <  1:2376242/149‑1 (MQ=255)
                                                                              cTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGt                                                                      <  1:1710301/149‑1 (MQ=255)
                                                                              cTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGt                                                                      >  2:3580075/1‑149 (MQ=255)
                                                                               tGGCGGTAGCAGTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGtt                                                                     <  2:5160189/149‑1 (MQ=255)
                                                                                ggCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTa                                                                    >  1:1259713/1‑149 (MQ=255)
                                                                                 gCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTAt                                                                   <  2:1781361/149‑1 (MQ=255)
                                                                                   ggTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGa                                                                 <  2:40608/149‑1 (MQ=255)
                                                                                    gTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGAc                                                                <  2:1218103/149‑1 (MQ=255)
                                                                                       gCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTa                                                             <  1:2003019/149‑1 (MQ=255)
                                                                                       gCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTa                                                             >  2:3096131/1‑149 (MQ=255)
                                                                                          ttCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTTATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTg                                                          <  2:4849890/149‑1 (MQ=255)
                                                                                           tctcGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGa                                                         <  1:5020979/149‑1 (MQ=255)
                                                                                             tcGGGGATAACACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGAtt                                                       <  2:1585158/149‑1 (MQ=255)
                                                                                             tcGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGAtt                                                       <  2:1935026/149‑1 (MQ=255)
                                                                                             tcGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGAtt                                                       <  2:4732919/149‑1 (MQ=255)
                                                                                               ggggATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGt                                                     <  2:438898/149‑1 (MQ=255)
                                                                                                    tAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGaaa                                                <  2:3719183/149‑1 (MQ=255)
                                                                                                    tAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGaaa                                                <  2:1123004/149‑1 (MQ=255)
                                                                                                     aaaaCCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAAc                                               <  1:145325/149‑1 (MQ=255)
                                                                                                     aaaaCCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAAc                                               <  2:3801408/149‑1 (MQ=255)
                                                                                                       aaCCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACgg                                             <  2:4943719/149‑1 (MQ=255)
                                                                                                       aaCCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACgg                                             >  2:2726009/1‑149 (MQ=255)
                                                                                                        aCCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGt                                            >  2:2498772/1‑149 (MQ=255)
                                                                                                        aCCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGt                                            <  1:745659/149‑1 (MQ=255)
                                                                                                        aCCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGt                                            >  1:1814668/1‑149 (MQ=255)
                                                                                                          ccacATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTg                                          <  2:365681/146‑1 (MQ=255)
                                                                                                             cATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCg                                       >  1:3300088/1‑149 (MQ=255)
                                                                                                               tgaatgaCCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTc                                     <  2:4259115/149‑1 (MQ=255)
                                                                                                               tgaatgaCCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTc                                     >  2:1704441/1‑149 (MQ=255)
                                                                                                               tgaatgaCCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTc                                     >  1:2102158/1‑149 (MQ=255)
                                                                                                                  atgaCCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCggg                                  >  2:230350/1‑149 (MQ=255)
                                                                                                                      ccTGGTGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGt                              <  1:2880180/149‑1 (MQ=255)
                                                                                                                       cTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTc                             >  2:2951928/1‑149 (MQ=255)
                                                                                                                        tGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATTACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCg                            >  1:327648/1‑149 (MQ=255)
                                                                                                                         ggCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCAACCTACGGTCATCATGGCGGTGCCAGTCACCGAGCGCGGGAAATAATGCTGTTATAAGAAACAGCGGGTTAATACATAGTTATTTTACAAACGCTTTGCGTCGGGCACTCgg                           >  1:2674869/1‑149 (MQ=255)
                                                                                                                         ggCGCGTGCCGAAGCGGCGTTTATTAGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCgg                           >  1:4202181/1‑149 (MQ=255)
                                                                                                                          gcgcGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGa                          <  2:216869/149‑1 (MQ=255)
                                                                                                                                gCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGa                    >  2:773932/1‑149 (MQ=255)
                                                                                                                                  cGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAg                  <  1:1599785/149‑1 (MQ=255)
                                                                                                                                   gAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATTCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGt                 >  2:470220/1‑149 (MQ=255)
                                                                                                                                   gAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGt                 >  2:3384858/1‑149 (MQ=255)
                                                                                                                                     aGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCg               <  1:896618/149‑1 (MQ=255)
                                                                                                                                        ggcgTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGccc            <  2:5105942/149‑1 (MQ=255)
                                                                                                                                        ggcgTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGccc            >  2:1369393/1‑149 (MQ=255)
                                                                                                                                            tttATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGACTTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCAt        >  2:300722/1‑149 (MQ=255)
                                                                                                                                            tttATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCAt        <  1:4035320/149‑1 (MQ=255)
                                                                                                                                            tttATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCAt        >  1:4216572/1‑149 (MQ=255)
                                                                                                                                              tATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCAtgg      >  2:3231188/1‑149 (MQ=255)
                                                                                                                                               aTTAGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCAtggt     <  2:3079022/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGAGGCCTTTGGCATGTTGTTGATTCGAGAGGGATTAAAGGTCGCGGTTATTGCGGTCGATCCCAGCAGCCCGGTCACTGGCGGTAGCATTCTCGGGGATAAAACCCGCATGAATGACCTGGCGCGTGCCGAAGCGGCGTTTATTCGCCCGGTACCATCCTCCGGTCATCTGGGCGGTGCCAGTCAGCGAGCGCGGGAATTAATGCTGTTATGCGAAGCAGCGGGTTATGACGTAGTGATTGTCGAAACGGTTGGCGTCGGGCAGTCGGAAACAGAAGTCGCCCGCATGGTGGA  >  CP009273/3056557‑3056851

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: