Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 334,849 G→T 5.3% E272* (GAA→TAA)  yahG → DUF1116 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273334,8490GT 5.3% 1397.0 / 32.1 350E272* (GAA→TAA) yahGDUF1116 family protein
Reads supporting (aligned to +/- strand):  ref base G (167/164);  new base T (14/5);  total (181/169)
Fisher's exact test for biased strand distribution p-value = 5.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 363 total aligned reads displayed.

TCGGCGAAATCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGG  >  CP009273/334700‑334998
                                                                                                                                                     |                                                                                                                                                     
 cGGCGAAATCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGa                                                                                                                                                      >  2:1999528/1‑150 (MQ=255)
 cGGCGAAATCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGa                                                                                                                                                      >  2:3708557/1‑150 (MQ=255)
  ggCGAAATCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGaa                                                                                                                                                     <  1:2399802/150‑1 (MQ=255)
  ggCGAAATCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGaa                                                                                                                                                     >  2:156766/1‑150 (MQ=255)
     gAAATCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAgtg                                                                                                                                                  >  1:673520/1‑150 (MQ=255)
      aaaTCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAgtgt                                                                                                                                                 >  2:3596454/1‑150 (MQ=255)
      aaaTCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGGTATTCCGTCGAGCACCCGCGCTAAgtgt                                                                                                                                                 >  1:354909/1‑150 (MQ=255)
        aTCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGttt                                                                                                                                               <  2:334605/150‑1 (MQ=255)
        aTCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGttt                                                                                                                                               <  2:3902663/150‑1 (MQ=255)
         tCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTg                                                                                                                                              <  2:3380678/150‑1 (MQ=255)
           gATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAg                                                                                                                                            >  1:1431019/1‑150 (MQ=255)
           gATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAg                                                                                                                                            <  2:5005466/150‑1 (MQ=255)
           gATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAg                                                                                                                                            <  2:618108/150‑1 (MQ=255)
            aTCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGt                                                                                                                                           >  1:1830481/1‑150 (MQ=255)
              cTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGttt                                                                                                                                         >  1:2568781/1‑150 (MQ=255)
                gcgcTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATACCCGCAATACCGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGt                                                                                                                                       <  2:4229586/150‑1 (MQ=255)
                 cgcTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTc                                                                                                                                      <  1:3654492/150‑1 (MQ=255)
                   cTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGc                                                                                                                                    >  1:2332963/1‑150 (MQ=255)
                   cTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGc                                                                                                                                    >  1:1507040/1‑150 (MQ=255)
                   cTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAAACGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGc                                                                                                                                    >  1:4876713/1‑150 (MQ=255)
                    ttAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGcc                                                                                                                                   >  1:1075664/1‑150 (MQ=255)
                    ttAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGcc                                                                                                                                   >  1:4809368/1‑150 (MQ=255)
                     tAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCca                                                                                                                                  <  1:2663073/150‑1 (MQ=255)
                      aaTGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCcag                                                                                                                                 <  2:3005066/150‑1 (MQ=255)
                         gCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCcagcag                                                                                                                              <  1:4749999/150‑1 (MQ=255)
                          cTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCagcagc                                                                                                                             <  1:2007533/150‑1 (MQ=255)
                             gcgcAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGac                                                                                                                          <  2:2741067/150‑1 (MQ=255)
                              cgcAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGact                                                                                                                         <  1:253520/150‑1 (MQ=255)
                              cgcAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGact                                                                                                                         >  2:3813047/1‑150 (MQ=255)
                               gcAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGacta                                                                                                                        >  2:5137301/1‑150 (MQ=255)
                                    cgcTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTtct                                                                                                                   <  1:1131497/150‑1 (MQ=255)
                                     gcTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTtctc                                                                                                                  >  1:3704803/1‑150 (MQ=255)
                                     gcTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTtctc                                                                                                                  >  2:2293855/1‑150 (MQ=255)
                                     gcTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTtctc                                                                                                                  <  1:1254271/150‑1 (MQ=255)
                                     gcTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTtctc                                                                                                                  <  2:110798/150‑1 (MQ=255)
                                     gcTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTtctc                                                                                                                  <  1:4922177/150‑1 (MQ=255)
                                       tGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCg                                                                                                                >  2:659618/1‑150 (MQ=255)
                                       tGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCg                                                                                                                >  2:13105/1‑150 (MQ=255)
                                        gCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCgg                                                                                                               >  2:4586856/1‑150 (MQ=255)
                                        gCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCgg                                                                                                               >  2:4067531/1‑150 (MQ=255)
                                        gCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCgg                                                                                                               >  2:2962181/1‑150 (MQ=255)
                                           tatGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGccc                                                                                                            <  1:2928391/150‑1 (MQ=255)
                                           tatGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGccc                                                                                                            <  1:3120463/150‑1 (MQ=255)
                                           tatGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGccc                                                                                                            >  2:4827688/1‑150 (MQ=255)
                                            atGGGCGATGAGTGCCAGAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCNg                                                                                                           <  1:2821216/150‑1 (MQ=255)
                                              gggCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGAc                                                                                                         <  1:4382541/150‑1 (MQ=255)
                                                gCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGt                                                                                                       <  1:3029393/150‑1 (MQ=255)
                                                   aTGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGa                                                                                                    <  1:2251213/150‑1 (MQ=255)
                                                     gAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATg                                                                                                  >  2:1353527/1‑150 (MQ=255)
                                                             tAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGAtgtg                                                                                          <  2:497180/150‑1 (MQ=255)
                                                             tAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGAtgtg                                                                                          <  1:4654821/150‑1 (MQ=255)
                                                                ccGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTaa                                                                                       <  2:184913/150‑1 (MQ=255)
                                                                   cAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCTAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAgc                                                                                    >  1:646207/1‑150 (MQ=255)
                                                                   cAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAgc                                                                                    >  1:2855719/1‑150 (MQ=255)
                                                                       aaCGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAgcggcg                                                                                <  1:5132161/150‑1 (MQ=255)
                                                                         cGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCgat                                                                              >  1:602155/1‑150 (MQ=255)
                                                                         cGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCgat                                                                              >  1:4973724/1‑150 (MQ=255)
                                                                         cGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCgat                                                                              >  1:4194706/1‑150 (MQ=255)
                                                                         cGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCgat                                                                              >  2:43549/1‑150 (MQ=255)
                                                                            cGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCgatgga                                                                           >  2:1218612/1‑150 (MQ=255)
                                                                             gggACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCTAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCgatggat                                                                          >  1:1318038/1‑150 (MQ=255)
                                                                               gacgacACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATgc                                                                        <  2:143907/150‑1 (MQ=255)
                                                                                 cgacACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATgcgg                                                                      >  2:2381173/1‑150 (MQ=255)
                                                                                      ctgctgATTCAGGCGCTGCCGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCAt                                                                 <  1:4026589/150‑1 (MQ=255)
                                                                                       tgctgATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATg                                                                <  2:1686883/150‑1 (MQ=255)
                                                                                          tgATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCa                                                             >  2:4274941/1‑150 (MQ=255)
                                                                                                   cgcTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACa                                                    <  1:3307159/150‑1 (MQ=255)
                                                                                                      tGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCa                                                 >  2:1159318/1‑150 (MQ=255)
                                                                                                        aCGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCAcc                                               >  1:3922218/1‑150 (MQ=255)
                                                                                                           ccGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTg                                            >  1:991660/1‑150 (MQ=255)
                                                                                                            cGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTgg                                           >  2:2448264/1‑150 (MQ=255)
                                                                                                            cGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTgg                                           >  1:1849491/1‑150 (MQ=255)
                                                                                                             ggggATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGt                                          >  2:16243/1‑150 (MQ=255)
                                                                                                               ggATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTca                                        <  1:682738/150‑1 (MQ=255)
                                                                                                                gATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTcac                                       <  1:4498780/150‑1 (MQ=255)
                                                                                                                    attCAGGCGGGTTATTCCGTCGAGCAACAGCGCTAAGTGTTTGAGTTTGTGGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCaccacc                                   >  2:2543154/1‑150 (MQ=255)
                                                                                                                    attCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCaccacc                                   <  2:805474/150‑1 (MQ=255)
                                                                                                                    attCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCaccacc                                   >  1:3212314/1‑150 (MQ=255)
                                                                                                                     ttCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCAccacca                                  <  2:1887050/150‑1 (MQ=255)
                                                                                                                      tCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCAt                                 <  1:1875201/150‑1 (MQ=255)
                                                                                                                       cAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATg                                <  2:2632581/150‑1 (MQ=255)
                                                                                                                       cAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATg                                <  1:1382703/150‑1 (MQ=255)
                                                                                                                        aGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATgg                               <  2:4710463/150‑1 (MQ=255)
                                                                                                                             ggTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGt                          <  1:16403/150‑1 (MQ=255)
                                                                                                                             ggTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGt                          >  1:92915/1‑150 (MQ=255)
                                                                                                                             ggTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGt                          >  2:4890571/1‑150 (MQ=255)
                                                                                                                                tATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAAc                       >  2:344258/1‑150 (MQ=255)
                                                                                                                                tATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAAc                       >  1:43198/1‑150 (MQ=255)
                                                                                                                                 aTTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAAca                      <  2:871794/150‑2 (MQ=255)
                                                                                                                                  ttCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACgg                     <  2:75388/150‑1 (MQ=255)
                                                                                                                                     cGTCGAGCAACAGCGCGAAGTGTTGGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGt                  <  2:2226775/150‑1 (MQ=255)
                                                                                                                                      gTCGAGCAACAGCGCTAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTc                 >  1:1471801/1‑150 (MQ=255)
                                                                                                                                        cGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGa               >  2:4303798/1‑150 (MQ=255)
                                                                                                                                         gAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAg              >  2:321855/1‑150 (MQ=255)
                                                                                                                                         gAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAg              >  1:2361818/1‑150 (MQ=255)
                                                                                                                                             aaCAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCg          >  1:797434/1‑150 (MQ=255)
                                                                                                                                              aCAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCgg         >  1:5048570/1‑150 (MQ=255)
                                                                                                                                                 gcgcGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCt      <  2:3285102/150‑1 (MQ=255)
                                                                                                                                                   gcgAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGc    >  1:2549554/1‑150 (MQ=255)
                                                                                                                                                     gAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCgg  >  1:2209538/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGCGAAATCGATCTGCGCTTAATGCTGGCGCAGGCGCTGCATATGGGCGATGAGTGCCATAACCGCAATAACGCCGGGACGACACTGCTGATTCAGGCGCTGACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGG  >  CP009273/334700‑334998

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: