Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,139,910 G→T 92.3% R1055L (CGA→CTA)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,139,9100GT92.3% 870.8 / 62.0 297R1055L (CGA→CTA) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (4/19);  new base T (149/125);  total (153/144)
Fisher's exact test for biased strand distribution p-value = 7.81e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 298 total aligned reads displayed.

CAGCCAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCGATT  >  CP009273/2139761‑2140058
                                                                                                                                                     |                                                                                                                                                    
   ccAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACt                                                                                                                                                   >  2:3324938/1‑150 (MQ=255)
     aGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCg                                                                                                                                                 >  1:35792/1‑150 (MQ=255)
     aGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCg                                                                                                                                                 >  2:2746581/1‑150 (MQ=255)
      ggTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCgg                                                                                                                                                <  2:939734/150‑1 (MQ=255)
      ggTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCgg                                                                                                                                                <  2:3782008/150‑1 (MQ=255)
         gggCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGAt                                                                                                                                             <  2:2690856/150‑1 (MQ=255)
         gggCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGAt                                                                                                                                             <  2:2690847/150‑1 (MQ=255)
         gggCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGAt                                                                                                                                             >  2:296664/1‑150 (MQ=255)
          ggCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATg                                                                                                                                            >  2:2882657/1‑150 (MQ=255)
                gATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAAcc                                                                                                                                      >  1:1635138/1‑150 (MQ=255)
                    tGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCg                                                                                                                                  >  1:2089219/1‑150 (MQ=255)
                    tGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCg                                                                                                                                  <  1:2477662/150‑1 (MQ=255)
                       aaaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCg                                                                                                                               >  2:4136206/1‑150 (MQ=255)
                       aaaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCg                                                                                                                               <  1:115327/150‑1 (MQ=255)
                         aTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCggg                                                                                                                             <  1:3604814/150‑1 (MQ=255)
                           cTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGcc                                                                                                                           >  1:3890305/1‑150 (MQ=255)
                            ttCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCa                                                                                                                          <  2:2497260/150‑1 (MQ=255)
                               aaCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTc                                                                                                                       <  1:548417/150‑1 (MQ=255)
                                 cTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGt                                                                                                                     >  1:2987628/1‑150 (MQ=255)
                                 cTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGt                                                                                                                     <  2:1743039/150‑1 (MQ=255)
                                  tCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtt                                                                                                                    <  2:283956/150‑1 (MQ=255)
                                    gCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtttt                                                                                                                  >  2:2344114/1‑150 (MQ=255)
                                      tGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAc                                                                                                                >  1:3081836/1‑150 (MQ=255)
                                      tGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAc                                                                                                                >  2:2661769/1‑150 (MQ=255)
                                       gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAcc                                                                                                               >  1:568385/1‑150 (MQ=255)
                                       gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAcc                                                                                                               >  1:1932441/1‑150 (MQ=255)
                                       gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAcc                                                                                                               >  1:498411/1‑150 (MQ=255)
                                       gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAcc                                                                                                               >  1:498262/1‑150 (MQ=255)
                                        aaaGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAccc                                                                                                              <  1:327936/150‑1 (MQ=255)
                                         aaGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCCGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCt                                                                                                             >  2:4490869/1‑150 (MQ=255)
                                          aGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCtt                                                                                                            >  1:1106147/1‑150 (MQ=255)
                                            cGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAg                                                                                                          >  2:3226314/1‑150 (MQ=255)
                                               atatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGa                                                                                                       <  1:4215196/150‑1 (MQ=255)
                                                tatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGAt                                                                                                      <  1:763550/150‑1 (MQ=255)
                                                 atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTCATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATg                                                                                                     >  2:4424893/1‑150 (MQ=255)
                                                 atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATg                                                                                                     >  2:1350276/1‑150 (MQ=255)
                                                 atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATg                                                                                                     >  2:3617560/1‑150 (MQ=255)
                                                      aGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAc                                                                                                <  1:814005/150‑1 (MQ=255)
                                                           aCCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCttt                                                                                           <  1:4401875/150‑1 (MQ=255)
                                                           aCCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCttt                                                                                           >  1:3545293/1‑150 (MQ=255)
                                                             cTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCt                                                                                         <  2:2034281/150‑1 (MQ=255)
                                                              tGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTg                                                                                        >  2:4304827/1‑150 (MQ=255)
                                                               gCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTgg                                                                                       <  1:1621627/150‑1 (MQ=255)
                                                                cTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGc                                                                                      >  1:1091119/1‑150 (MQ=255)
                                                                 tACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCa                                                                                     >  2:3226634/1‑150 (MQ=255)
                                                                 tACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCa                                                                                     >  2:1371103/1‑150 (MQ=255)
                                                                   cTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATc                                                                                   >  1:2420249/1‑150 (MQ=255)
                                                                      gATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGc                                                                                >  2:763152/1‑150 (MQ=255)
                                                                      gATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGc                                                                                <  1:4114787/150‑1 (MQ=255)
                                                                      gATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGc                                                                                <  1:720789/150‑1 (MQ=255)
                                                                         ggTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTc                                                                             >  2:1623134/1‑150 (MQ=255)
                                                                            gAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGAt                                                                          <  2:1332101/150‑1 (MQ=255)
                                                                            gAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGAt                                                                          <  1:758650/150‑1 (MQ=255)
                                                                            gAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGAt                                                                          <  2:1541469/150‑1 (MQ=255)
                                                                              gTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGACGATCt                                                                        >  2:3764668/1‑150 (MQ=255)
                                                                                 aTGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAt                                                                     <  1:1749875/150‑1 (MQ=255)
                                                                                   gcgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAttt                                                                   <  1:3324938/150‑1 (MQ=255)
                                                                                    cgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTa                                                                  >  2:3431431/1‑150 (MQ=255)
                                                                                    cgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTa                                                                  <  1:4327776/150‑1 (MQ=255)
                                                                                     gcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTAt                                                                 <  2:172750/150‑1 (MQ=255)
                                                                                     gcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTAt                                                                 <  2:1034049/150‑1 (MQ=255)
                                                                                     gcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTAt                                                                 >  2:2625123/1‑150 (MQ=255)
                                                                                        aaCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggt                                                              >  1:16719/1‑150 (MQ=255)
                                                                                           gTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgag                                                           >  1:2738427/1‑150 (MQ=255)
                                                                                             gCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggt                                                         <  1:934205/150‑1 (MQ=255)
                                                                                              cAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggtg                                                        >  1:3938089/1‑150 (MQ=255)
                                                                                                     aTTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCg                                                 >  1:3218881/1‑150 (MQ=255)
                                                                                                       ttGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGAt                                               >  1:2061557/1‑150 (MQ=255)
                                                                                                       ttGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGAt                                               <  1:4087098/150‑1 (MQ=255)
                                                                                                        tGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATg                                              <  2:3851046/150‑1 (MQ=255)
                                                                                                            ggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCa                                          <  2:1008914/150‑1 (MQ=255)
                                                                                                            ggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCa                                          >  1:3741885/1‑150 (MQ=255)
                                                                                                             gatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCaa                                         >  2:4275742/1‑150 (MQ=255)
                                                                                                              atgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAAc                                        >  1:4215664/1‑150 (MQ=255)
                                                                                                              atgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAAc                                        <  2:593624/150‑1 (MQ=255)
                                                                                                                      cTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAt                                <  2:3814725/150‑1 (MQ=255)
                                                                                                                       tGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAtt                               <  2:3078873/150‑1 (MQ=255)
                                                                                                                       tGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAtt                               <  2:3077756/150‑1 (MQ=255)
                                                                                                                        gATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAttt                              >  2:1227811/1‑150 (MQ=255)
                                                                                                                        gATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAttt                              >  2:736399/1‑150 (MQ=255)
                                                                                                                        gATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAttt                              >  2:2817216/1‑150 (MQ=255)
                                                                                                                         aTTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTg                             <  2:2440800/150‑1 (MQ=255)
                                                                                                                           tACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCt                           <  2:4201802/150‑1 (MQ=255)
                                                                                                                            aCGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTg                          <  2:1852363/150‑1 (MQ=255)
                                                                                                                              gATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGCTCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGt                        >  2:3987978/1‑150 (MQ=255)
                                                                                                                                ttattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGa                      <  1:2848624/150‑1 (MQ=255)
                                                                                                                                  attCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGTAt                    >  1:3425963/1‑150 (MQ=255)
                                                                                                                                     cAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTGTCTGGCATCGGCGTCGATCTGAGTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagt                 >  2:1249521/1‑150 (MQ=255)
                                                                                                                                      aGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagta                <  1:296488/150‑1 (MQ=255)
                                                                                                                                       ggggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagtag               >  2:3152300/1‑150 (MQ=255)
                                                                                                                                       ggggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagtag               <  1:2356778/150‑1 (MQ=255)
                                                                                                                                       ggggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagtag               >  2:175017/1‑150 (MQ=255)
                                                                                                                                        gggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATagtagt              <  2:3270718/150‑1 (MQ=255)
                                                                                                                                          gCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTa            <  2:4102231/150‑1 (MQ=255)
                                                                                                                                          gCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTa            >  1:1483364/1‑150 (MQ=255)
                                                                                                                                          gCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTa            <  1:222659/150‑1 (MQ=255)
                                                                                                                                          gCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTa            >  1:3769774/1‑150 (MQ=255)
                                                                                                                                               cTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTcg       <  2:4044236/150‑1 (MQ=255)
                                                                                                                                                 cAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTcgcg     >  1:4193629/1‑150 (MQ=255)
                                                                                                                                                    cGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCGAtt  <  2:3465803/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGCCAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCGATT  >  CP009273/2139761‑2140058

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: