Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,755,723 G→T 78.3% E61* (GAA→TAA)  rhsA → Rhs family protein, putative polymorphic toxin; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,755,7230GT78.3% 119.7 / 61.7 97E61* (GAA→TAA) rhsARhs family protein, putative polymorphic toxin; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base T (27/49);  total (39/58)
Fisher's exact test for biased strand distribution p-value = 8.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGC  >  CP009273/3755575‑3755770
                                                                                                                                                    |                                               
cATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGa                                                >  1:2378734/1‑150 (MQ=34)
 aTGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaa                                               <  2:1206932/150‑1 (MQ=21)
  tGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaaa                                              >  1:819569/1‑150 (MQ=21)
  tGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaaa                                              >  2:3885707/1‑150 (MQ=21)
  tGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaaa                                              >  2:3273509/1‑150 (MQ=21)
   gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGTCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAAc                                             >  2:3620225/1‑150 (MQ=21)
   gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAAc                                             <  1:2896416/150‑1 (MQ=21)
   gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAAc                                             <  1:2735426/150‑1 (MQ=21)
   gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCAGTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAAc                                             <  1:4753452/150‑1 (MQ=21)
    aCGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAAcc                                            <  1:578727/150‑1 (MQ=21)
    aCGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAAcc                                            <  2:1231709/150‑1 (MQ=21)
    aCGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAcc                                            >  1:4051930/1‑150 (MQ=34)
    aCGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCAACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAAcc                                            <  2:2669156/150‑1 (MQ=21)
    aCGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCAACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAAcc                                            <  2:2140627/150‑1 (MQ=21)
    aCGCAGTATGGCGGGAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAcc                                            <  2:4500719/150‑1 (MQ=21)
     cGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCg                                           <  2:1350881/150‑1 (MQ=21)
       cAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTACAAAGGTCCTTCCCGGTGAAACCGAc                                         >  1:3819999/1‑150 (MQ=21)
         gTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACAt                                       <  2:1533097/150‑1 (MQ=21)
         gTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACAt                                       <  2:149093/150‑1 (MQ=21)
         gTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACAt                                       <  2:5166575/150‑1 (MQ=21)
         gTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACAt                                       <  1:2153767/150‑1 (MQ=21)
         gTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACAt                                       <  2:133858/150‑1 (MQ=21)
         gTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACAt                                       <  1:565859/150‑1 (MQ=21)
         gTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACAt                                       >  2:3869543/1‑150 (MQ=34)
          tATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATc                                      <  2:4344183/150‑1 (MQ=21)
          tATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATc                                      <  1:3474402/150‑1 (MQ=21)
          tATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATc                                      <  1:3474512/150‑1 (MQ=21)
          tATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATc                                      >  2:4751969/1‑150 (MQ=21)
          tATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATc                                      <  2:4344432/150‑1 (MQ=21)
          tATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATc                                      <  2:2179604/150‑1 (MQ=21)
           aTGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCg                                     <  2:2401883/150‑1 (MQ=21)
           aTGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCg                                     <  1:4465506/150‑1 (MQ=34)
            tGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGc                                    >  2:1739592/1‑150 (MQ=21)
               cGGTAGCATTGTTCAGGGTTCTGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCt                                 <  1:642011/150‑1 (MQ=21)
               cGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTATGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCt                                 <  1:4406781/150‑1 (MQ=21)
                ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTg                                >  1:2376496/1‑150 (MQ=21)
                ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTg                                >  1:4146608/1‑150 (MQ=21)
                ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTg                                >  1:2376551/1‑150 (MQ=21)
                ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTg                                <  2:4882497/150‑1 (MQ=34)
                ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTg                                >  1:4518706/1‑150 (MQ=34)
                 gTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgc                               >  1:4769391/1‑150 (MQ=34)
                  tAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcc                              >  1:1963720/1‑150 (MQ=21)
                  tAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcc                              >  1:1963817/1‑150 (MQ=21)
                   aGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgccc                             >  1:3296925/1‑150 (MQ=21)
                   aGCATTGGTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgccc                             >  2:2352048/1‑150 (MQ=21)
                    gCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccg                            <  1:1223480/150‑1 (MQ=21)
                     cATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccgg                           >  1:3569497/1‑150 (MQ=21)
                      aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggc                          >  1:4820702/1‑150 (MQ=21)
                      aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggc                          >  1:3860798/1‑150 (MQ=21)
                      aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggc                          <  2:3122141/150‑1 (MQ=34)
                      aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggc                          >  2:3195118/1‑150 (MQ=34)
                      aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggc                          >  2:2893299/1‑150 (MQ=34)
                      aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggc                          >  1:172428/1‑150 (MQ=34)
                        tGTTCAGGGTTCAGCCGGGGTGCTCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggccc                        <  1:2704076/150‑1 (MQ=21)
                        tGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggccc                        <  2:1698263/150‑1 (MQ=21)
                        tGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggccc                        <  1:2554790/150‑1 (MQ=21)
                        tGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggccc                        >  2:5291331/1‑150 (MQ=34)
                         gTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggcccg                       <  2:4718089/150‑1 (MQ=34)
                           tCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCt                     <  2:3901567/150‑1 (MQ=21)
                           tCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCt                     <  1:4585423/150‑1 (MQ=34)
                            cAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTg                    <  1:2474189/150‑1 (MQ=21)
                             aGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGc                   >  1:2351464/1‑150 (MQ=21)
                             aGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGc                   >  1:1768803/1‑150 (MQ=21)
                             aGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTACGTGTGCCCCGGCGGAGTGACGTCCGGCCAACCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCTCCCTGCCCGGCCCGCTGc                   >  1:2351484/1‑150 (MQ=11)
                              gggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGcc                  <  2:2174823/150‑1 (MQ=21)
                              gggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGcc                  >  1:2820784/1‑150 (MQ=21)
                              gggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGcc                  >  1:2820785/1‑150 (MQ=21)
                              gggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGGTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCAACCCGTCAATCCCCAGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACCTCGCCCTCCCCGGCCCGCTGcc                  >  1:1910613/1‑150 (MQ=11)
                               ggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCg                 <  1:4169621/150‑1 (MQ=21)
                               ggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGACAAACCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCg                 >  2:1260559/1‑150 (MQ=255)
                                 ttCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGtt               >  1:1877386/1‑150 (MQ=34)
                                  tCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTc              <  1:4296090/150‑1 (MQ=21)
                                  tCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTc              <  1:1592718/150‑1 (MQ=21)
                                  tCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGGTc              <  1:4296184/150‑1 (MQ=21)
                                  tCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTc              <  2:1949954/150‑1 (MQ=34)
                                   cAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCa             <  2:844571/150‑1 (MQ=21)
                                   cAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCa             <  2:186443/150‑1 (MQ=21)
                                      ccGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATcc          <  1:35223/150‑1 (MQ=21)
                                       cGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCct         <  2:2094948/150‑1 (MQ=21)
                                       cGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCct         <  1:929869/150‑1 (MQ=21)
                                        ggggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCGGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctc        >  2:363473/1‑150 (MQ=21)
                                        ggggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctc        <  2:409022/150‑1 (MQ=21)
                                        ggggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctc        <  1:557732/150‑1 (MQ=21)
                                        ggggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctc        <  1:557891/150‑1 (MQ=21)
                                        ggggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctc        >  2:1181029/1‑150 (MQ=21)
                                        ggggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctc        >  2:3857820/1‑150 (MQ=21)
                                        ggggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGGCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCCGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATccac        >  1:5402/1‑148 (MQ=255)
                                         gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctct       <  2:4089509/150‑1 (MQ=21)
                                         gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctct       <  2:4171692/150‑1 (MQ=21)
                                         gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctct       <  1:544568/150‑1 (MQ=21)
                                           gTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTcc     <  2:1391070/150‑1 (MQ=21)
                                            tGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTccc    <  1:790744/150‑1 (MQ=18)
                                            tGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTccc    <  1:2264199/150‑1 (MQ=18)
                                            tGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTccc    <  1:1129789/150‑1 (MQ=18)
                                            tGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTccc    <  1:1076719/150‑1 (MQ=18)
                                             gcgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCg   <  1:3047584/150‑1 (MQ=17)
                                              cgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGc  >  2:119576/1‑150 (MQ=17)
                                                                                                                                                    |                                               
CATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGC  >  CP009273/3755575‑3755770

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: